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The canonical example of a quantum mechanical two-level system is spin. The simplest picture 
of spin is a magnetic moment pointing up or down. The full quantum properties of spin become 
apparent in phenomena such as superpositions of spin states, entanglement among spins and 
quantum measurements. Many of these phenomena have been observed in experiments performed 
on ensembles of particles with spin. Only in recent years systems have been realized in which 
individual electrons can be trapped and their quantum properties can be studied, thus avoiding 
unnecessary ensemble averaging. This review describes experiments performed with quantum dots, 
which are nanometer-scale boxes defined in a semiconductor host material. Quantum dots can hold 
a precise, but tunable number of electron spins starting with 0, 1, 2, etc. Electrical contacts can be 
made for charge transport measurements and electrostatic gates can be used for controlling the dot 
potential. This system provides virtually full control over individual electrons. This new, enabling 
technology is stimulating research on individual spins. This review describes the physics of spins 
in quantum dots containing one or two electrons, from an experimentalist's viewpoint. Various 
methods for extracting spin properties from experiment are presented, restricted exclusively to 
electrical measurements. Furthermore, experimental techniques are discussed that allow for: (1) 
the rotation of an electron spin into a superposition of up and down, (2) the measurement of the 
quantum state of an individual spin and (3) the control of the interaction between two neighbouring 
spins by the Heisenberg exchange interaction. Finally, the physics of the relevant relaxation and 
dephasing mechanisms is reviewed and experimental results are compared with theories for spin- 
orbit and hyperfine interactions. All these subjects are directly relevant for the fields of quantum 
information processing and spintronics with single spins (i.e. single-spintronics). 
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I. INTRODUCTION 

The spin of an electron remains a somewhat mysteri- 
ous property. The first derivations in 1925 of the spin 
magnetic moment, based on a rotating charge distribu- 
tion of finite size, are in conflict with special relativity 
theory. Pauli advised the young Ralph Kronig not to 
publish his theory since "it has nothing to do with real- 
ity" . More fortunate were Samuel Goudsmit and George 
Uhlcnbcck, who were supervised by Ehrcnfcst: "Publish, 
you are both young enough to be able to afford a stupid- 
ity!" 1 . It requires Dirac's equation to find that the spin 
eigenvalues correspond to one-half times Planck's con- 
stant, h, while considering the electron as a point parti- 
cle. The magnetic moment corresponding to spin is really 
very small and in most practical cases it can be ignored. 
For instance, the most sensitive force sensor to date has 
only recently been able to detect some e ffect from the 
magnetic moment of a single electron spin ( Rugar et al\ . 



I2004D . In solids, spin can apparently lead to strong ef- 
fects, given the existence of permanent magnets. Curi- 
ously, this has little to do with the strength of the mag- 
netic moment. Instead, the fact that spin is associated 
with its own quantum number, combined with Pauli's 
exclusion principle that quantum states can at most be 
occupied with one fermion, leads to the phenomenon of 
exchange interaction. Because the exchange interaction 
is a correction term to the strong Coulomb interaction, it 
can be of much larger strength in solids than the dipolar 
interaction between two spin magnetic moments at an 
atomic distance of a few Angstroms. It is the exchange 
interaction that forces the electron spins in a collective 
align ment, together yielding a ma croscopic magnetiza- 
tion ( Ashcroft and Merminl IT974 ). It remains striking, 



1 See http:/ /www. lorentz.leidenuniv.nl/history/spin/goudsmit. html 



that an abstract concept as (anti-)symmetrization in the 
end gives rise to magnets. 

The magnetic state of solids has found important appli- 
cations in electronics, in particular for memory devices. 
An important field has emerged in the last two decades 
known as spintronics. Phenomena like Giant Magneto 
Resistance or Tunneling Magneto Resistance form the 
basis for magnetic heads for reading out the magnetic 
state of a memory cell. Logic gates have b een realized 
based on magnetoresist ance effects as well ( Wolf et~al\ , 
l200ll: IZut ic et al. 1. 12004 . In addition to applications, im- 
portant scienti fic discoveries have been made in the field 
of spintronics ( Awschalom and Flattd . l2007f ). including 
magne tic semiconducto rs ( Ohnol . 1998f ) and the spin Hall 
effect ( Sih et al\ , 120051 ). It is important to note that all 
the spintronics phenomena consider macroscopic num- 
bers of spins. Together these spins form things like spin 
densities or a collective magnetization. Although the 
origin of spin densities and magnetization is quantum 
mechanical, these collective, macroscopic variables be- 
have entirely classically. For instance, the magnetization 
of a micron-cubed piece of Cobalt is a classical vector. 
The quantum state of this vector dephases so rapidly 
that quantum superpositions or entanglement between 
vectors is never observed. One has to go to systems 
with a small number of spins, for instance in magnetic 
molecules, in order to find quantum effects in the be- 
haviour of the collective magnetizatio n (for an overview, 
see e.g.. lGunther and Barbara! (|l994l )). 

The technological drive to make electronic devices con- 
tinuously smaller has some interesting scientific conse- 
quences. For instance, it is now routinely possible to 
make small electron "boxes" in solid state devices that 
contain an integer number of conduction electrons. Such 
devices are usually operated as transistors (via field-effect 
gates) and are therefore named single electron transis- 
tors. In semiconductor boxes the number of trapped 
electrons can be reduced all the way to zero, or one, two, 
etc. Such semiconduc tor single electron t r ansist ors are 
called quantum dots ( Kouwenhoven et "all . l200l[ ). Elec- 
trons are trapped in a quantum dot by repelling electric 
fields imposed from all sides. The final region in which 
a small number of electrons can still exist is typically at 
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the scale of tens of nanometers. The eigenenergies in such 
boxes are discrete. Filling these states with electrons fol- 
lows the rules from atomic physics, including Hund's rule, 
shell filling, etc. 

Studies with quantum dots have been performed suc- 
cessfully during the nineties. By now it has become stan- 
dard technology to confine single electron charges. Elec- 
trons can be trapped as long as one desires. Changes in 
charge when one electron tunnels out of the quantum dot 
can be measured on a microsecond timescale. Compared 
to this control of charge, it is very difficult to control 
individual spins and measure the spin of an individual 
electron. Such techniques have been developed only over 
the past few years. 

In this review we describe experiments in which indi- 
vidual spins are controlled and measured. This is mostly 
an experimental review with explanations of the under- 
lying physics. This review is strictly limited to experi- 
ments that involve one or two electrons strongly confined 
to single or double quantum dot devices. The experi- 
ments show that one or two electrons can be trapped in 
a quantum dot; that the spin of an individual electron 
can be put in a superposition of up and down states; 
that two spins can be made to interact and form an en- 
tangled state such as a spin singlet or triplet state; and 
that the result of such manipulation can be measured on 
individual spins. 

These abilities of almost full control over the spin of 
individual electrons enable the investigation of a new 
regime: single spin dynamics in a solid state environ- 
ment. The dynamics are fully quantum mechanical and 
thus quantum coherence can be studied on an individ- 
ual electron spin. The exchange interaction is now also 
controlled on the level of two particular spins that are 
brought into contact simply by varying some voltage 
knob. 

In a solid the electron spins are not completely decou- 
pled from other degrees of freedom. First of all, spins and 
orbits are coupled by the spin-orbit interaction. Second, 
the electron spins have an interaction with the spins of 
the atomic nuclei, i.e. the hyperfine interaction. Both 
interactions cause the life time of a quantum superposi- 
tion of spin states to be finite. We therefore also describe 
experiments that probe spin-orbit and hyperfine interac- 
tions by measuring the dynamics of individual spins. 

The study of individual spins is motivated by an in- 
terest in fundamental physics, but also by possible ap- 
plications. First of all, miniaturized spintronics is de- 
veloping towards single spins. In this context, this 
field can be denoted as single-spintronics 2 in analogy 
to single-electronics. A second area of applications is 
quantum information science. Here the spin states form 
the q ubits. The original proposa l by Loss and DiVin- 
cenzo ( Loss and DiVincenzol . Il998l) has been the guide in 



this field. In the context of quantum information, the ex- 
periments described in this review demonstrate that the 
five DiVincenzo criteria for universal quantum computa- 
tion using si ngle electron spins have been fulfilled to a 
large extent ( DiVincenzol . |2000| ): initialization, one- and 
two-qubit operations, long coherence times and readout. 
Currently, the state of the art is at the level of single and 
double quantum dots and much work is required to build 
larger systems. 

In this review the system of choice are quantum dots 
in GaAs semiconductors, simply because these have been 
most successful. Nevertheless, the physics is entirely gen- 
eral and can be fully applied to new material systems 
such as silicon based transistors, carbon nanotubes, semi- 
conductor nanowires, graphene devices, etc. These other 
host materials may have advantageous spin properties. 
For instance, carbon-based devices can be purified with 
the isotope 12 C in which the nuclear spin is zero, thus 
entirely suppressing spin dephasing by hyperfine interac- 
tion. This kind of hardware solution to engineer a long- 
lived quantum system will be discussed at the end of this 
review. Also, we here restrict ourselves exclusively to 
electron transport measurements of quantum dots, leav- 
ing out optical spectroscopy of quantum dots, which is a 
very active field in its own 3 . Again, much of the physics 
discussed in this review also applies to optically measured 
quantum dots. 

Section[TT]starts with an introduction on quantum dots 
including the basic model of Coulomb blockade to de- 
scribe the relevant energies. These energies can be visu- 
alized in transport experiments and the relation between 
experimental spectroscopic lines and underlying energies 
are explained in section IIIII This spectroscopy is specif- 
ically applied to spin states in single quantum dots in 
section IIVI Section [V] introduces a charge-sensing tech- 
nique that is used in section I VII to read out the spin 
state of individual electrons. Section IVIII provides an 
extensive description of spin-orbit and hyperfine interac- 
tions. In section IVIII1 spin states in double quantum 
dots are introduced and the important concept of Pauli 
spin blockade is discussed. Quantum coherent manipula- 
tions of spins in double dots are discussed in section HXl 
Finally, a perspective is outlined in section X. 



II. BASICS OF QUANTUM DOTS 

A. Introduction to quantum dots 

A quantum dot is an artificially structured system that 
can be filled with electrons (or holes). The dot can be 
coupled via tunnel barriers to reservoirs, with which elec- 
trons can be exchanged (see Fig.rjJ. By attaching current 



2 Name coined by Stu Wolf, private communication. 



3 see e.g. Atature et al\ teOOrjt); iBerezovskv et al. I 1 120061) : 

IGreilich et all l|2006ar i: I Krenner et al .1 l|2006ft and references 
therein 
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and voltage probes to these reservoirs, we can measure 
the electronic properties. The dot is also coupled capaci- 
tively to one or more 'gate' electrodes, which can be used 
to tune the electrostatic potential of the dot with respect 
to the reservoirs. 




FIG. 1 Schematic picture of a quantum dot in (a) a lat- 
eral geometry and (b) in a vertical geometry. The quantum 
dot (represented by a disk) is connected to source and drain 
reservoirs via tunnel barriers, allowing the current through 
the device, I, to be measured in response to a bias voltage, 
Vsd and a gate voltage, Vg- 



Because a quantum dot is such a general kind 
of system, there exist quantum dots of many differ- 
ent sizes and materials: for instance sin g le mo lecules 
trapped be tween electrodes ([Park et all 120021) . nor- 
mal metaldPetta and Ralph! 12001 ^ su percon ducting 



(jvon Delft and Ralph! . 2001; Ralph et al 




romagnetic nanoparticles (IGueron et all 

assembled quan tum dots (IKlein et alT 19961) . semicon- 
ducto r lateral feouwenhoven et all Il997t) or vertical 
dots ( Kouwenhoven et all 200lt ). and also semiconduct- 
ing nanowire s or carbon nan otubes ( Biork et all 12004 
iDekkerl . I1999L iMcEuenl . l2000h . 

The electronic properties of quantum dots are domi- 
nated by two effects. First, the Coulomb repulsion be- 
tween the electrons on the dot leads to an energy cost 
for adding an extra electron to the dot. Due to this 
charging energy, tunneling of electrons to or from the 
reservoirs can be dramatically suppressed at low tem- 
peratures; this phenome na is called Coulomb blockade 
van Houten et a7] . ll992h . Second, the confinement in all 



three directions leads to quantum effects that strongly 
influence the electron dynamics. Due to the resulting 
discrete energy spectrum , quantum dots behave in m any 
ways as artificial atoms ( Kouwenhoven et all l200lh . 



The physics of dots contai ning more than two el e ctrons 
has been reviewed before (|Kouwenhoven et al. |, 119971: 
iReimann and Manninenl . I2002D . Therefore, we focus 
on single and coupled quantum dots containing only 
one or two electrons. These systems are particularly 
important as they constitute the building blocks of 



cessors dBvrd and Lidai 


, 20021 iDiVincenzo et all 2000 


Hanson and Burkard , 


20071 Kvriakidis and Pennev 


2005|: Lew 


, 20021 Loss and DiVincenzo. 1998 


Meier et all 


2003: ITavlor et al, 120051: IWu and Lidar 


2OO2a0). 



B. Fabrication of gated quantum dots 

The bulk of the experiments discussed in this review 
was performed on electrostatically defined quantum dots 
in GaAs. These devices are sometimes referred to as 
"lateral dots" because of the lateral gate geometry. 

Lateral GaAs quantum dots are fabricated from het- 
erostructures of GaAs and AlGaAs grown by molecular 
beam epitaxy, (see Fig. [5]) . By doping the AlGaAs layer 
with Si, free electrons are introduced. These accumulate 
at the GaAs/AlGaAs interface, typically 50-100 nm be- 
low the surface, forming a two-dimensional electron gas 
(2DEG) - a thin (^10 nm) sheet of electrons that can 
only move along the interface. The 2DEG can have a 
high mobility and relatively low electron density (typ- 
ically 10 5 - 10 7 cm 2 /Vs and - 1 - 5 x 10 15 m" 2 , re- 
spectively). The low electron density results in a large 
Fermi wavelength (~ 40 nm) and a large screening length, 
which allows us to locally deplete the 2DEG with an 
electric field. This electric field is created by applying 
negative voltages to metal gate electrodes on top of the 
heterostructure (see Fig. [2^) . 

Electron-beam lithography enables fabrication of gate 
structures with dimensions down to a few tens of nanome- 
ters (Fig. [2]), yielding local control over the depletion 
of the 2DEG with roughly the same spatial resolution. 
Small islands of electrons can be isolated from the rest of 
the 2DEG by choosing a suitable design of the gate struc- 
ture, thus creating quantum dots. Finally, low-resistance 
(Ohmic) contacts are made to the 2DEG reservoirs. To 
access the quantum phenomena in GaAs gated quantum 
dots, they have to be cooled down to well below 1 K. 
All experiments that are discussed in this review are per- 
formed in dilution refrigerators with typical base temper- 
atures of 20 mK. 

In so-called vertical quantum dots, control over the 
number of el ectrons down to zero was already achieved 
in the 1990s ( Kouwenhoven et "oLl. l200lh . In lateral gated 
dots this proved to be more difficult, since reducing the 
electron number by driving the gate voltage to more neg- 
ative values tends to decrease the tunnel coupling to the 
leads. The resulting current through the dot can then be- 
come unmeasurably small before the few-electron regime 
is reached. However, by proper design of the surface gate 
geometry the decrease of the tunnel coupling can be com- 
pensated for. 

In 2000, Ciorga et al. reported mea surements on 
the fi rst lateral few-electron quantum dot ( Ciorga et all 
I2000D . Their device, shown in Fig. f2J}, makes use of two 
types of gates specifically designed to have different func- 
tionalities. The gates of one type are big and largely 
enclose the quantum dot. The voltages on these gates 
mainly determine the dot potential. The other type of 
gate is thin and just reaches up to the barrier region. 
The voltage on this gate has a very small effect on the 
dot potential but it can be used to set the tunnel barrier. 
The combination of the two gate types allows the dot po- 
tential (and thereby electron number) to be changed over 
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FIG. 2 Lateral quantum dot device denned by metal surface 
electrodes, (a) Schematic view. Negative voltages applied 
to metal gate electrodes (dark gray) lead to depleted regions 
(white) in the 2DEG (light gray) . Ohmic contacts (light gray 
columns) enable bonding wires (not shown) to make electri- 
cal contact to the 2DEG reservoirs, (b)-(c) Scanning elec- 
tron micrographs of a few-electron single-dot device (b) and a 
double-dot device (c), showing the gate electrodes (light gray) 
on top of the surface (dark gray). The white dots indicate the 
location of the quantum dots. Ohmic contacts are shown in 
the corners. White arrows outline the path of current Idot 
from one reservoir through the dot(s) to the other reservoir. 
For the device in (c) , the two gates on the side can be used to 
create two quantum point contacts, which can serve as elec- 
trometers by passing a current Iqpc- Note that this device 
can also be used to define a single dot. Image in (b) courtesy 
of A. Sachrajda. 



a wide range while keeping the tunnel rates high enough 
for measuring electron transport through the dot. 

Applying the same gate design principle to a dou- 
ble quantum dot, Elzerman et al. demonstrated in 
2003 control over the electron number in both dots 
while maintaining tunable t unnel coupling to the reser- 
voir (Elze rman et al I I2003D . Their design is shown in 
Fig. [5t (for more details on design consi deration s and r e- 
lated versions of this gate design, see lHanso 3 (120051) 1. 
In addition to the coupled dots, two quantum point con- 
tacts (QPCs) are incorporated in this device to serve as 
charge sensors. The QPCs are placed close to the dots, 
thus ensuring a good charge sensitivity. This design has 



become the standard for lateral coupled quantum dots 
and is used with m i nor adaptions by several researc h 
groups (jPetta et all |2004 iPioro-Ladriere et all 12005ft : 
one noticable improvement has been the electrical iso- 
lation of the charge sensing part o f the circuit from th e 
reservoirs that connect to the dot ( Hanson et all F2005). 



C. Measurement techniques 

In this review, two all-electrical measurement tech- 
niques are discussed: i) measurement of the current due 
to transport of electrons through the dot, and ii) detec- 
tion of changes in the number of electrons on the dot with 
a nearby electrometer, so-called charge sensing. With the 
latter technique, the dot can be probed non-invasively in 
the sense that no current needs to be sent through the 
dot. 

The potential of ch arge sensing wa s first demonstrated 
in the early 1990s (|Ashoori et all Il992t iField et all 
1993). But whereas current measurements were al- 
ready used e xtensively in the first experi ments on quan- 
tum dots ( Kouwenhoven et Mi 11997ft . charge sens- 



ing has only re cently been full y devel o ped as a spec- 
trosco pic tool ( Elzerman et all l2004at iJohnson et all 



2005a). Several implementations of electrometers cou- 
pled to a quantum dot have been demonstrated: a single 
electr on transistor fabricat e d on top o f the heterostruc- 
ture ( Ashoori et all Il992t iLu et all 12003ft . a second 



electr ostatically defined q uantum dot ( Fuiisawa et all 
l2004lHofmamT et all 1 1995ft and a quantum point contact 
(QPC) ([Field et all 119931 ISprinzak eTail . 12002ft . The 
QPC is the most widely used because of its ease of fabri- 
cation and experimental operation. We discuss the QPC 
operation and charge sensing techniques in more detail 
in section [V] . 

We briefly compare charge sensing to electron trans- 
port measurements. The smallest currents that can be 
resolved in optimized setups and devices are roughly 
10 fA, which sets a lower bound of order 10 fA/e ~ 
100 kHz on the tunnel rate to the reservoir, T, 
for which transpor t expe riments are possible (see e.g. 
IVandersvpen et al\ (2004) for a discussion on noise 
sources). For T < 100 kHz the charge detection technique 
can be used to resolve electron tunneling in real time. Be- 
cause the coupling to the leads is a source of decoherence 
and relaxation (most notably via cotunneling) , charge 
detection is preferred for quantum information purposes 
since it still functions for very small couplings to a (sin- 
gle) reservoir. 

Measurements using either technique are conveniently 
understood with the Constant Interaction model. In the 
next section we use this model to describe the physics of 
single dots and show how relevant spin parameters can 
be extracted from measurements. 
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D. The Constant Interaction model 

We briefly outline the main ingredients of the Con- 
st an^JntCTaction_jrwdejj_for more extensive discussions 
van Houten et all ( 19921 ); iKouwenhoven et ai (| Oil 



see 



119971 ). The nodel is based on two assumptions. First, the 
Coulomb interactions among electrons in the dot, and be- 
tween electrons in the dot and those in the environment, 
are parameterized by a single, constant capacitance, C. 
This capacitance is the sum of the capacitances between 
the dot and the source, Cs, the drain, Cd, and the gate, 
C G : C — Cs + Cd + Cq. (In general, capacitances to 
multiple gates and other parts of the 2DEG will also play 
a role; they can simply be added to C). The second as- 
sumption is that the single-particle energy level spectrum 
is independent of these interactions and therefore of the 
number of electrons. Under these assumptions, the to- 
tal energy U (N) of a dot with N electrons in the ground 
state, with voltages Vs, Vd and V G applied to the source, 
drain and gate respectively, is given by 



U(N) 



-\e\{N-N ) + C S V S + C D V D + C G V G }< 



2C 



N 



(i) 



where — |e| is the electron charge, N \e\ is the charge 
in the dot compensating the positive background charge 
originating from the donors in the heterostructure, and 
B is the applied magnetic field. The terms CsVs, CdVd 
and CgVg can be changed continuously and represent 
an effective induced charge that changes the electrostatic 
potential on the dot. The last term of Eq. [1] is a sum 
over the occupied single-particle energy levels, E n (B), 
which depend on the characteristics of the confinement 
potential. 

The electrochemical potential n(N) of the dot is de- 
fined as: 

(i(N) = U(N) - U(N-l) = 

(N-N -l)E c -^{CsVs+C D V D + C G V G )+E N (2) 
2 e 



Note that the electrochemical potential depends lin- 
early on the gate voltage, whereas the energy has a 
quadratic dependence. In fact, the dependence is the 
same for all N and the whole 'ladder' of electrochemical 
potentials can be moved up or down while the distance 
between levels remains constant 4 . It is this property that 
makes the electrochemical potential the most convenient 
quantity for describing electron tunneling. 

The electrochemical potentials of the transitions be- 
tween successive ground states are spaced by the so-called 
addition energy: 

E add (N) = fi(N+l) - n(N) = E C + AE (4) 

The addition energy consists of a purely electrostatic 
part, the charging energy E G , plus the energy spacing 
between two discrete quantum levels, AE. Note that AE 
can be zero, when two consecutive electrons are added to 
the same spin-degenerate level. 

Electron tunneling through the dot critically depends 
on the alignment of electrochemical potentials in the dot 
with respect to those of the source, /is, and the drain, 
Hd- The application of a bias voltage Vsd = Vs — Vd 
between the source and drain reservoir opens up an en- 
ergy window between /is and [Id of [is — fiD = — \e\VsD- 
This energy window is called the bias window. For en- 
ergies within the bias window, the electron states in one 
reservoir are filled whereas states in the other reservoir 
are empty. Therefore, if there is an 'appropiate' electro- 
chemical potential level within the bias window, electrons 
can tunnel from one reservoir onto the dot and off to the 
empty states in the other reservoir. Here, 'appropriate' 
means that the electrochemical potential corresponds to 
a transition that involves the current state of the quan- 
tum dot. 

In the following, we assume the temperature to be 
negligible compared to the energy level spacing AE (for 
GaAs dots this roughly means T <0.5 K). The size of 
the bias window then separates two regimes: the low- 
bias regime where at most one dot level is within the 
bias window (— |e|Vsr) < AE, E a dd), and the high-bias 
regime where multiple dot levels can be in the bias win- 
dow {-\e\VsD > and/or -\e\V SD > E add ). 



where E G = e 2 /C is the charging energy. The elec- 
trochemical potential contains an electrostatic part (first 
two terms) and a chemical part (last term). Here, fi(N) 
denotes the transition between the iV-electron ground 
state, GS(N), and the (N — l)-electron ground state, 
GS(N-l). When also excited states play a role, we have 
to use a more explicit notation to avoid confusion: the 
electrochemical potential for the transition between the 
(N — l)-electron state |o) and the TV-electron state \b) 
is then denoted as [i a ^b, and is defined as the difference 
in total energy between state \b), Ub{N), and state \a), 
U a (N-l): 



E. Low-bias regime 

For a quantum dot system in equilibrium, electron 
transport is only possible when a level corresponding to 
transport between successive ground states is in the bias 
window, i.e. fis > f*{N) > [Id for at least one value of TV. 



fx a ^ b = U b (N)-U a (N-l) 



(3) 



4 Deviations from this model are sometimes observed in systems 
where the source-drain voltage and gate voltage are varied over 
a very wide range; one notable example being single molecules 
trapped between closely-spaced electrodes, where the capaci- 
tances can depend on the electron state. 
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If this condition is not met, the number of electrons on 
the dot remains fixed and no current flows through the 
dot. This is known as Coulomb blockade. An example of 
such a level alignment is shown in Fig. 



a) 



I h U(A/+-I) |^ 

m 

1 



WN_1)| | 



b) 




Is 








U(N+1) 




IZ 




|i(N) t _ 








U(N-1) 





C) £ 




c 

3 



o 















H 


A/-1 


H 


A/ + 1 


A/-3J 




\ 1 




i ► 



gate voltage \/ G 

FIG. 3 (a)-(b)Schematic diagrams of the electrochemical po- 
tential levels of a quantum dot in the low-bias regime, (a) 
If no level in the dot falls within the bias window set by fis 
and fiD, the electron number is fixed at N— 1 due to Coulomb 
blockade, (b) The fJ,(N) level is in the bias window, so the 
number of electrons can alternate between N — 1 and N, re- 
sulting in a single-electron tunneling current. The magnitude 
of the current depends on the tunnel rate between the dot 
and the reservoir on the l eft, Fs, and on the right, Td (see 
iKouwenhoven et al\ (1997) for details), (c) Schematic plot of 
the current Idot through the dot as a function of gate volt- 
age Vg- The gate voltages where the level alignments of (a) 
and (b) occur are indicated. 

Coulomb blockade can be lifted by changing the volt- 
age applied to the gate electrode, as can be seen from 
Eq. [21 When /x(iV) is in the bias window one extra elec- 
tron can tunnel onto the dot from the source (see Fig.[3jD), 
so that the number of electrons increases from N — 1 to 
N. After it has tunneled to the drain, another electron 
can tunnel onto the dot from the source. This cycle is 
known as single- electron tunneling. 

By sweeping the gate voltage and measuring the cur- 
rent through the dot, Idot, a trace is obtained as shown 
in Fig. [3b. At the positions of the peaks in Idot, an elec- 
trochemical potential level corresponding to transport 
between successive ground states is aligned between the 
source and drain electrochemical potentials and a single- 
electron tunneling current flows. In the valleys between 
the peaks, the number of electrons on the dot is fixed due 
to Coulomb blockade. By tuning the gate voltage from 
one valley to the next one, the number of electrons on the 
dot can be precisely controlled. The distance between 
the peaks corresponds to E a dd (see Eq. 2|) , and therefore 
provides insight into the energy spectrum of the dot. 
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FIG. 4 Schematic diagrams of the electrochemical potential 
levels of a quantum dot in the high-bias regime. The level 
in grey corresponds to a transition involving an excited state, 
(a) Here, Vsd exceeds AE and electrons can now tunnel via 
two levels, (b) Vsd exceeds the addition energy for N elec- 
trons, leading to double-electron tunneling. 



F. High-bias regime 

We now look at the regime where the source-drain bias 
is so high that multiple dot levels can participate in elec- 
tron tunneling. Typically the electrochemical potential 
of only one of the reservoirs is changed in experiments, 
and the other one is kept fixed. Here, we take the drain 
reservoir to be at ground, i.e. [Id — 0. When a negative 
voltage is applied between the source and the drain, /is 
increases (since us = — l e | Vsd)- The levels of the dot 
also increase, due to the capacitive coupling between the 
source and the dot (see Eq.[2|). Again, a current can flow 
only when a level corresponding to a transition between 
ground states falls within the bias window. When Vsd 
is increased further such that also a transition involving 
an excited state falls within the bias window, there are 
two paths available for electrons tunneling through the 
dot (see Fig. In general, this will lead to a change 
in current, enabling us to perform energy spectroscopy of 
the excited states. How exactly the current changes de- 
pends on the tunnel coupling of the two levels involved. 
Increasing Vs d even more eventually leads to a situation 
where the bias window is larger than the addition energy 
(see Fig. 2b>). Here, the electron number can alternate 
between N — 1, N and N+l, leading to a double-electron 
tunneling current. 

We now show how the current spectrum as a function 
of bias and gate voltage can be mapped out. First, the 
electrochemical potentials of all relevant transitions are 
calculated by applying Eq. [3] For example, consider two 
successive ground states, GS(iV) and GS(iV+l), and the 
excited states ES(iV) and ES(iV+l), which are separated 
from the GSs by AE(N) and AE(N+1) respectively (see 
Fig. . The resulting electrochemical potential ladder 
is shown in Fig. (we omit the transition between the 
two ESs). Note that the electrochemical potential of the 
transition ES(7V)^GS(7V + 1) is lower than that of the 
transition between the two ground states. 

The electrochemical potential ladder is used to define 
the gate voltage axis of the (— |e| V^cVg) plot, as in 
Fig. [5b. Here, each transition indicates the gate voltage 
at which its electrochemical potential is aligned with fis 
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ES(A/)<-> GS(N)o GS(A/)<-> v g 
GS(A/+1) GS(/V+1) ES(A/+1) 



FIG. 5 (Color in online edition) (a) Energies for N electrons, U(N), and for N+l electrons, U(N + 1)- Possible transitions 
are indicated by arrows, (b) The electrochemical potential ladder for the transitions depicted in (a), (c) Schematic plot of the 
differential conductance dloor / 'dVsD as a function of — |e| Vsd and Vo- At several positions the level alignment is indicated 
with schematic diagrams. 



and /id at Vsd = 0. Analogous to Fig.[3]>d, sweeping the 
gate voltage at low bias will show electron tunneling only 
at the gate voltage indicated by GS(N) <h> GS{N+1). For 
all other gate voltages the dot is in Coulomb blockade. 

Next, for each transition a V-shaped region is out- 
lined in the (— |e| Vs_D,V(3)-plane, where its electrochem- 
ical potential is within the bias window. This yields 
a plot like Fig. EJ:. The slopes of the two edges of 
the V-shape depend on the capacitances; for Vn = 0, 
the two slopes d(— |e| VsD)/dVo are —Cg/(C — Cs) and 
+Cg/Cs- The transition between the A^-electron GS and 
the (N+ l)-electron GS (black solid line) defines the re- 
gions of Coulomb blockade (outside the V-shape) and 
tunneling (within the V-shape). The other solid lines 
indicate where the current changes due to the onset of 
transitions involving excited states. 

The set of solid lines indicate all the values in the pa- 
rameter space spanned by Vs d and Vq where the current 
Idot changes. Typically, the differential conductance 
dlDOT / dVsD is plotted, which has a nonzero value only 
at the solid lines 5 . 

A general 'rule of thumb' for the positions of the lines 
indicating finite differential conductance is this: if a line 
terminates at the N -electron Coulomb blockade region, 
the transition necessarily involves an N -electron excited 



state. This is true for any N. As a consequence, no lines 
terminate at the Coulomb blockade region where iV=0, 
as there exist no excited state for ^=0 6 . For a transition 
between two excited states, say ES(iV) and ES(iV-f-l), the 
position of the line depends on the energy level spacing: 
for AE(N+1) > AE(N), the line terminates at the (N+ 
l)-electron Coulomb blockade region, and vice versa. 

A measurement as shown in Fig. [5J: is very useful for 
finding the energies of the excited states. Where a line 
of a transition involving one excited state touches the 
Coulomb blockade region, the bias window exactly equals 
the energy level spacing. Figure 0; shows the level dia- 
grams at these special positions for both ES(A^)<->GS(A L f 
1) and GS(iV)<->ES(./V-r-l). Here, the level spacings can 
be read off directly on the — |e| Vsc-axis. 

We briefly discuss the transition ES(A^)^ES(A^+ 1), 
that was neglected in the discussion thus far. The visi- 
bility of such a transition depends on the relative mag- 
nitudes of the tunnel rates and the relaxation rates. 
When the relaxation is much faster than the tunnel rates, 
the dot will effectively be in its ground state all the 
time and the transition ES(-/V)^->ES(iV + 1) can there- 
fore never occur. In the opposite limit where the relax- 
ation is much slower than the tunneling, the transition 
ES(A^)<->ES(Af+I) participates in the electron transport 



In practice, a dependence of the tunnel couplings on Vq and Vg 
may result in a nonzero value of dlpoT/dVgD throughout the 
region where current flows. Since this "background" of nonzero 
dlDOT /dVs d is typically more uniform and much smaller than 
the peaks in dIr>oT/dVsD at the solid lines, the two are easily 
distinguished in experiments. 



Note that energy absorption from the environment can lead to 
exceptions: photon- or phonon-assisted tunneling can give rise 
to lines ending in the N=0 Coulomb blockade region. However, 
many experiments are performed at very low temperatures where 
the number of photons and phonons in thermal equilibrium is ex- 
tremely small. Therefore, these processes are usually negligible. 
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and will be visible in a plot like in Fig. [5]:. Thus, the 
visibility of transitions can give inform ation on the re- 
laxatio n rates between different levels ( Fujisawa et all . 
l2002bh . 

If the voltage is swept across multiple electron tran- 
sitions and for both signs of the bias voltage, the 
Coulomb blockade regions appear as diamond shapes 
in the (— |e| Vsd,Vg) -plane. These are the well-known 
Coulomb diamonds. 



III. SPIN SPECTROSCOPY METHODS 

In this section, we discuss various methods for get- 
ting information on the spin state of the electrons on a 
quantum dot. These methods make use of various spin- 
dependent energy terms. First, each electron spin is in- 
fluenced directly by an external magnetic field via the 
Zeeman energy Ez = S z gfi B B where S z is the spin z- 
component. Moreover, the Pauli exclusion principle for- 
bids two electrons with equal spin orientation to occupy 
the same orbital, thus forcing one of the electrons into a 
different orbital. This generally leads to a state with a 
different energy. Finally, the Coulomb interaction leads 
to an energy difference (the exchange energy) between 
states with symmetric and anti-symmetric orbital wave- 
functions. Since the total wavefunction of the electrons 
is anti-symmetric, the symmetry of the orbital part is 
linked to that of the spin. 



A. Spin filling derived from magnetospectroscopy 



Assuming S z only changes by ±i, the possible out- 
comes and the corresponding filling schemes are 



dE add {N) 
dB 



: T,T or 4,4 
-ffMB : i,T) 



where the first (second) arrow depicts the spin added 
in the N - 2 -> N - 1 (N-l -> N) electron transi- 



tion. Spin filling of both vertical (ISasaki et al 



1998) 



2000: 



and lateral GaA s quantum dots dDuncan et al.. . , 
iLindemann et all . 120021 : iPotok et att , 20031 ) hasbeen de- 
termined using this method, showing clear deviations 
from a simple "Pauli" filling (S z alternating between 
and i). Note that transitions where S z of the ground 
state changes by more than ~, which can occur due 
to many-body interactions in the dot, can lead to a 



spin blockade of the c urrent ( Korkusihski et all . |2004 
IWeinmann et aLl . ll995l ). 

In circularly symmetric few-electron vertical dots, spin 
states have been determined from the evolution of orbital 
states in a magnetic field perpendicular to the plane of 
the dots. This indirect determination of the spin state 
has allowed the observation of a two-electron singlet-to- 
triplet ground state transition and a four-electron spin 
filling following H und's rule. For a re view on these 
experiments, see IKouwenhoven et all ( 200lh . Simi- 
lar techniques were also used in experiments on few- 
electr on lateral dot s in both weak and strong m agnetic 
fields ( Ciorga et all . 120001 : iKvriakidis et "oil. 1200211 . 



The spin filling of a quantum dot can be derived from 
the Zeeman energy shift of the Coulomb peaks in a mag- 
netic field. (An in-plane magnetic field orientation is fa- 
vored to ensure minimum disturbance of the orbital lev- 
els). On adding the iVth electron, the z-component S z 
of the spin on the dot is either increased by 1/2 (if a 
spin-up electron is added) or decreased by 1/2 (if a spin- 
down electron is added). This change in spin is reflected 
in the magnetic field dependence of the electrochemical 
potential n(N) via the Zeeman term 



9[ i B B [S Z {N) - S„{N-1)] = g[i B B [AS Z (N)] . (5) 

As the g- factor in GaAs is negative (see Appendix [A"|l . 
addition of a spin-up electron (AS z (N)—+l/2) results 
in /z(iV) decreasing with increasing B. Spin-independent 
shifts of n(N) with B (e.g. due to a change in con- 
finement potential) are removed by lookin g at the de- 
pend ence of the addition energy E ac id on B ( Weis et all . 
11993b : 



dE add (N) = d f i(N) dfi(N-l) 
dB dB dB 

= gt i B [AS z (N)-AS z (N-l)]. (6) 



B. Spin filling derived from excited-state spectroscopy 

Spin filling can also be deduced from excited- 
state spectroscopy witho ut changing the magnetic 
field i Cobden et all . fl998). provided the Zeeman en- 
ergy splitting AEz = 2\E%\ = g[i B B between spin-up 
and spin-down electrons can be resolved. This power- 
ful method is based on the simple fact that any single- 
particle orbital can be occupied by at most two electrons 
due to Pauli's exclusion principle. Therefore, as we add 
one electron to a dot containing N electrons, there are 
only two scenarios possible: either the electron moves 
into an empty orbital, or it moves into an orbital that 
already holds one electron. As we show below, these 
scenarios always correspond to ground state filling with 
spin-up and spin-down, respectively. 

First consider an electron entering an empty orbital 
with well-resolved spin splitting (see Fig. [6ji). Here, ad- 
dition of a spin-up electron corresponds to the transition 
GS(N) <-> GS(N+1). In contrast, addition of a spin- 
down electron takes the dot from GS(N) to ES(N+1), 
which is AEz higher in energy than GS(N+1). Thus we 
expect a high-bias spectrum as in Fig. [5^. 

Now consider the case where the (N + l)th electron 
moves into an orbital that already contains one electron 



10 




GS(A/)^GS(A/+1) GS(A/)^ES(A/+1) 




ES(A/)^GS(A/+1) GS(A/)^GS(A/+1) 



FIG. 6 Spin filling deduced from high-bias excited-state spec- 
troscopy. Shown are schematic diagrams of dlooT /dVsn in 
the {Vsd, VbO-prane. (a) Ground-state filling is spin-up: a 
line corresponding to an (iV+l)-electron ES, separated from 
GS(7V+1) by AE Z , terminates at the edge of the (N+Vj- 
electron Coulomb blockade region (point P). (b) Ground-state 
filling is spin-down: a line corresponding to an iV-electron ES, 
separated from GS(iV) by AEz, terminates at the iV-electron 
Coulomb blockade region (point Q). 



(see Fig.[B}3). The two electrons need to have anti-parallel 
spins, in order to satisfy the Pauli exclusion principle. 
If the dot is in the ground state, the electron already 
present in this orbital has spin-up. Therefore, the elec- 
tron added in the transition from GS(N) to GS(N + 1) 
must have spin-down. A spin-up electron can only be 
added if the first electron has spin-down, i.e. when 
the dot starts from ES(N), AEz higher in energy than 
GS(N). The high-bias spectrum that follows is shown 
schematically in Fig. [6Jd. 

Comparing the two scenarios, we see that the spin fill- 
ing has a one-to-one correspondence with the excited- 
state spectrum: if the spin ES line terminates at the 
(A r +l)-electron Coulomb blockade region (as point P in 
Fig. [6^), a spin- up electron is added to the GS; if however 
the spin ES line terminates at the TV-electron Coulomb 
blockade region (as point Q in Fig. [Us), a spin-down elec- 
tron is added to the GS. 

The method is valid regardless of the spin of the ground 
states involved, as long as the addition of one electron 



changes the spin ^-component of the ground state by 
| AS, | = 1/2. If | AS* | > 1/2, the (N + l)-electron GS 
cannot be reached from the ./V-electron GS by addition 
of a single electron. This would caus e a spin blockade 
of ele ctron transport through the dot ( Weinmann et all 
Il995h . 



C. Other methods 

If the tunnel rates for spin-up and spin-down are not 
equal, the amplitude of the current can be used to de- 
termine the spin filling. This method has been termed 
spin-blockade spectroscopy. This name is slightly mis- 
leading as the current is not actually blocked, but rather 
assumes a finite value that depends on the spin orienta- 
tion of the transported electrons. This method has been 
demonstrated and utilized in the quantum Hall regime, 
where the spatial separation of spin-split edge channels 
induces a large difference in the tunnel rates o f spin- 
up and spin-down ele ctrons ( Ciorga et all l2002l l2000t 
iKupidura et all [20061 ). Spin-polarized leads can also be 
obtained in moderate magnetic fields by changing the 
electron density near the dot with a gate. This concept 
was used to perform spin spectroscopy on a quantum dot 
conn ected to gate-tunab le quasi-one-dimensional chan- 
nels dHitachi et a/.l . l2006l ). 

Care must be taken when inferring the spin filling from 
the amplitude of the current as other factors, such as 
the orbital spread of the wavefunction, can have a large, 
even dominating influence on the current amplitude. A 
prime example is the difference in tunnel rate between 
the two-electron spin singlet and triplet states due to 
the different orbital wavefunctions of these states. In 
fact, this difference is large enough to allow single-shot 
readout of the two-electron spin state, as will be discussed 
in Section |VLCl 

In zero magnetic field, a state with total spin S is 
(2S + l)-fold degenerate. This degeneracy is reflected 
in the current if the dot has strongly asymmetric barri- 
ers. As an example, in the transition from a one-electron 
S'=l/2 state to a two-electron S=0 state, only a spin-up 
electron can tunnel onto the dot if the electron that is al- 
ready on the dot is spin-down, and vice-versa. However, 
in the reverse transition (S—Q to S=l/2), both electrons 
on the dot can tunnel off. Therefore, the rate for tun- 
neling off the dot is twice the rate for tunneling onto the 
dot. In general, the ratio of the currents in opposite bias 
directions at the GS(N) <-* GS(N+l) transition is, for 
spin-independent tunnel rates and for strongly as ymmet- 
ric ba rriers, given by [2S(N+1) + 1]/[2S(N) + 1] (jAkeral . 
Il999h . Here, S(N) and S{N+1) denote the total spin of 
GS(N) and GS(N+l) respectively. This relation can be 
used in experiments to determine the ground s tate total 
spin (jCobden et all Il998t iHavashi et aZ.l . l2003h . 

Information on the spin of the ground state 
can also be found from (inelastic) cotunneling 
currents (K ogan et all 12004 ) or the current due 
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FIG. 7 (a)-(f) Excited-state spectroscopy on two devices: device A is fabricated on a heterostructure with the 2DEG at 90 nm 
below the surface and device B with the 2DEG at 60 nm below the surface. Differential conductance dloor / dVsD is plotted 
as a function of Vsd and gate voltage near the 0*-+l electron transition, for in-plane magnetic fields (as indicated in top right 
corners). Darker corresponds to larger dlnor / 'dVsD- Data on device A shows spin splitting in both the orbital ground and first 
excited state; data on device B only displays the orbital ground state, (g) Zeeman splitting AEz as a function of B extracted 
from (a)-(f) and similar measurements. Gra y solid lines are fits to the data. T h e dashed line shows A Ez expected for the bulk 
GaAs g-factor of -0.44. Data adapted from Wi llems van Beveren et al\ (|2005l ); lHanson et~al\ ||2003f ) 



to a Kondo resonance (ICronenwett et al. 

I, 119981: 

iGoldhaber-Gordon et all Il998f) . If a magnetic field 
B drives the onset of these currents to values of 
Vsd = ±9HbB/ |e|, it follows that the ground state has 
nonzero spin. Since the processes in these currents can 
change the spin z-component by at most 1, the absolute 
value of the spin can not be deduced with this method, 
unless the spin is zero. 

We end this section with some remarks on spin filling. 
First, the parity of the electron number can not be in- 
ferred from spin filling unless the sequence of spin filling 
is exactly known. For example, consider the case where 
the electron added in the GS(iV)^GS(iV+l) transition 
has spin-down. Then, if the dot follows an alternating 
(Pauli) spin filling scheme, TV is odd. However, if there is 
a deviation from this scheme such that GS(N) is a spin 
triplet state (total spin S—l), then N is even. 

Second, spin filling measurements do not yield the 
absolute spin of the ground states, but only the 
change in ground state spin. However, by starting 
from zero electrons (and thus zero spin) and track- 
ing the change in spin at subsequent electron transi- 
tions, t he total spin of the ground state can be deter- 
mined ( Willems van Beveren et oil 12005b . 



potentials: 



IV. SPIN STATES IN A SINGLE DOT 

A. One-electron spin states 

The simplest spin system is that of a single electron, 
which can have one of only two orientations: spin-up or 
spin-down. Let E^q and £^0 {Ef t i and i?x,i) denote 
the one-electron energies for the two spin states in the 
lowest (first excited) orbital. With a suitable choice of the 
zero of energy we arrive at the following electrochemical 
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where AE or b is the orbital level spacing. 

Figures show excited-state spectroscopy measure- 
ments on two devices, A and B, via electron transport 
at the N—0 <-> 1 transition, at different magnetic fields 
Bn applied in the plane of the 2DEG. A clear split- 
ting of both the orbital ground and first excited state 
is obs erved, which increases wit h increasing magnetic 
field (IWillems van Beveren et'al. I. 120051: lHanson et all 
l2003t lK5nemann et all 120051: iPotok et all I2003D . The 
orbital level spacing AE or i> in device A is about 1.1 meV. 
Comparison with Fig. [5] shows that a spin-up electron is 
added to the empty dot to form the one-electron ground 
state, as expected. 

In Fig.[7g the Zeeman splitting AEz is plotted as func- 
tion of Bn for the same two devices, A and B, which 
are made on different heterostructures. These measure- 
ments allow a straightforward determination of the elec- 
tron g- factor. The measured g-factor can be affected 
by: (i) extension of the electron wave function into the 
Al n . 3 Gap. 7 As r egion where g = +0.4 ( Salis et all , l200lt 



ISnelling et all Il99lt) . (ii) thermal nuclear polarization 
which decreases the eff ective magnetic field thro u gh th e 
hyperfine interaction (jMeier and Zakharchenval 1984b . 
(iii) dynamic nuclear polarization due to electron- 
nuclear flip-flop processes in the dot, which enhances 
the e ffective magnetic field (jMeier and Zakharchenval 
Il984l ). (iv) the nonparabolicity of the GaAs conduc- 
tion b and (ISnelling et ~al. , [1991). (v) the spin-orbit cou- 
pling (jFalko et all 120051 ) . and (vi) the confinement 
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FIG. 8 (Color in online edition) (a) Energy diagram schematically showing the energy levels of the one- and two-electron states. 
The allowed transitions between these levels are indicated by arrows, (b) Electrochemical potential ladder corresponding to the 
transitions shown in (a), using the same color coding. Changing the gate voltage shifts the ladder as a whole. Note that the 
three triplet states appear at only two values of the electrochemical potential, (c) Energetically allowed l<-»2 electron transitions 
as a function of Vsd and Vg- The lines corresponding to j^S 1 outline the region of transport; outside this region, where lines 
are dashed, the dot is in Coulomb blockade, (d) dloor / dVsD as a function of Vg and Vsd around the l<-»2 electron transition 
at B a — 12 T in device A. The regions labelled with letters A- F correspond well to those in (c). In the region labeled A only 
spin-down electrons pass through the dot. Data adapted from Han son et al\ (|2004aT ). 



poten tial ( Biork et all 120051 : iHermann and Weisbuchl . 
Il977l ). The effect of the nuclear field on the measured 
(/-factor is discussed in more detail in Appendix [X] More 
experiments are needed to separate these effects, e.g. by 
measuring the dependence of the g-factor on the orien- 
tation of the in-plane magneti c field with respect to the 



crystal axis (jFalko et al 



aagnetic 



B. Two-electron spin states 

The ground state of a two-electron dot in zero magnetic 
field is always a spin singlet (total sp in quantum number 
5 = 0) (jAshcroft and Merminl . Il974f) . formed by the two 
electrons occupying the lowest orbital with their spins 
anti-parallel: | S) = (| tl ) - | IT »/v / 2. The first excited 
states are the spin triplets (S =1), where the antisymme- 
try of the total two-electron wave function requires one 
electron to occupy a higher orbital. Both the antisym- 
metry of the orbital part of the wavefunction and the 
occupation of different orbitals reduce the Coulomb en- 
ergy of the triplet states with respect to the singlet with 



two e lectrons in the same orbital (jKouwenhoven et all 
l200ll ). We include this change in Coulomb energy by the 
energy term Ek- The three triplet states are degenerate 
at zero magnetic field, but acquire different Zeeman en- 
ergy shifts Ez in finite magnetic fields because their spin 
z-components differ: S z = +1 for | T + ) = | f f ) , S z = for 
|T ) = (|U) + UT»/>/2 and S z = -1 for |T_) = |U>- 

Using the Constant Interaction model, the energies of 
the states can be expressed in terms of the single-particle 
energies of the two electrons plus a charging energy Eq 
which accounts for the Coulomb interactions: 

Us =E^ t o + E^q + Ec = 2E^q + AEz + Ec 



Ut + — E-^fl + E 



J 1A 



E 



K 



E, 



c 



= 2E^o + AEorf, — Ek + Ec 

= 2Ei > o + Est +Ec 
U To =E^ q +E^ + E S t + Ec 

= 2E^ + AE orb — E K + AE z + E c 

= 2E^Q-\-EsT~\~AEz J rEc 
Vt_ = 2E^q + Est + Ec 

= 2E^q + AE or i, — Ek + 2AEz + Ec 
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= 2E h0 + EsT + 2AE z +Ec, 

with Est denoting the singlet-triplet energy difference 
in the absence of the Zeeman splitting AEz'- Est = 

AE or b — Ek ■ 

We first consider the case of an in-plane magnetic field 
Bn. Here, Est is almost independent of Bn and the 
ground state remains a spin singlet for all fields attainable 
in the lab. The case of a magnetic field perpendicular to 
the plane of the 2DEG will be treated below. 

Fig. [8^1 shows the possible transitions between the 
one-electron spin-split orbital ground state and the two- 
electron states. The transitions f<-> T_ and J.<-> T + 
are omitted, since these require a change in the spin 
z-component of more th an | and are thus spin-blocked 
(jWeinmann et all Il995l ). .^rom the energy diagram the 
electrochemical potentials can be deduced (see Fig. [8p): 
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Note that ^t,o*-»t+ = /*|,o«r„ and ^,o^r = AH,owr_- 
Consequently, the three triplet states change the first- 
order transport through the dot at only two values of 
Vsd- The reason is that the first-order transport probes 
the energy difference between states with successive elec- 
tron number. In contrast, the onset of second-order 
(cotunneling) currents is governed by the energy differ- 
ence between states with the same number of electrons. 
Therefore, the triplet states change the second-order (co- 
tunneling) currents at three values of V s d if the ground 
state is a singlet 7 ([Paaske et 

In Fig. [SJ: we map out the positions of the electro- 
chemical potentials as a function of Vc and Vsd- For 
each transition, the two lines originating at Vsd = 
span a V-shaped region where the corresponding elec- 
trochemical potential is in the bias window. In the re- 
gion labeled A, only transitions between the one-electron 
ground state, | T , ) , and the two-electron ground state, 
\S), are possible, since only /if,o<-»s is positioned inside 
the bias window. In the other regions several more tran- 
sitions are possible which leads to a more complex, but 
still understandable behavior of the current. Outside the 
V-shaped region spanned by the ground state transition 
/i|fl«s, Coulomb blockade prohibits first order electron 
transport. 



Experimental results from device A, shown in Fig. [8ji, 
are in excellent agreement with the predictions of Fig.[5J;. 
Comparison of the data with Fig. [6] indicates that indeed 
a spin-down electron is added to the one-electron (spin- 
up) ground state to form the two-electron singlet ground 
state. From the data the singlet-triplet energy difference 
Est is found to be w 520 ^eV. The fact that Est is about 
half the single-particle level spacing (AE or b = 1 meV) 
indicates the importance of Coulomb interactions. The 
Zeeman energy, and therefore the g-factor, is found to 
be the same for the one-electron states as for the two- 
electron states (within the measurement accuracy of w 
5%) on both device A and B. We note that the large 
variation in differential conductance observed in Fig. [5Ji, 
can be explained by a sequential tunneli ng model with 
spin- a nd orbital-dependent tunnel rates ( Hanson et all 
20043). 
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7 If the ground state is a triplet, the cotunneling current only 
changes at two values of Vsd (0 and AEz/\e\), due to the spin 
selection rules. 



FIG. 9 Single-triplet ground state transition in a two-electron 
quantum dot. (a) Differential conductance cIIdot / 'cLVsd ver- 
sus gate voltage, Va, and perpendicular magnetic field, B±. 
Dark (light) corresponds to high (low) value for cIIdot /dVsD- 
Within the stripe of finite conductance, set by the source- 
drain bias voltage, the evolution of the energy difference be- 
tween the singlet state (ground state at zero field) and the 
triplet state is visible. At around 1.1 T the singlet and triplet 
states cross and the ground state becomes a spin triplet, (b) 
Energy difference between the singlet and the triplet states, 
Est, as a function of B±, ex tracted from (a). Data adapted 
from iKvriakidis et all <|2002t) . 

By applying a large magnetic field perpendicular to 
the plane of the 2DEG a spin singlet-triplet ground state 
transition can be induced, see Fig. [9j This transition 
is driven by two effects: (i) the magnetic field reduces 
the energy spacing between the ground and first ex- 
cited orbital state and (ii) the magnetic field increases 
the Coulomb interactions which are larger for two elec- 
trons in a single orbital (as in the singlet state) than 
for two electrons in different orbitals (as in a triplet 
state). Singl et-triplet transitions were first observed i n 
vertical dots ( Kouwenhoven et "oD . 120011 : ISu et gj. , 1992). 
In lateral dots, the gate-voltage dependence of the con- 
finement potential has allowed electrical tuning of the 
singlet-triplet t ransit ion field ( Kvriakidis et all l2002t 
IZumbiihl efal\ , \2004 . 

In very asymmetric lateral confining potentials with 
large Coulomb interaction energies, the simple single- 
particle picture breaks down. Instead, the two electrons 
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in the ground state spin singlet in such dots will tend to 
avoid each other spatially, thus forming a quasi-double 
dot state. Experiments and calculations indicating this 
double-d ot-like behaviour in asymmetric dots have bee n 
reported (|Ellenberger et aZj . l2006tlZumbuhl eFaZl . l2004l ). 



C. Quantum dot operated as a bipolar spin filter 

If the Zeeman splitting exceeds the width of the en- 
ergy levels (which in most cases is set by the thermal en- 
ergy), electron transport through the dot is (for certain 
regimes) spin -polarized and the dot can be operated a s 
a spin filter (jHanson et all l2004al iRecher et all . l2000j b 
In particular, the electrons are spin-up polarized at the 
N = <-> 1 transition when only the one-electron spin- 
up state is energetically accessible, as in Fig.fTOb. At the 
N = 1 <-> 2 transition, the current is spin-down polarized 
if no excited states are accessible (region A in Fig. [5b), 
see Fig. [TUb . Thus, the polarization of the spin filter can 
be reversed electrically, by tuning the dot to the relevant 
transition. 




b) 



T<->S 



FIG. 10 Few-electron quantum dot operated as a bipolar spin 
filter. Schematic diagrams show the level arrangement for 
ground state transport at (a) the 1 electron transition, 
where the dot filters for spin-up electrons, and (b) at the 
1 <-> 2 electron transition, where the dot only transmits spin- 
down electrons. 

Spectroscopy on dots containing more than two elec- 
trons has shown important deviations from an alternat- 
ing spin filling scheme. Already for four electrons, a spin 
ground state with total spin S=l i n zero magnetic field 



has b een observed in both v ertical ([K ouwcnhovc n et al. 
2001) and lateral dots (jWillems van Beveren et al. 
2005). 



V. CHARGE SENSING TECHNIQUES 

The use of local charge sensors to determine the num- 
ber of electrons in single or double quantum dots is a re- 
cent technological improvement that has enabled a num- 
ber of experiments that would have been difficult, or im- 
possible to per form using s tanda rd electrical transport 
measurements ( Field etal\ , Il993l ). In this section, we 
briefly discuss relevant measurement techniques based on 
charge sensing. Much of the same information as found 
by measuring the current can be extracted from a mea- 
surement of the charge on the dot, Qdot, using a nearby 



electrometer, such as a quantum point contact (QPC). 
In contrast to a measurement of the current through the 
dot, a charge measurement can be also used if the dot is 
connected to only one reservoir. 
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FIG. 11 Quantum point contact operated as an electrometer. 
A typical device, with the current paths through the dot and 
through the QPC, is shown in Fig. [5fc. (a) QPC conductance 
Gqpc vs. gate voltage on one of the two gates that defines 
the QPC, Vg,qpc- Halfway the last conductance step, at 
Gqpc ~ e 2 /h (indicated by a cross), the QPC is very sensitive 
to the charge on the dot. (b) Direct comparison between 
current measurement (top p anel) and charge sensin g (bottom 
panel). Data adapted from Elzerman 

The conductance Gqpc through a QPC is quan- 
tized (jvan Wees et all Il988t IWharam etal\ . 1 19881) . At 
the transitions between quantized conductance plateaus, 
Gqpc is very sensitive to the electrostatic environment 
including the number of electrons iV on a nearby quan- 
tum dot (see Fig. [TTh). This property can be ex- 
ploited to determine the absolu te number of electrons 
in si ngle (|Sprinzak et all l2002t) and coupled quantum 



dots ( Elzerman et al. , 20031 ). even when the tunnel cou- 



pling is so small that no current through the dot is 
detected. Figure II lb shows measurements of the cur- 
rent and of dGc-pc /dVc over the same range of Vq- 
Dips in dGQpc I dVc coincide with the current peaks, 
demonstrating the validity of charge sensing. The sign 
of dGQpe /dVc is understood as follows. On increasing 
Vg, an electron is added to the dot. The electric field cre- 
ated by this extra electron reduces the conductance of the 
QPC, and therefore dGo-pc/dVc dips. The sensitivity of 
the charge sensor to changes in the dot char ge can be op- 
timiz ed using an appropriate gate design ( Zhang et all 
I2004D . 

We should mention here that charge sensing fails when 
the tunnel time is much longer than the measurement 
time. In this case, no change in electron number will be 
observed when the gate voltage i s swept and the e quilib - 
rium charge can not be probed ( Rushforth et all |2004[) . 
Note that a quantum dot with very large tunnel barriers 
can trap electrons for minutes or even hour s under non- 
equilibrium conditions ( Cooper et aZ.I . [2000h . This again 
emphasizes the importance of tunable tunnel barriers (see 
Section Hl.B[) . Whereas the regime where the tunnel time 
largely exceeds the measurement time is of little interest 
for this review, the regime where the two are of the same 
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order is actually quite useful, as we explain below. 

We can get information on the dot energy level spec- 
trum from QPC measurements, by monitoring the aver- 
age charge on the dot while applying short gate voltage 
pulses that bring the dot out of its charge equilibrium. 
This is the case when the voltage pulse pulls fi(N) from 
above to below the electrochemical potential of the reser- 
voir fires- During the pulse with amplitude Vp > 0, the 
lowest energy state is GS(N), whereas when the pulse is 
off (V P = 0), the lowest energy state is GS(N - 1). If 
the pulse length is much longer than the tunnel time, the 
dot will effectively always be in the lowest-energy charge 
configuration. This means that the number of electrons 
fluctuates between N — 1 and N at the pulse frequency. 
If, however, the pulse length is much shorter than the 
tunnel time, the equilibrium charge state is not reached 
during the pulse and the number of electrons will not 
change. Measuring the average value of the dot charge 
as a function of the pulse length thus yields information 
on the tunnel time. In between the two limits, i.e. when 
the pulse length is comparable to the tunnel time, the av- 
erage value of the dot charge is very sensitive to changes 
in the tunnel rate. 
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FIG. 12 Excited-state spectroscopy on a one-electron dot us- 
ing charge sensing, (a) SGqpc at / = 385 Hz versus Vg, with 
Vp = 6 mV. Here, V « 2.4 kHz. (b) Schematic electrochemi- 
cal potential diagrams for the case that only the GS is pulsed 
across fires- (c) Idem when both the GS and an ES are pulsed 
across fires- (d) Derivative of SGqpc with respect to Vg plot- 
ted as a function of Vg and Vp. Note that here Vp is negative, 
and therefore the region of tunneling extends to more positive 
gate voltage as \Vp\ is increased. Th e curve in (a) is taken a t 
the dotted line. Data adapted from lElzerman" et al\ 11200481 ). 

In this situation, excited-state spectroscopy 
can be performed by raising the pulse amp litude 
Vp ([Elzerman et all l2004al I Johnson et all l2005a|) . For 
small pulse amplitudes, at most one level is available for 
tunneling on and off the dot, as in Fig. [12b . Whenever 
Vp is increased such that an extra transition becomes 
energetically possible, the effective tunnel rate changes 
as in Fig[l"2"b. This change is reflected in the average 



value of the dot charge and can therefore be measured 
using the charge sensor. 

The signal-to-noise ratio is enhanced significantly by 
lock-in detection of Gqpc at the pulse frequency, thus 
measuring the avera ge change in Gnpc when a volt- 
age pulse is applied (|Sprinzak et all I2002T ) ■ We denote 
this quantity by SGqpc- Figure [l"2"h shows such a mea- 
surement of SGqpc, lock-in detected at the pulse fre- 
quency, as a function of Vc around the <-> 1 elec- 
tron transition. The different sections of the dip corre- 
spond to FigsIT2"b and c as indicated, where GS (ES) 
is the electrochemical potential of the <->t (0 
transition. Figure [T2H shows a plot of the derivative 
of SGqpc with respect to Vg in the (Vp, Vc<)-plane, 
where the one-electron Zeeman splitting is clearly re- 
solved. This measurement is analogous to increasing 
the source-drain bias Vsd in a transport measurement, 
and therefore leads to a similar plot as in Fig. with 
Vsd replaced b y Vp, and dlpgr I 'dVsp replaced by 
d(SGq pc)ldVa (jElzerman et all [2004al iFuiisawa et all 
l2002bi 
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FIG. 13 (a) Measured changes in the QPC current, AIqpc, 
when a pulse is applied to a gate, near the degeneracy point 
between and 1 electrons on the dot (bias voltage across 
the QPC is 1 mV). The pulse of positive voltage increases 
the QPC current due to the capacitive coupling between the 
pulsed gate and the QPC. Shortly after the start of the pulse, 
an electron tunnels onto the dot and the QPC current de- 
creases. When the pulse has ended, the electron tunnels off 
the dot again, (b) Average of 286 traces as in (a). The top 
and bottom panel are taken with a different gate settings, and 
therefore different tunnel rates are observed. The damped 
oscillation following the pulse edges is d ue to the 8th-order 
40 kH z filter used. Data adapted from IVandersypen et al\ 
(|2004f ). 

The QPC response as a function of pulse length is 
a unique function of tunnel rate. Therefore, compar- 
ison of the obtained response function with the theo- 
rctica l function yield s an accurate v a lue o f the tunnel 
rate (jElzerman et qL, 2004a; H anso nl l2005h . In a dou- 
ble dot, charge sensing can be us ed to quantitat i vely set 
the ratio of the tunnel rates (see lJohnson et all ( 2005bl) 
for det ails), and also to obser ve the direction of tunnel 
events ( Fuiisawa et all l2006bl ). 

Electron tunneling can be observed in real time if 
the time between tunnel events is longer than the time 
needed to determine the number of electrons on the dot 
- or equivalently: if the bandwidth of the charge de- 
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tection exceeds the tunnel rate and the signal from a 
single elect ron charge exceeds the noise level ove r that 
bandwidth ( Lu et fflfl. l200fllScnoelkor)f et ojlll998T i. Fig- 
ure 113a shows gate-pulse-induced electron tunneling in 
real time. In Fig. [T3b . the average of many such traces 
is displayed; from the exponential decay of the signal the 
tunnel rate can be accurately determined. 

Optimized charge sensing setups typically have a band- 
width that allow s tunneling to b e obser v ed on a m ic rosec- 
ond timescale dFul isawa et al 
Sch leser et 



. .20041: lLu~et all 120031: 

2004IVandersvpen et aLl . l2004r i. If the re- 
laxation of the electron spin occurs on a longer timescale, 
single-shot readout of the spin state becomes possible. 
This is the subject of the next section. 



VI. SINGLE-SHOT READOUT OF ELECTRON SPINS 

A. Spin-to-charge conversion 

The ability to measure individual quantum states in 
a single-shot mode is important both for fundamental 
science and for possible applications in quantum infor- 
mation processing. Single-shot immediately implies that 
the measurement must have high fidelity (ideally 100%) 
since only one copy of the state is available and no aver- 
aging is possible. 

Because of the tiny magnetic moment associated 
with the electron spin it is very challenging to measure 
it directly. However, by correlating the spin states 
to different charge states and subsequently measuring 
the cha rge on the dot, the s pin s tate can be deter- 
mined (Loss and DiVincenzol . [l998). This way, the 
measurement of a single spin is replaced by the mea- 
surement of a single charge, which is a much easier task. 
Several schemes for such a spin - to- charge conversi on 
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Vandersvpen et all l2002h . Two methods, both outlined 



energy-selective readout (E-RO), the spin levels are posi- 
tioned around the electrochemical potential of the reser- 
voir fj, res (see Fig. 1 14a ). such that one electron can tunnel 
off the dot from the spin excited state, \ES), whereas 
tunneling from the ground state, \GS), is energetically 
forbidden. Therefore, if the charge measurement shows 
that one electron tunnels off the dot, the state was \ES), 
while if no electron tunnels the state was \GS). This 
readout concept was pioneered by Fujisawa et al. in 
a series of transport experiments, where the measured 
current reflected the a yerage state of t he elec tron after 
a pulse sequence (see iFui isawa et Ml (l2006ah for a re- 
view). Using this pump-probe technique, the orbital re- 
laxation time and a lower bound on the spin relaxation 
time in f ew-electron ver t ical and later al dots was de- 
term ined (jFuiisawa et a/.l . [2002all2001allb1 : lHanson et all 
I2003D . A variation of E-RO can be used for reading out 
the two-electron spin states in a double dot (see Sec- 
tion EnEE}. 

Alternatively, spin-to-charge conversion can be 
achieved by exploiting the difference in tunnel rates of 
the different spin states to the reservoir. We outline the 
concept of this tunnel-rate-selective readout (TR— RO) 
in Fig. [14b . Suppose that the tunnel rate from \ES) to 
the reservoir, Tes, is much higher than the tunnel rate 
from \GS), Tqs, i- e - ^es ^> ^gs- Then, the spin state 
can be read out as follows. At time t=0, the levels of 
both \ES) and \GS) are positioned far above n res , so 
that one electron is energetically allowed to tunnel off the 
dot regardless of the spin state. Then, at a time t = r, 
where 3> r ^> r^ji, an electron will have tunneled 
off the dot with a very high probability if the state was 
\ES), but most likely no tunneling will have occurred if 
the state was \GS). Thus, the spin information is con- 
verted to charge information, and a measurement of the 
number of electrons on the dot reveals the original spin 
state. The TR-RO can be used in a similar way if T es 
is much lower than Tqs- A conceptually similar scheme 
has allowed single-shot readout of a superconducting 
charge qubit ( Astafiev et all l2~004r i . 



in Fig. 1141 have been experimentally demonstrated. 
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FIG. 14 Energy diagrams depicting two different methods 
for spin-to-charge conversion: (a) Energy-selective readout 
(E-RO) and (b) Tunnel-rate-selective readout (TR-RO). 

In one method, a difference in energy between the 
spin states is used for spin-to-charge conversion. In this 



B. Single-shot spin readout using a difference in energy 

Single-shot readout of a single electron spin has 
first been demonstra t ed us ing the E-RO tech- 
nique (|Elzerman et all l2004bl) . In this section we 
discuss this experiment in more detail. 

A quantum dot containing zero or one electrons is tun- 
nel coupled to a single reservoir and electrostatically cou- 
pled to a QPC that serves as an electrometer. The elec- 
trometer can determine the number of electrons on the 
dot in about 10 [is. The Zeeman splitting is much larger 
than the thermal broadening in the reservoir. The read- 
out configuration therefore is as in Fig. EH, with the 
•f->t transition as \GS) and the *-»J. transition as 
\ES). In the following, we will also use just f and j to 
denote these transitions. 
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FIG. 15 Scheme for E-RO of a single electron spin, (a) Two- 
level voltage pulse scheme. The pulse level is 10 mV during 
twait and 5 mV during t rea d (which is 0.5 ms for all measure- 
ments), (b) Schematic response of the QPC if the injected 
electron has spin-j (solid line) or spin- J. (dotted line; the dif- 
ference with the solid line is only seen during the read-out 
stage). Arrows indicate the moment an electron tunnels into 
or out of the quantum dot. (c) Energy diagrams for spin-up 
(Ei) and spin-down (E±) during the different stages of the 
pulse. If the spin is up at the start of the read-out stage, 
no change in the charge on the dot occurs during t rea£ j. In 
contrast, if the spin is down, the electron can escape and be 
replaced by a spin-up electron. 



To test the single-spin measurement technique, the fol- 
lowing three-stage procedure is used: 1) empty the dot, 
2) inject one electron with unknown spin, and 3) measure 
its spin state. The different stages are controlled by gate 
voltage pulses as in Fig. fT5"b , which shift the dot's en- 
ergy levels as shown in Fig.[l5b. Before the pulse the dot 
is empty, as both the spin-up and spin-down levels are 
above the electrochemical potential of the reservoir \x re . s . 
Then a voltage pulse pulls both levels below fj, res . It is 
now energetically allowed for one electron to tunnel onto 
the dot, which will happen after a typical time ss 
That particular electron can have spin-up or spin-down 
as shown in the lower and upper diagram respectively. 
During this stage of the pulse, lasting t wa u, the electron 
is trapped on the dot and Coulomb blockade prevents 
a second electron to be added. After t wa u the voltage 
pulse is reduced, in order to position the energy levels 
in the readout configuration. If the electron has spin-up, 
its energy level is below [i res , so the electron remains on 
the dot. If the electron has spin-down, its energy level is 
above fi res , so the electron tunnels to the reservoir after 



a typical time « T -1 . Now Coulomb blockade is lifted 
and an electron with spin- up can tunnel onto the dot. 
Effectively, the spin on the dot has been flipped by a sin- 
gle electron exchange with the reservoir. After t rea d, the 
pulse ends and the dot is emptied again. 

The expected QPC-response, AIqpc, to such a two- 
level pulse is the sum of two contributions (Fig. 115b ). 
First, due to a capacitive coupling between pulse-gate 
and QPC, AIqpc will change proportionally to the pulse 
amplitude. Second, AIqpc tracks the charge on the dot, 
i.e. it goes up whenever an electron tunnels off the dot, 
and it goes down by the same amount when an elec- 
tron tunnels onto the dot. Therefore, if the dot contains 
a spin-down electron at the start of the readout stage, 
a characteristic step appears in AIqpc during t rea d for 
spin-down (dotted trace inside grey circle). In contrast, 
AIqpc is flat during t rea d for a spin- up electron. Mea- 
suring whether a step is present or absent during the 
readout stage constitutes the spin measurement. 

Fig. [T6b shows experimental traces of the pulse- 
response at an in-plane field of 10 T. The expected two 
types of traces are indeed observed, corresponding to 
spin-up electrons (as in the top panel of Fig. [TBb ) , and 
spin-down electrons (as in the bottom panel of Fig. [T6k). 
The spin state is assigned as follows: if AIqpc crosses 
a threshold value (grey line in Fig. \Wk ) , the electron is 
declared 'spin-down'; otherwise it is declared 'spin-up'. 

As t wa it is increased, the number of 'spin-down' traces 
decays exponentially (see Fig. fTHb) , precisely as expected 
because of spin relaxation to the ground state. This con- 
firms the validity of the spin readout procedure. The spin 
decay time T± is plotted as a function of B in the inset 
of Fig. 116b . The processes underlying the spin relaxation 
will be discussed in section IViTl 

The fidelity of the spin measurement is characterized 
by two error probabilities a and (3 (see inset to Fig. [16b) . 
Starting with a spin-up electron, there is a probability 
a that the measurement yields the wrong outcome 'J.'. 
Similarly, f3 is the probability that a spin-down electron 
is mistakenly measured as 'f- These error probabili- 
ties ca n be determined from complementary measure- 
ments (|Elzerman et all l2004bl) . Both a and (3 depend 
on the value of the threshold as shown in Fig. [T6b for 
data taken at 10 T. The optimal value of the threshold 
is the one for which the visibility 1 — a — (3 is maximal 
(vertical line in Fig. [T6b). For this setting, a=0.07 and 
/?=0.28, so the measurement fidelity for the spin-up and 
the spin-down state is 0.93 and 0.72 respectively. The 
measurement visibility in a single-shot measurement is 
thus 65%, and the fidelity (1 - (a + /3)/2) is 82%. Signifi- 
cant improvements in the spin measurement visibility can 
be made by lowering the electron temperature (smaller a) 
and by making the charge measurement faster (smaller 
/?)■ 

The first all-electrical single-shot readout of an electron 
spin has thus been performed using E-RO. However, this 
scheme has a few drawbacks: (i) E-RO requires an energy 
splitting of the spin states larger than the thermal energy 
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FIG. 16 Experimental results of single-shot readout of a single electron spin, (a) The two types of measurement outcomes, 
corresponding to a spin-up electron (upper panel) and spin-down electron (lower panel); see Fig. 115b for comparison, (b) 
Dependence of the fraction of spin-down electrons on the waiting time, showing a clear exponential decay. Red line is a fit to 
the data. Inset: the spin relaxation time Ti as a function of B. (c) Determination of the readout fidelity. Inset: definition of 
the readout error probabilities a and (3. Main figure: experimentally determi ned error probab i lities a t B—10 T. At the vertical 
line, the visibility l-a-/3 reaches a maximum of 65%. Data reproduced from lElzerman et al\ (|2004bh . 



of the electrons in the reservoir. Thus, for a single spin 
the readout is only effective at very low electron tem- 
perature and high magnetic fields (fceT <C AEz). Also, 
interesting effects occurring close to degeneracy, e.g. near 
the singlet-triplet crossing for two electrons, can not be 
probed, (ii) Since the E-RO relies on precise position- 
ing of the spin levels with respect to the reservoir, it is 
very sensitive to fluctuations in the e lectrostatic poten- 
tial. Background charge fluctuations ( Jung et all . l2004f ) 
can easily push the levels out of the readout configura- 
tion, (iii) High-frequency noise can spoil the E-RO by 
inducing photon-assis ted tunneling from the spin ground 
state to the reservoir ( Onac et aZ.L l2006h . Since the QPC 
is a source of shot noise, this limits the current through 
the QPC and thereby the band width of the charge detec- 
tion ( Vandersvpen et all , 12004 ) . These constraints have 
motivated the search for a different method for spin-to- 
charge conversion, and have led to the demonstration 
of the tunnel-rate-selective readout (TR-RO) which we 
treat in the next section. 



C. Single-shot spin readout using a difference in tunnel rate 

The main ingredient necessary for TR-RO is a spin de- 
pendence in the tunnel rates. To date, TR-RO has only 
been demonstrated for a two-electron dot, where the elec- 
trons are either in the spin-singlet ground state, denoted 



by IS*), or in a spin-triplet state, denoted by |T). In \S), 
the two electrons both occupy the lowest orbital, but 
in |T) one electron is in the first excited orbital. Since 
the wave function in this excited orbital has more w eight 
near the edge of the dot ( Kouwenhoven et aLl . l200lh . the 
coupling to the reservoir is stronger than for the lowest 
orbital. Therefore, the tunnel rate from a triplet state to 
the reservoir Ft is much large r than the rate from the 
singlet state T s , i.e. T T > T s (|Hanson et all l2004b l ) . 



The TR-RO is tested experimentally in lHanson et al\ 
(2005) by applying gate voltage pulses as depicted in 
Fig. [TTa . Figure [TTb shows the expected response of 
Iqpc to the pulse, and Fig.[T7c depicts the level diagrams 
in the three different stages. Before the pulse starts, there 
is one electron on the dot. Then, the pulse pulls the levels 
down so that a second electron can tunnel onto the dot 
(N = 1 — >2), forming either a singlet or a triplet state with 
the first electron. The probability that a triplet state is 
formed is given by 3Tt/(Ts + 3Tt), where the factor of 
3 is due to the degeneracy of the triplets. After a vari- 
able waiting time t wa a the pulse ends and the readout 
process is initiated, during which one electron can leave 
the dot again. The rate for tunneling off depends on the 
two-electron state, resulting in the desired spin-to-charge 
conversion. Due to the direct capacitive coupling of the 
pulse gate to the QPC channel, AIqpc follows the pulse 
shape. Tunneling of an electron on or off the dot gives 
an additional step in AIqpc &s indicated by the arrows 
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FIG. 17 Single-shot readout of two-electron spin states using 
TR-RO. (a) Voltage pulse waveform applied to one of the gate 
electrodes, (b) Response of the QPC current to the waveform 
of (a), (c) Energy diagrams indicating the positions of the 
levels during the three stages. In the final stage, spin is con- 
verted to charge information due to the difference in tunnel 
rates for states \ S) and \T). 



in Fig. [T3>. 

In the experiment, T$ is tuned to 2.5 kHz, and Tt is 
« 50 kHz. The filter bandwidth is 20 kHz, and there- 
fore many of the tunnel events from \T) are not re- 
solved, but the tunneling from \S) is clearly visible. Fig- 
ure [T8a shows several traces of AIqpc, from the last part 
(0.3 ms) of the pulse to the end of the readout stage (see 
inset), for a waiting time of 0.8 ms. In some traces, there 
are clear steps in AIqpc, due to an electron tunneling 
off the dot. In other traces, the tunneling occurs faster 
than the filter bandwidth. In order to discriminate be- 
tween \S) and \T), the number of electrons on the dot 
is determined at the readout time (vertical dashed line 
in Fig. fl"8a ) by comparing AIqpc to a threshold value 
(as indicated by the horizontal dashed line in the bottom 
trace of Fig. [T5kV If AIqpc is below the threshold, it 
means N = 2 and the state is declared 'S'. If AIqpc is 
above the threshold, it follows that N = 1 and the state 
is declared 'T'. 

To verify that 'T' and 'S' indeed correspond to the spin 
states |T) and \ S), the relative occupation probabilities 
are changed by varying the waiting time. As shown in 
Fig. [18b . the fraction of 'T' indeed decays exponentially 
as t wa it is increased, due to relaxation, as before. The er- 
ror probabilities are found to be a = 0.15 and 0= 0.04, 
where a {(3) is the probability that a measurement on the 
state \S) (\T) yields the wrong outcome 'T' ('S'), The 
single-shot visibility is thus 81% and the fidelity is 90%. 
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FIG. 18 (a) Real-time traces of AIqpc during the last part 
of the waveform (dashed box in the inset), for t wa it = 0.8 ms. 
At the vertical dashed line, N is determined by comparison 
with a threshold (horizontal dashed line in bottom trace) and 
the spin state is declared 'T' or ' S' accordingly, (b) Fraction 
of 'T' as a function of waiting time at Bn = 0.02 T, showing a 
single-exponential decay with a time constant T\ of 2.58 ms. 
(c) Normalized fraction of 'T' vs. twait for different values of 
B/i. The singlet-triplet splitti ng Est in thi s expe riment is 
1 meV. Data reproduced from [Hanson et al\ (|2005T ). 



These numbers agree very well wit h the values predicte d 
by a simple rate-equation model ( Hanson et all l2005h . 
Figure [TBb shows data at different values of the mag- 
netic field. These results are discussed in more detail in 
section I VIII 

A major advantage of the TR-RO scheme is that it 
does not rely on a large energy splitting between the 
spin states. Furthermore, it is robust against back- 
ground charge fluctuations, since these cause only a 
small variat i on in the tunnel rates (of order 10~ 3 in Ref. 
( Jung et ~al\ , \2004) ) . Finally, photon-assisted tunneling is 
not harmful since here tunneling is energetically allowed 
regardless of the initial spin state. Thus, TR-RO over- 
comes several constraints of E-RO. However, TR-RO can 
only be used if there exist state-dependent tunnel rates. 
In general, the best choice of readout method will depend 
on the specific demands of the experiment and the nature 
of the states involved. 



It is interesting to think about a measurement proto- 
col that would leave the spin state unaffected, a so-called 
Quantum Non-Demolition (QND) measurement. With 
readout schemes that make use of tunneling to a reser- 
voir as the ones described in this section, QND mea- 
surements are not possible because the electron is lost 
after tunneling; the best one can do in this case is to re- 
initialize the dot electrons to the state they were in before 



20 



the tunneling occured ( Meunier et all feOQGf ) . However, 
by making th e electron tunnel not to a reservoir , but t o 
a second dot (|Engel et all 12004 lEngel and Losl l2005h . 
the electron can be preserved and QND measurements 
are in principle possible. One important example of such 
a scheme is the readout of double-dot singlet and triplet 
states using Pauli blockade that we will discuss in Sec- 
tion Euro 



VII. SPIN-INTERACTION WITH THE ENVIRONMENT 

The magnetic moment of a single electron spin, 
/xb=9.27x 10~ 24 J/T, is very small. As a result, electron 
spin states are only weakly perturbed by their magnetic 
environment. Electric fields affect spins only indirectly, 
so generally spin states are only weakly influenced by 
their electric environment as well. One notable excep- 
tion is the case of two-electron spin states - since the 
singlet-triplet splitting directly depends on the Coulomb 
repulsion between the two electrons, it is very sensiti ve 
to electric field fluctuations (|Hu and Das Sarmal l2006h - 
but we will not discuss this further here. 

For electron spins in semiconductor quantum dots, the 
most important interactions with the environment occur 
via the spin-orbit coupling, the hyperfine coupling with 
the nuclear spins of the host material and virtual ex- 
change processes with electrons in the reservoirs. This 
last process can be efficiently suppressed by reducing the 
dot-reservoir tunnel coupling or creating a large gap be- 
twee n the eletrochemical po tentials in the dot and in the 
lead ( Fuiisawa et aZ.l . r2002af ). and we will not further con- 
sider it in this section. The effect of the spin-orbit and 
hyperfine interactions can be observed in several ways. 
First, the spin eigenstates are redefined and the energy 
splittings are renormalized. A good example is the fact 
that the (/-factor of electrons in bulk semiconductors can 
be very different from 2, due to the spin-orbit interaction. 
In bulk GaAs, for instance, the ^-factor is —0.44. Sec- 
ond, fluctuations in the environment can lead to phase 
randomization of the electron spin, by convention charac- 
terized by a time scale T2 . Finally, electron spins can also 
be flipped by fluctuations in the environment, thereby 
exchanging energy with degrees of freedom in the envi- 
ronment. This process is characterized by a timescale 



A. Spin-orbit interaction 

1. Origin 

The spin of an electron moving in an electric field 
E experiences an internal magnetic field, proportional 
to E X p, where p is the momentum of the electron. 
This is the case, for instance, for an electron "orbiting" 
about a positively charged nucleus. This internal mag- 
netic field acting on the spin depends on the orbital the 



electron occupies, i.e. spin and orbit are coupled. An 
electron moving through a solid also experiences electric 
fields, from the charged atoms in the lattice. In crys- 
tals that exhibit bulk inversion asymmetry (BIA), such 
as in the zinc-blende structure of GaAs, the local elec- 
tric fields lead to a net contribution to the spin-orbit in- 
teraction (which generally becomes stronger for heavier 
elements). This effect is known as th e Dresselhaus con- 
tribution to the spin-orbit interact i onrtDresselhaus . Il955t 
iDvakonov and KachorovskiiL Il986t I Wrinkled . 120031 1. 



In addition, electric fields associated with asymmetric 
confining potentials also give rise to a spin-orbit interac- 
tion (SIA or structural inversion asymmetry). This oc- 
curs for instance in a 2DEG formed at a GaAs/AlGaAs 
heterointerface. It is at first sight surprising that there 
is a net spin-orbit interaction: since the state is bound 
along the growth direction, the average electric field in 
the conduction band must be zero (up to a correction due 
to the effective mass discontinuity at the interface, which 
results in a small force that is balanced by a small average 
electric field). The origin of the net spin-orbit interaction 
lies in mixing with other bands, mainly the valen ce band , 
which contribu t e a no n-zero average electric field ( Pfefferl . 
Il999l : IWrinklerl [2003h . Only in symmetric quantum wells 
with symmetric doping, these other contributions are 
zero as well. The spin-orb it contribution from SIA is 
known as the Rashba term ( Bvchkov and Rashbal 11984 
lRashbal . fl960h . 



2. Spin-orbit interaction in bulk and 2D 

In order to get insight in the effect of the Dresselhaus 
spin-orbit interaction in zin c-blende crystals, we start 
from the bulk Hamiltonian ( Dvakonov and Perell 1 1972t 
IWrinklerL[200l . 



W D oc [p x (p - p z )a x +p y {p z - p x )(T y +p z (p x 



(11) 

where x, y and z point along the main crystallographic 
directions, (100), (010) and (001). 

In order to obtain the spin-orbit Hamiltonian in 2D 
systems, we integrate over the growth direction. For 
2DEGs grown along the (001) direction, (p z ) — 0, and 
(p 2 z ) is a heterostructure dependent but fixed number. 
The Dresselhaus Hamiltonian then reduces to 
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The first two terms are the linear Dresselhaus terms and 
the last two are the cubic terms. Usually the cubic terms 
are much smaller than the linear terms, since (p 2 ) 3> 
p x ,p 2 due to the st rong c onfinement along z. We then 
retain (jDresselhausl . Il955| ) 
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where (3 depends on material properties and on (p 2 ). It 
follows from Eq. [13] that the internal magnetic field is 
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aligned with the momentum for motion along (010), but 
is opposite to the momentum for motion along (100) (see 
Fig. USt). 






Py 







FIG. 19 The small arrows indicate the orientation of the ap- 
parent magnetic field acting on the electron spin as a result 
of (a) the Dresselhaus and (b) the Rashba spin-orbit interac- 
tion when the electron travels through a GaAs crystal with 
momentum p. 

Similarly, we now write down the spin-orbit Hamilto- 
nian for the Rashba contribution. Assuming that the 
confining electric field is along the z-axis, we have 

Hr oc [E x p] a = E z (-p y a x + p x <r y ) , (14) 



or 



'Hr = ot{-pyO x + p x a y 



(15) 



with a a number that is material-specific and also de- 
pends on the confining potential. Here the internal mag- 
netic field is always orthogonal to the momentum (see 
Fig. H9t>). 

We point out that as an electron moves ballistically 
over some distance, the angle by which the spin is 
rotated, whether through Rashba or linear Dresselhaus 
spin-orbit interaction, is independent of the velocity with 
which the electron moves. The faster the electron moves, 
the faster the spin rotates, but the faster the electron 
travels over the distance I as well. In the end, the rota- 
tion angle is determined by I and the spin-orbit strength 
only. A useful quantity is the distance associated with 
a 7r rotation, known as the spin-orbit length, lso- I n 
GaAs, estimates for (3 vary from 10 3 m/s to 3 x 10 3 m/s, 
and it follows that the spin-orbit length, lso — fi/{/3m*) 
is 1 — lOjum, in agreement with e xperimentally mea- 
sured values ( Zumbuhl et all 120021 ). The Rashba con- 
tribution can be smaller or larger than the Dresselhaus 
contribution, depending on the structure. From Fig. 1191 
we see that the Rashba and Dresselhaus contributions 
add up for motion along the (110) direction and oppose 
each other along (110), i.e. the spi n-orbit interaction is 
anisotropic (jKonemann et all , 12005). 

In 2DEGs, spin-orbit coupling (whether Rashba or 
Dressel haus) can lead to sp in relaxation via several mech- 
anisms (IZutic et all 2004h. The D'yakonov-Perel mech- 
anism ([Dvakonov and PerelLll972tlWrinkler[|2003f) refers 
to spin randomization that occurs when the electron 



follows randomly oriented ballistic trajectories between 
scattering events (for each trajectory, the internal mag- 
netic field is differently oriented). In addition, spins can 
be fli pped upon scattering, via the Elliot- Yafet mecha- 
nism (lElliottj. 11954 lYafedll963l) or the Bir-Aronov-Pikus 



mechanism (jBir et all 



3. Spin-orbit interaction in quantum dots 

From the semi-classical picture of the spin-orbit inter- 
action, we expect that in 2D quantum dots with dimen- 
sions much smaller than the spin-orbit length lso, the 
electron spin states will be hardly affected by the spin- 
orbit interaction. We now show that the same result 
follows from the quantum-mechanical description, where 
the spin-orbit coupling can be treated as a small pertur- 
bation to the discrete orbital energy level spectrum in 
the quantum dot. 

First, we note that stationary states in a quantum dot 
are bound states, for which (p x ) = (p y ) = 0. This leads 
to the important result that 

(nl i \H so \nl |) cx (nl\ Px . y \nl)(i \a x , y \1) = , (16) 

where n and I label the orbitals in the quantum 
dot, and Hso stands for the spin-orbit Hamiltonian, 
which consists of terms of the form p X:V a x ^ y both 
for the Dresselhaus and Rashba contributions. Thus, 
the spin-orbit interaction does not directly couple the 
Zeeman-split sublevels of a quantum dot orbital. How- 
ever, the spin-orbit Hamiltonian does couple states 
that contain both different orbita l and different spin 
parts! Khaetskii and Nazarovi , l2000f) . As a result, what 
we usually call the electron spin states 'spin-up' and 
'spin-down' in a quantu m dot, are in reality admixture s 
of spin and orbital states ([Khaetskii and Nazarovi [20011) . 
When the Zeeman splitting is well below the orbital level 
spacing, the perturbed eigenstates can be approximated 
as 



Mt) (1) = MT) 



■ E 

n'l'^nl 



E n i — En'v — AE? 



I»U> (1) = Mi) 



n'V ^nl 
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(the tru e eigenstates can b e obtained via exact diagonal- 
ization ( Cheng et aLl . l2004l )). Here AEz refers to the un- 
perturbed spin splitting (in the remainder of the review, 
AEz refers to the actual spin splitting, including all per- 
turbations). The energy splitting between the spin-up 
and spin-down states will be renormalized accordingly, 
AE^ ] = e\ 1] - e\ 1] (|d"e Sousa and Das Sarmal . 12003d: 
IStano and Fabianl . 120051) (see also Fig. H^a)). In GaAs 
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few-electron quantum dots, the measured g-factor in ab- 
solute value is usually in the range of 0.2 — 0.4, and is 
sometimes magnetic field dependent (see Fig. [7]) . A sim- 
ilar beh avior of the q- f actor was found in GaAs/AlGaAs 
2DEGs (|Dobers et aLl . Il988h . 

In contrast to single-electron spin states in a quan- 
tum dot, the lowest two-electron spin states, sin- 
glet and triplet, are coupled directly by the spin- 
orbit interaction (except for To and S, which are 
not coupled to lowest order in t he spin-orbit interac- 
tion, due to spin selection rules (Tciimente et 



2007 



Dickmann and Hawrvlak, 2003; Florescu and Hawrylak, 
20061 : iGolovach et all 120071 ; ISasaki et all 120061) ). This 



is not so surprising since the singlet and triplet states 
themselves involve different orbitals. Nevertheless, cou- 
pling to two-electron spin states composed of higher or- 
bitals needs to be i ncluded as well, as their effect is gener- 
ally n ot negligible ( Climente et a/.l . [2007tlGolovach et all 
120071) . The leading order correction to the two-electron 
wavefunction is then given by 



ISq)^ = \Sq) 



{ T ±q '\H so \ Sq ) lT±q/) 

77 E T ±q> - £j S q 
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q^q E T q'-^T ± q 

where g is shorthand for the quantum numbers n\l\nili 
that label the orbital for each of the two electrons. It can 
be seen from inspection of the spin-orbit Hamiltonian and 
the form of the wavefunctions that many of the matrix 
elements in these expressions are zero. A more detailed 
discussio n is beyond the scope of this review but can be 
found in lClimente et all (|2007fk IGolovach et all (|2007h . 



4. Relaxation via the phonon bath 

Electric fields cannot cause transitions between pure 
spin states. However, we have seen that the spin-orbit 
interaction perturbs the spin states and the eigenstates 
become admixtures of spin and orbital states, see Eqs.[T7l 
|2T1 These new eigenstates can be coupled by electric 
fields (see Fig. [2H|) . and electric field fluctuatio n s can 
lead to spin relaxation ([Khaetskii and Nazarovl l2000l 
l200lt IWoods et al\ . 120021 ). As we will see, this indirect 
mechanism is not very efficient, and accordingly, very 
long spin relaxation times have been observed experi 



mentally (lAmasha et all, l2006t Elzerman et al. 



iFuiisawa et all l2002al l2001bt lHanson et all \20oA 12003 ; 



2004b; 



Kroutvar et al. I, 12004 iMeunier et all , 120071 : ISasaki et all , 
20061) . 

In general, fluctuating electric fields could arise 
from many sources, including fluctuations in the 

gate potentials, background c harge fluctuati ons or 

other electrical noise sources dBorhani et all , l2006t 



iMarquardt and Abalmasso\il2005l ). However, as we shall 
see, it appears that in carefully designed measurement 
systems, the electric field fluctuations of these extrane- 
ous noise sources is less important than those caused by 
the phonon bath. Phonons can produce electric field fluc- 
tuations in two ways. First, so-called deformation poten- 
tial phonons inhomogeneously deform the crystal lattice, 
thereby altering the bandgap in space, which gives rise to 
fluctuating electric fields. This mechanism occurs in all 
semiconductors. Second, in polar crystals such as GaAs, 
also homogeneous strain leads to electric fields, through 
the piezo-electric effect (piezo-electric phonons). 

The phonon-induced transition rate between the renor- 
malized states |n, Z,T) and |n,/,J.) ' is given by 
Fermi's golden rule (an analogous expression can be de- 
rived for relaxation from triplet to singlet states, or be- 
tween other spin states): 



2tt 



r = — \ (1) (nl T \H e , ph \nl i) {1) \ 2 D(AE 



(22) 



where D(E) is the phonon density of states at energy E. 
H e p h is the electron-phonon coupling Hamiltonian, given 
by' 



n 



(23) 



where M^j is a measure of the electric field strength of 
a phonon with wavevector <f and phonon branch j (one 
longitudinal and two transverse modes), r is the position 
vector of the electron and btj and b%j are the phonon 
creation and annihilation operators respectively. 

The relaxation rate thus depends on the phonon den- 
sity of states at the spin-flip energy (the phonons have to 
carry away the energy) , and on how strongly the electron- 
phonon coupling connects the spin-orbit perturbed spin 
states (see Eq. |2"2"|) . The latter in turn depends on (i) 
the degree of admixing between spin and orbital states 
(see Eqs. \T7[ [T8|). (ii) the electric field strength of a sin- 
gle phonon (oc Mq>j), (iii) the effectiveness of phonons 
at coupling different dot orbitals, via e iqr , and (iv) the 
phonon occupation, via (bt. + bq-j). In addition, (v), we 
will see that an external magnetic field is necessary for 
spin relaxation to occur. We next discuss each of these 
elements separately, starting with the density of states. 

The phonons are taken to be bulk phonons in most dis- 
cussions of relaxation in GaAs quantum dots. This may 
not be fully accurate since the dot is formed at a hetero- 
interface, and furthermore is in close vicinity to the sur- 
face (e.g. surface acoustic waves have also been observed 
to couple to dots, when explicitly excited ( Naber et all 
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FIG. 20 Two views on spin relaxation due to the spin-orbit 
interaction and the phonon bath, (a) The electron phonon 
interaction doesn't couple pure spin states, but it does couple 
the spin-orbit perturbed spin states, labeled with superscripts 
1 '. Here TLso and H e ph axe treated sequentially (although 
they don't commute), (b) The combined electron-phonon and 
spin-orbit Hamiltonian couples pure spin states. 



2006)). Nevertheless, this simplification has worked rea- 
sonably well so far for explaining observations of relax- 
ation in dots (lAmasha etal\, 120061: iFuiisawa et ctZI . 11998k 



iKroutvar et all |2004 iMeunier et all 120071) . In what fol- 
lows, we therefore consider bulk phonons only. Further- 
more, we only include acoustic phonons, as the energies 
for optical p honons are much highe r than typical spin- 
flip energies ( Ashcroft and Merminl . I1974T) . Since bulk 
acoustic phonons have a linear dispers ion relation at low 
energies ( Ashcroft and Merminl . Il974l ) , the phonon den- 
sity of states increases quadratically with energy. 

(i) The degree of admixing between spin and orbitals 
obviously scales with the spin-orbit coupling parameters, 
a and /?. Since the spin-orbit interaction is anisotropic 
(a and j3 can add up or cancel out depending on the 
magnetic field orientation with respect to the crystal 
axis), the admixin g and hence the relaxation rate are 
anisotropic as well (Fal ko et aU 12005). Furthermore, the 
admixing depends on how close together in energy the 
relevant orbitals are (see Eqs. [17]- [21]). At an avoided 
crossing of two levels caused by the s pin-orbit interac- 
tion the admixing will be complet e (|Bulaev and Lossl . 
120051: IStano and Fabianl . l2005l l2006h . 

(ii) The electric field associated with a single phonon 
scales as 1/ ^fq for piezo-electric phonons and as yfq for 
deformation potential phonons, where q is the phonon 
wavenumber. This difference can be understood from 
the fact that small phonon energies correspond to long 
wavelengths, and therefore nearly homogeneous crystal 
strain, which can only create electric fields through the 
piezo-electric effect. At sufficiently small energies (below 
~ 0.6 meV in GaAs), the effect of piezo-electric phonons 
thus dominates over the effect of deformation potential 
phonons. As the phonon energy increases, deformation 
potential phonons become more important than piezo- 
electric phonons. 

(iii) How effectively different orbitals are coupled by 
phonons, i.e. the size of the matrix element (nl 
\e %qr \n'V t), depends on the phonon wavelength and the 



dot size ( Bockelmannl . Il994) (this matrix element is ob- 
tained when substituting Eqs. [T7]- HH into Eq. [22]). In 
GaAs, the speed of sound c p h is of the order of 4000 
m/s, so the phonon wavelength is hc p h/ ' E p h, which gives 
~ 16 nm for a 1 meV phonon. For phonon wavelengths 
much shorter than the dot size (phonon energies much 
larger than a few hundred /ieV) , the electron-phonon in- 
teraction is averaged away (the matrix element vanishes). 
Also for phonon wavelengths much longer than the dot 
size, the electron-phonon coupling becomes inefficient, 
as it just shifts the entire dot potential uniformly up 
and down, and no longer couples different dot orbitals 
to each other (this is the regime where the often-used 
dipole approximation applies, where the matrix clement 
(nl t \e %qr \n'l' |) is taken to be cx q). When the phonon 
wavelength is comparable to the dot size, the phonons 
can most efficie ntly couple the orbitals , and spin relax- 
ation is fastest dBulaev and Lossl 120051 : iGolovach et all 
120041 : 1 Woods et al\ l2002h . This role of the phonon wave- 
length (convoluted with the other effects discussed in 
this section) has been clearly observed experimentally, 
see Fig. ED 




FIG. 21 The relaxation rate from two-electron triplet to sin- 
glet states, as a function of the singlet-triplet energy splitting 
(measured with TR-RO, see section IVI.C[) . The relaxation 
rate shows a maximum when the wavelength of the phonons 
with the right energy matches the size of the dot. The en- 
ergy splitting was varied by a magnetic field with a compo- 
nent perpendicular to the 2DEG (circles) and via the gate 
voltages that control the dot potential landscape (triangles). 
We note that the relaxation rate goes down near the singlet- 
triplet crossing (because of the long phonon wavelength and 
vanishing phonon density of states), even though the spin- 
orbit admixing o f singlet and tr i plet i s maximum here. Data 
reproduced from iMeunier et all (|2007l ) . 



(iv) A finite phonon occupation N p h leads to stimu- 
lated emission. This is accounted for by multiplying the 
relaxation rate by a factor (1 + N p h). N p h is given by the 
Bose-Einstein distribution, and can be approximated by 
k B T/E z when k B T > E z . 

(v) The last necessary ingredient for spin-orbit in- 
duced spin relaxation is a finite Zeeman splitting. With- 
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out Zeeman splitting, the various terms that are ob- 
tained when expanding Eg. |2"2"1 using Eqs. [T7J- [TH can- 
cel out (Khactskii and Nazarov , l200lh (this is known as 
van Vleck cancellation; it is a consequence of Kramer's 
theorem). A similar cancellation occurs for spin states 
of two or more electrons. We can understand the need 
for a magnetic field intuitively from the semiclassical dis- 
cussion of the spin-orbit interaction in Section IVII.A.21 
A phonon produces an electric field that oscillates along 
a certain axis, and this electric field will cause an elec- 
tron in a quantum dot to oscillate along the same axis. 
In the absence of any other terms in the Hamiltonian 
acting on the electron spin, the spin-orbit induced ro- 
tation that takes place during half a cycle of the elec- 
tric field oscillation will be reversed in the next half 
cycle, so no net spin rotation takes place. This is di- 
rectly connected to the fact that the spin-orbit interac- 
tion obeys time-reversal symmetry. In contrast, in the 
presence of an external magnetic field, the spin rotation 
(about the sum of the external and spin-orbit induced 
magnetic field) during the first half period doesn't com- 
mute with the spin rotation during the second half pe- 
riod, so that a net spin rotation results. Theory predicts 
that this e ffect leads to a B 2 dependence of the relaxation 
rate 1/Ti dBulaev and Lossl 120051: T Golovach et ~a!\ , l2007t 
IKhaetskii and Nazarovl 200lh ! A clear Bo dependence 
was indeed seen experimentally, see Fig. 

When two phonons are involved, a net spin rotation 
can be obtained even at zero field. Here, the electron will 
in general not just oscillate back and forth along one line, 
but instead describe a closed trajectory in two dimen- 
sions. Since the spin rotations induced during the various 
legs along this t rajectory generally do n't commute, a net 
rotation results ( San- Jose et aLll2006f ). Such two-phonon 
relaxation processes become relatively more imporant at 
very low magnetic fields, where single-phonon relaxation 
becomes very i neffic ient (jKhaetskii and Nazarovl l200lt 
ISan-Jose et oLI . l2006h (Fig. [22]). 

Putting all these elements together, we can predict the 
Bq dependence of the relaxation rate l/Ti beween Zee- 
man split sublevels of a single electron as follows. First, 
the phonon density of states increases with AE\. Next, 
the electric field amplitude from a single phonon scales 
as y/q oc y/AEz for deformation potential phonons and 
as l/yfq oc \j\fAEz for piezo-electric phonons. Fur- 
thermore, for Bq up to a few Tesla, AEz is well be- 
low the cross-over point where the dot size matches 
the phonon energy (several 100 fJ.eV), so we are in the 
long- wavelength limit, where the matrix element (nl j 
\e iqr \n'l' 1) scales as q oc AEz- Finally, due to the effect 
of the Zeeman splitting, the matrix element in Eq. [22] 
picks up another factor of AEz (assuming only single- 
phonon processes are relevant). Altogether, and taking 
into account that the rate is proportional to the ma- 
trix element squared, T-f" is predicted (at low temper- 
ature) to vary with AE% for coupling to piezo-electric 
phonons ( Khaetskii and Nazarovl l200l[ ) , and as AE Z for 
coupling to deformation potential phonons. At high tem- 
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FIG. 22 (Color in online edition) Relaxation rate between 
the two-electron triplet and singlet states in a single dot, as 
a function of in-plane magnetic field B/ /. The in-plane mag- 
netic field doesn't couple to the orbitals and therefore hardly 
modifies the triplet-singlet energy splitting (AEst ~ 1 meV, 
whereas the Zeeman splitting is only ~ 20/j,eV/T in GaAs 
quantum dots). Nevertheless, and as expected, the experi- 
mentally measured rate 1/Ti at first markedly decreases as 
B// decreases, before saturating as B// approaches zero and 
two-phonon relaxation mechanisms set in. The solid line is 
a second-order polynomial fit to the data. For comparison, 
lines with linear, quadratic, and cubic B// dependences are 
shown. The data are extracted f rom Fig. 118b and are repro- 
duced from (|Hanson et all 120051 ) 



perature, there is an extra factor of AEZ 1 . 

We can similarly work out the 1/Ti dependence on the 
dot size, I, or equivalently, on the orbital level spacing, 
AE or b oc l~ 2 (in single dots, AE or b can only be tuned 
over a small range, but in double dots, the splitting be- 
tween bonding and antibonding orbital s can be modi- 
fied o ver several orders of magnitudes (jWang and Wu, 
2006)). The degree of admixing of spin and orbital states 
by Ti-so contributes a factor I 2 to the rate via the nu- 



merator in Eqs. \I7\ - [TBI and another factor of ~ l A via 
E n i — E n 'i' (the dominant part in the denominator in 
Eqs. [T?l fTS]) . Taking the long- wavelength limit as be- 



fore, \(nl 1 



i'l' t)| 2 contributes a factor I 2 . We thus 



arrive at 1/Ti oc I 8 oc AE~ rb . 

In summary, the relaxation rate from spin down to 
spin up (for an electron in the ground state orbital of a 
quantum dot) scales as 



1/Ti oc 



AE%_ 



at temperatures low compared to AEz/kB, and as 

AE 4 z k B T 



1/Tioc 



(24) 



(25) 



at temperatures much higher than AEz/ks- 

Experimentally measured values for T\ between Zee- 
man sublevels in a one-electron GaAs quantum dot are 
shown in Fig. 1231 The relaxation times range from 120 
/is at 14 T to 170 ms at 1.75 T, about 7 orders of mag- 
nitude longer than the relaxation rate between dot or- 
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bitals ( Fuiisawa et all l2002a[ ). The expected B 5 depen- 
dence of T-f 1 is nicely observed over the applicable mag- 
netic field range. A similar depe ndence was observed i n 
optically measured quantum dots ( Kroutvar et q?J .|2004). 
In that system, t he 1/T temperatu re dependence of T\ 
was also verified ( Heiss et all . 120051 ). There are no sys- 
tematic experimental studies yet of the dependence on 
dot size. 
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FIG. 23 (Color in online edition) Relaxation rate between the 
Zeeman split sublevels of the ground state orbital in a quan- 
tum dot (measured with E- RO, see sectio n lVI.BI) . The square 
data points are taken from lAmasha et al\ (I2006T) ; the round 
datapoints are reproduced from lElzerman et al\ (| 2004bl 'l. The 
fact that the two datasets don't connect is explained by a pos- 
sible difference in orbital spacing, crystal orientation etc. For 
comparison, a solid line with a B 5 dependence is shown. 



So far we have first considered the effect of Hso on 
the eigenstates and then looked at transitions between 
these new eigenstates, induced by H e:P h (Fig. l2"0T a)). We 
point out that it is also possible to calculate the matrix 
element between the unperturbed spin states directly, for 
Hso and H e ^ p h together, for instance as 

{nll\(H SO + H e<ph )\in) 

for Zeeman split states of a single orbital (Fig. I2"07 b)). 

Finally, we remark that whereas at first sight phonons 
cannot flip spins by themselves as there are no spin 
operators in the phonon Hamiltonian, 7i e ,p/t, this is 
not strictly true. Since phonons deform the crys- 
tal lattice, the g-tensor may be modulated, and this 
can in fact lead to electron spin flips directly (when 
phonons modulate only the magnitude of the g-factor 
but not the anisotropy of the g-tensor, the electron spin 
phase gets randomized without energy exchange with the 
bath ( Semenov and Kiml l2004f )). Furthermore, the elec- 
tron spin could flip due to the direct relativistic cou- 
pling of the electron spin to the electric field of the emit- 
ted phonon. However, both mechanisms have been esti- 
mated to be much less efficient than the mechanism via 



admixing of spin and orbitals by the spi n-orbit interac- 
tion (jKhaetskii and Nazarovl I20001 l200lh . 



5. Phase randomization due to the spin-orbit interaction 

We have seen that the phonon bath can induce tran- 
sitions between different spin-orbit admixed spin states, 
and absorb the spin flip energy. Such energy relaxation 
processes (described by a time constant T\) unavoidably 
also lead to the loss of quantum coherence (described by 
a time constant T2). In fact, by definition T2 < 2Ti. 

Remarkably, in leading order in the spin-orbit interac- 
tion, there is no pure phase randomizat ion of the elec- 
tron spin, such that in fact T2 = 2Ti ( Golovach et all 
I2004D . For a magnetic field perpendicular to the plane of 



the 2DEG, this can be understood from the form of the 
spin-orbit Hamiltonian. Both the Dresselhaus contribu- 
tion, Eq. [T31 and the Rashba contribution, Eq. [TSJ only 
contain a x and a y terms. With B along z, these terms 
lead to spin flips but not to pure phase randomization. 
However, this intuitive argument doesn't capture the full 
story: for Bq along x, one would expect the a x term to 
contribute to pure phase randomization, but surprisingly, 
in leading order in the spin-orbit interaction, there is still 
no p ure randomization eve n with an in-plane magnetic 
field dGolovach et alJ . l200l . 



B. Hyperfine interaction 

1. Origin 

The spin of an electron in an atom can interact with 
the spin of "its" atomic nucleus through the hyperfine 
coupling. An electron spin in a quantum dot, in con- 
trast, may interact with many nuclear spins in the host 
material fFig. |2"4")) . The Hamiltonian for the Fermi con- 
tact hyperfine interaction is then given by 



N 



H.RF — 22 AkhS 



(26) 



where 1^ and S are the spin oper ator for nucleus 
k and the electron sp i n respectively ([Abragaml 1961 ; 
Abragam and Bleanevl . 119861 : iMeier and Zakharchenva . 
1984 ISlichterl . U990T ). Since the electron wavefunction is 



inhomogeneous, the coupling strength, A^, between each 
nucleus k and the electron spin varies, as it is propor- 
tional to the overlap squared between the nucleus and 
the electron wavefunction. 

This asymmetric situation combined with fast electron 
spin dynamics and slow nuclear spin dynamics, gives rise 
to a subtle and complex many-body quantum mechanical 
behavior, whereby the nuclear spins affect the electron 
spin time evolution, and the electron spin in turn acts 
back on the dynamics of each of the nuclei. 

Since both the nuclear spins and the localized elec- 
tron spin are quantum objects, the hyperfine coupling 
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FIG. 24 One electron spin interacts with (a) a single nuclear 
spin in an atom, versus (b) many nuclear spins in a semicon- 
ductor quantum dot. 



could in principle create entanglement between them (if 
both the electron spin and the nuclear spins ha d a suf - 
ficiently pure initial state, see iBraunstein etail I 1999h . 
For the electron spin, this interaction with uncontrolled 
degre es of freedom in the environment leads to decoher- 
ence (ICoish and Lossl . l2003 iKhaetskii et all l200l l2003t 
iMerkulov et a/.l . l2002f ). This implies that an electron spin 
starting off in a pure state will evolve to a statistical mix- 
ture of several states, i .e. to one of several states , each 
with some probability ( Nielsen and Chuangl . l2000l ) . 

An alternative and very useful description of the effect 
of the nuclei on the electron spin, is to treat the ensemble 
of nuclear spins as an apparent magnetic field, B^. This 
nuclear field, also known as the Overhauser field, acts on 
the electron spin much like an external magnetic field: 



N 



I k \S = gn B B N S 



(27) 



When this nuclear field assumes a random, unknown 
value, the electron spin will subsequently evolve in a ran- 
dom way and thus end up in a statistical mixture of states 
as well, just like in the quantum mechanical description. 

The semiclassical description of the nuclear spins yields 
an intuitive picture of the electron-nuclear dynamics and 
is sufficient to explain all the experimental observations 
discussed in this review. However, we note that the full 
quantum description is required to analyze the correla- 
tions between the microscopic nuclear spin states and the 
single electron spin state, as e.g. in a study of the entan- 
glement between electron and nuclear spins. 

The magnitude of the nuclear field, Bn = 
J2k Ak Ik/ig^Et) is maximum when all nuclear spins 
ar e fully pola r ized. In GaAs, B^ tmax is about 5 
T (|Paget et all . Il977h . For any given host material, this 
value is independent of the number of nuclei, N, that the 
electron overlaps with - for larger numbers of nuclei, the 
contribution from each nuclear spin to Bn is smaller (the 
typical value for A k is proportional to 1 /N) . 

This is distinctly different in the case of (nearly) unpo- 
larized nuclear spins, for instance nuclear spins in ther- 
modynamic equilibrium under typical experimental con- 
ditions. First there is a small average nuclear polariza- 
tion, oriented along the external magnetic field and with 
an amplitude given by the Boltzman distribution (see 
Appendix [A"]l . In addition, there is a statistical fluctua- 
tion about the average, analogous to the case of N coin 
tosses. For an electron spin interacting with N nuclear 
spins 1/2, the root mean square value of the statisti- 



cal fluctuation will be 2?jy.m ax / 

12002k IMerkulov et all . [20021 1. This quantity has recently 
been measured in various semiconductor quantum dots, 
both optica lly dBraun et all, 120051 ; iDutt et dj,l2005l) an d 
electrically (| Johnson et all , 12005a : iKoppens et aLl . I2005T ) . 
giving values in the range of a few mT, as expected 
since N w 10 6 in these dots. Similar values were ob- 
tained earlier for electrons bound to shallow donors in 
GaAs (jDzhioev et aZj . l2002h . 

A few comments on the importance of the host material 
are in order. First, the value of A k is typically smaller for 
lighter nuclei. Second, for particles in p- like orbitals, such 
as holes in GaAs, the wavefunction has almost no overlap 
with the nuclei (only s orbitals have a finite amplitude 
at the nucleus), so the Fermi contact hyperfine coupling 
constant, Ak, will be very small. Third, if a fraction x 
of the nuclei in the host material has zero nuclear spin, 
BN,max is scaled down with a factor 1 — x. The number 
of nuclei contributing to the statistical fluctuations in 
the nuclear field also scales down by 1 — x. As a result, 
the r.m.s. value of the nuclear field scales with y/l — x. 
While x = in GaAs, a fraction x w 0.95 of the nuclei 
( 28 Si) is non-magnetic in natural silicon, and x ~ 0.99 in 
carbon ( 12 C). Furthermore, both for silicon and carbon, 
purification to nearly 100% zero-spin isotopes is possible, 
so really small nuclear fields can be obtained 

Finally, we point out that since the r.m.s. value of the 
statistically fluctuating hyperfine field scales with l/ViV, 
it is much stronger for electrons localized in dots or bound 
to impurities than for electrons with extended wavefunc- 
tions, for instance in 2DEGs, where the electron wave- 
function overlaps with a very large number of nuclei. This 
is in sharp contrast to the effect of the spin-orbit inter- 
action, which becomes suppressed when the electron is 
confined to dimensions shorter than the spin-orbit length, 
such as in small quantum dots. 



2. Effect of the Overhauser field on the electron spin time 
evolution 

The electron spin will precess about the vector of the 
total magnetic field it experiences, here the vector sum 
of the externally applied magnetic field Bq and the nu- 
clear field Bpj- The longitudinal component of Bn, i.e. 
the component oriented parallel or opposed to Bq, di- 
rectly changes the precession frequency by g^sBN, irre- 
spective of the strength Bo (Fig. |2"5r ). Throughout this 
section, we shall call the longitudinal component B Z N . 



For B Z N = 1 mT, the precession rate is increased by 
about 6 MHz (taking g — —0.44), and the electron spin 
picks up an extra phase of 180° in just 83 ns. The effect 
of the transverse components of the nuclear field, B^ y , 
strongly depends on the strength of B , For B <C B x ^ v , 
the electron spin will precess about an axis very close 
to B X N V (Fig. dSb). For B > B X N V ', in contrast, the 
transverse components of the nuclear field only have a 
small effect: a change in the electron spin precession rate 
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by 



g/is-B^r / '(2Bq), and a tilt of the rotation axis by 



1 T and Bf^ = 1 



arctan(£jv/£o) (Fig.|2Sfc). Taking B 
mT, the precession frequency is shifted by just 3 kHz, 
causing an extra phase of 180° only after 166 ms; the 
rotation axis is then tilted by w 0.06°. For external mag- 
netic fields above say 100 mT, we are therefore mainly 
concerned with the longitudinal nuclear field. 
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FIG. 25 Longitudinal magnetic field fluctuations, Bff, add 
directly to the external field, -Bo, whereas transverse fluctua- 
tions, B^ y , change the total field only in second order when 



If the nuclear field, Bn, were fixed and precisely 
known, it would affect the electron spin dynamics in a 
systematic and known way. In this case, there would be 
no contribution to decoherence. However, the orienta- 
tion and magnitude of the nuclear field change over time. 
First, the hyperfine field or Overhauser field, Bn, will 
change if the local nuclear polarization, ^2 k Ik, changes. 
This can occur for instance through dynamic nuclear po- 
larization. Second, Bm can also change while the net nu- 
clear polarization remains constant. This happens when 
two nuclei with different A k flip-flop with each other, 
such that J2k Aklk changes. 

At any given time, the nuclear field thus assumes a ran- 
dom and unknown value and orientation, and this ran- 
domness in the nuclear field directly leads to a random- 
ness in the electron spin time evolution. During free evo- 
lution, the electron spin will thus pick up a random phase, 
depending on the value of the nuclear field, i.e. the single- 
spin coherence decays. The shape of the decay (exponen- 
tial, power-law, etc.) is determined by the distribution 
of nuclear field values. For a longitudinal nuclear field, 
B Z N that is randomly drawn from a Gaussian distribution 
of nuclear fields with standard deviation yj {(B Z N ) 2 ) (e.g. 
when every nuclear spin had equal probabilities for being 
up or down), the decay would be Gaussian as well , i.e. o f 
the form exp(-£ 2 /(T 2 *) 2 ), where (|Merkulov et all |2002|) 



t: 



(28) 



For y/((B z N ) 2 ) = 1 mT, T 2 * would be as short as 30 ns. 

The timescale T 2 * can be measured as the decay time 
of the electron spin signal during free evolution, aver- 
aged over the nuclear field distribution(Fig. I2"6"hh The 
free evolution of the spin can be measured by tipping the 
spin into the x — y plane, subsequently allowing the spin 
to freely evolve about B (assumed to be along z), and 
recording the magnetization in this plane as a function of 
the free evolution time interval (in NMR, this is known as 



the free induction decay or FID). If the spin can only be 
measured in the ±1 basis, a so-called Ramsey experiment 
should be performed instead. It starts off like an FID but 
the magnetization is rotated back to the z axis after the 
free evolution time interval, so it can be measured along 
z. Averaging each data point in an FID or Ramsey ex- 
periment over a sufficiently long time, is equivalent to 
averaging over a large number of uncorrelated nuclear 
field values (assuming the system is ergodic). Such ex- 
periments have recently been performed (see Section UX]) . 
and gave the expected short timescales for T 2 *. 
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FIG. 26 (a) Amplitude of the x-component of the electron 
spin as a function of time, under free precession about (Bo + 
Bff) for three different values of B%, (dotted and dashed lines). 
Also shown is the average of the three oscillations, which is 
seen to rapidly decay (solid line). If an average was taken over 
many more values of B Z N taken from a Lorentzian distribution 
in the same range, the envelope of the oscillation would decay 
with a single exponent (solid lines). A Gaussian distribution 
would yield a Gaussian decay (see text), (b) Representation of 
the rotating magnetic field, Bi, and the longitudinal nuclear 
field, B^, in a reference frame rotating about z at the same 
rate as Bi. The electron spin will precess about the vector 
sum of these two fields, rather than about the axis defined by 
Si. 

Also in the case of driven evolution, the nuclear field 
will affect the time evolution. In a spin resonance ex- 
periment, for instance, a rotating magnetic field B\ is 
applied with frequency g/isBo/h (on-resonance with the 
Zeeman splitting), and perpendicular to Bq (see also Sec- 
tion HX.Ap . In the usual rotating reference frame, this 
corresponds to a rotation about the Bi-axis. However, a 
longitudinal nuclear field will shift the electron spin res- 
onance frequency, so that B\ is no longer on-resonance 
with the electron spin precession. In the rotating frame, 
the spin will then rotate about the vector sum of B\ and 
B Z N (Fig. [Jib) j which may be a rather different rotation 
than intended. In fact, the nuclear field has been the 
main limitation on the fidelity of spin rotations in recent 
electron spin resonance experiments in a quantum dot 
(see Section |LX7A"|) . 

We have so far focussed on the effect of the nuclear field 
(mainly B^) on the electron spin phase. We now briefly 
turn to electron spin flips caused by the nuclear field 
(mainly by B^j y ). The semi-classical picture says that 
for Bq <C B^ v , the electron spin will rotate about B^ v , 
i.e. it is changed from spin-up to spin-down and back, 
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while for Bq ^> B^ v , this hardly occurs (see Fig.[25f . We 
arrive at a similar conclusion from the quantum mechan- 
ically picture: the hyperfine Hamiltonian, Eql26[ permits 
direct electron-nuclear flip-flops only when the two rele- 
vant electron spin states are very close in energy (since 
a nuclear spin flip can absorb only a small amount of 
energy). This effect has been observed in recent experi- 
ments on two-ele ctron singlet and tripl et states in a dou- 
ble quantum dot ||johnson_e£aD, \2QQM) (Section lVni.D|) . 



There still is another contribution from B^ y to spin 
flips (i.e. to Ti). Since the nuclear field strength and 
orientation depend on the collection of nuclear spins 
that the electron wavefunction overlaps with, Bn de- 
pends on the orbital the electron occupies. As a re- 
sult, like the spin-orbit interaction (see section IVII.AI) . 
the hyperfine interaction also leads to admixing of spin 
and orbital states. Here too, phonons can induce tran- 
sitions between the perturbed spin states, and absorb 
the spin-flip energy (lAbalmassov and Ma rquardt. 2004 



Erlingsson and Nazarovl I2002L 12004 lErlingsson et 
200f|) ■ Whereas the transition amplitude due to the spin- 
orbit interaction vanishes in lowest order at Bo = (see 
discussion in section lVII.A.4[) , this is not the case for hy- 
perfine mediated transitions, so the hyperfine mechanism 
will be relatively more important at low magnetic fields. 



3. Mechanisms and timescales of nuclear field fluctuations 

We have seen that the nuclear field only leads to a 
loss of spin coherence because it is random and unknown 
- if Bn were fixed in time, we could simply determine 
its value and the uncertainty would be removed. We 
here discuss on what timescale the nuclear field actually 
fluctuates. This timescale is denoted by t nuc . 

The study of nuclear dynamics due to internuclear 
and electron-nuclear interactions has a long and rich 



history (lAbragaml 196ll : Abragam and Bleanevl . 1 19861 
iMeier and Zakharchenval Il984 ISlichterl . Il99dh . When 
applied to quantum dots, theory predicts that the two 
most important mechanisms responsible for fluctua- 
tions in the nuclear field are the internuclear magnetic 
dipole - di pole interaction Jde Sousa and Das Sarmal . 
2003a. cj; IWitzel and Rogerio de Sousal . 12005c lYao et all 
20051 ) and the elec tr on-nu cl ear hyperfine interac- 



tion (ICoish and Lossl. 12004 IKhaetskii et all l2002t 



Sh envi and Rogerio de SousaTl2005h . 

The Hamiltonian describing magnetic dipole-dipole in- 
teractions between neighbouring nuclei is of the form 



n 



DD 



E 



It ■ h 



' r ij){Ij ' r ij) 



where /j,q is the permeability of free space, gi 



(29) 
is the 



g-factor of nucleus i, (In the nuclear magneton, and 
flj is the vector connecting the two nuclei. The 
strength of the effective magnetic dipole-dipole interac- 
tion between neighbouring nuclei in GaAs is about (100 



/is)~ 1 (IShulman et all [l958). In strong magnetic fields, 
we can discard the non-secular part of this Hamiltonian, 
and retain 
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for the coupling term between nuclear spins i and j of 
the same species (for coupling between spins of different 
isotopes, only the IfH term survives at high field). Here 
1^ are the nuclear spin raising and lowering operators. 
The If I? terms in Eq. [3D] are responsible for changing 
B^ and B V N . This may occur on the 100/xs timescale. 
The flip-flop terms, If If +IfT^, affect B Z N , but the flip- 
flop rate between nuclei i and i + 1 may be suppressed, 
namely when \Ai — Aj+i | is greater tha n the internuclear 
coupling strength (jDeng and Hd . I20061 ) (as this causes an 
energy mismatch). Thus, when we consider the dipolar 
interaction only, B^ y evolves on a 100 fis timescale, but 
B Z N may evolve more slowly. 

We now turn to the hyperfine interaction. So far 
we only considered its effect on the localized electron 
spin, but naturally the interaction (Eq. I26|) works both 
ways and the nuclei evolve about the electron spin, 
just like the electron spin evolves about the nuclear 
field. The apparent field experienced by the nuclei is 
called the Knight shift, and it has a strength A k ~ 
(TOusT 1 dCoish and Lossl |2004 IKhaetskii et al. I, 120021 
iMerkulov et aLl . l2002h (the iV nuclei with which the elec- 
tron wavefunction overlaps "share" the total coupling 
strength A). When we look at the hyperfine interaction 
Hamiltonian, Eq. 1261 we see that similar to the internu- 
clear dipole-dipole interaction, it contains If S z terms as 
well as flip-flop terms: 



JV 



Hhf = J2 A k(i'kS x + i y k sy + r k s z ) 

k 
N 

= J2 A « + I k ~S- + 2I z k S z ) 12 (31) 



where S ± are the electron spin raising and lowering op- 
erators. The transverse component of the nuclear field, 
B^ v , will evolve due to the I Z S Z terms, on a 10/is 
timescale. B Z N will change on the same timescale only 
near Bq = 0, due to the electron-nuclear flip-flop com- 
ponents in Eq. 1311 At finite B , the energy mismatch 
between the electron and nuclear Zeeman energies sup- 
presses electron-nuclear flip-flops, so here B^ cannot 
change by direct electron- nuclear flip-flops. 

The hyperfine interaction can also affect B Z N indirectly. 
Two virtual electron-nuclear flip-flops (between one nu- 
cleus and the electron and between the electron and an- 
other nu cleus) can together lead to a nuc l ear-nuclear 
flip-fl op ( Shenvi and Rogerio de Sousal . 120051 : lYao et all . 
2005). Such a flip-flop process between two nuclei i 
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and j modifies B Z N whenever A, 7^ Aj. This virtual 
electron-nuclear flip-flop process continues to be effective 
up to much higher Bo than real electron- nuclear flip-flops. 
Eventually it is suppressed at high Bo as well. 

Altogether the dipole-dipole and hyperfine interactions 
are expected to lead to moderate timescales (10-100 /is) 
for B^ v fluctuations. At low -Bo, the timescale for 
B Z N fluctuations is similar, but at high B , B Z N fluctu- 
ations are very slow. Here t nuc is certainly longer than 
10 — 100/iS and perhaps longer than a second. This still 
needs to be confirmed experimentally, but an indication 
that t nuc may indeed be very long is that the decay (due 
to spin diffusion) of nuclear polarization built up locally 
at a quantum dot o r impurity, occurs on a timescale of 
seconds to minute s (jHuttel et all 120041 : iKoppens et all 
120051: iPaeetl Il982ft (see Section IVIII.DI) . 



4. Electron spin decoherence in a fluctuating nuclear field 



In section IVII.B. 21 we saw that we lose our knowledge 
of the electron spin phase after a time T 2 , in case the 
nuclear field orientation and strength are unknown. Now 
suppose that we do know the orientation and strength of 
the nuclear field exactly at time t = , but that the nu- 
clear spin bath subsequently evolves in a random fashion 
on a timescale t nuc , as described in the previous sub- 
section. On what timescale, T 2 , will the phase of the 
electron spin then be randomized? 

It may come as a surprise at first that T 2 is not sim- 
ply the same or even of the same order as t nuc . The 
reason is that T 2 depends not only on the timescale of 
the nuclear field fluctuations (t nuc ), but also on the am- 
plitude and stochastics of the fluctuations. The typical 
amplitude is given by the width of the nuclear field distri- 
bution, which can be expressed in terms of l/Tj* . Exam- 
ples of different stochastic models include Gaussian noise 
and Lorentzian noise, which lead to distinct decoherence 
chara cteristics ([Klauder and Anderson! . Il962t Ide Sousal . 
200G). The actual value of T 2 is difficult to calculate 
exactly, but can be estimated in various regimes to be 
1 - 100/xs. 

The timescale T 2 is also hard to obtain experimen- 
tally. In principle T 2 could be determined by record- 
ing an FID or Ramsey decay, whereby B Z N is reset to 
the same initial value for every datapoint. This may re- 
quire measuring B Z N accurately and quickly, i.e. better 
than the initial unc ertainty in B*L and i n a time much 
shorter than t nun dGiedke et all 120061 : iKlauser et all 
120061: IStepanenko et all feOod TAlternativelv, T 2 could 
be obtained by recording all the datapoints needed to 
construct an FID or Ramsey experiment within a time 
short compared to t nuc . 

Experimentally, it may be much easier to obt ain a spin- 
echo decay time, Tprhm well-known from NMR ( Freeman! 
Il997t IVandersvpen and Chuand . |2004 ). In its simplest 
form, the Halm-echo, the random time evolution that 
takes place during a certain time interval r is reversed 



during a second time interval of the same duration, by 
applying a so-called echo pulse (180° rotation) in between 
the two time intervals. Importantly, this unwinding of 
random dephasing only takes place to the extent that 
the random field causing it is constant for the duration 
of the entire echo sequence. Thus, the slow time evolu- 
tion of the nuclear field implies that the echo will not be 
complete. We call the timescale of the remaining loss of 
phase coherence T ec ^ Q . 

Like T 2 , T ec h is much longer than T 2 but also much 
shorter than t nuc . For example, if the nuclear field fluc- 
tuations had Gaussian noise characteristics, the electron 
spin coherence in a H ahn echo experiment wo uld decay 
as exp[-t 3 /ft,„,.r 2 * 2 )] (|Herzog and Hahil Il956h . Taking 
T 2 = 10 ns and t nuc — 10 s, we would obtain a T 2 of 
10/tzs, much faster than t nuc itself. The nuclear field fluc- 
tuations may not be characterized exactly by Gaussian 
noise, but nevertheless, predictions for T 2 still range 
from lyits to 100/is, with contr ibutions from the internu- 
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2003bl: IWitzel and Rogerio de Sousal . l2005t lYao et al 



20051 ) , the electr o n-nucl ea r hyperfine interac 



tion ( Coish and Lossl 12004 iKhaetskii et al. I, I2002D . 



and indirect nuclear-nuclear interactions, mediated b 



tne n ypernne c o uplin g (|bnenvi ana Kogeno de bousal . 
120051 : lYao et all I2005T) . Also, contrary to the usual 
case, the echo decay is not well described by a single 
exponential. The predicted form for the echo decay 
depends on the magnetic field strength, and on what 
terms in the Ha miltonian are then important, but can 
be very complex ( Coish and Lossl 12004 Ide Sousal . 120061 : 
lYao et q/.l . l2005h . 

The usefulness of the echo technique for effectively 
obtaining an extended coherence time has been demon- 
strated experimentally with two electron spins in a dou- 
ble quantum dot, whereby a lower bound on T e c ^ of 
1/is was o btained at 100 mT (jPetta et all l2005bl ) (Sec- 
tion [DCbJi. Similar echo- like decay times were observed 
in optical measurements on an ensemble o f quantum dots 
that e ach contain a single electron spin ( Greilich et all 
l2006bD . 

We note that sometimes a distinction is made be- 
tween "dephasing" , referring to a loss of phase coher- 
ence that can be reversed with echo techniques, and 
"decoherence", referring to a loss of phase coherence 
that cannot be reversed. While this distinction is use- 
ful in practice, we note that it is also somewhat artifi- 
cial, in the sense that any time evolution can in princi- 
ple be reversed by a sufficie ntly rapid sequence of multi- 
ple generalized echo pulses (lAugustine and Hahnl . Il997t 
IViola et all 119991 : IViola and Llovdl . Il998h . 

Finally, we point out that it may be possible to 
extend t nuc , i.e. to (almost) freeze the nuclear field 
fluctuations. One possibility to do this is to fully 
polarize the nuclear spins: if all nuclear spins point 
the same way, nuclear-nuclear flip-flop processes can 
no longer take place and also electron- nuclear flip-flops 
can only have a very small effect ( Khaetskii et all 
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I2002L l2003t ISchliemann et all l2002h . For this approach 
to be effective, the nucle ar spin polarization mu st be 
really very close to 100% ( Schliemann et aZ.1 . l2002h : just 
90% polarization hardly helps. At present the highest 
nuclear spin polarizations reach ed in quantum dot s 
are 60 %, via optical pumping (jBracker et Ml, l2005h . 
Certainly other mechanisms for freezing the nuclear spin 
fluctuations could be considered, but no such effect has 
been demonstrated to date. 



VIII. SPIN STATES IN DOUBLE QUANTUM DOTS 

In this section, we discuss the spin physics of double 
quantum dots. We start by describing the properties of 
"spinless" electrons. Then, we show how the spin selec- 
tion rules can lead to a blockade in electron transport 
through the double dot. Finally, we describe how this 
spin blockade is influenced by the hyperfine interaction 
with the nuclear spins, and discuss the resulting dynam- 
ics. 



A. Electronic properties of electrons in double dots 



C. Summary of mechanisms and timescales 

Our present understanding of the mechanisms and 
timescales for energy relaxation and phase randomization 
of electron spins in few-electron quantum dots is summa- 
rized as follows (as before, most numbers are specific to 
GaAs dots, but the underlying physics is similar in other 
dot systems). 

Energy relaxation is dominated by direct electron- 
nuclear flip-flops near zero field (or whenever the relevant 
electron spin states are degenerate). In this case, T\ is 
as low as 10-100 ns. As Bo increases, electron-nuclear 
flip-flops become suppressed, and energy must be dissi- 
pated in the phonon bath. Spin-phonon coupling is in- 
efficient, and occurs mostly indirectly, either mediated 
by the hyperfine interaction or by spin-orbit interaction. 
As a result Ti rapidly increases with Bo, and at 1.75 T, 
Ti has been measured to be 170 ms. As Bo further in- 
creases, the phonon density of states increases and the 
phonons couple more efficiently to the dot orbitals (the 
phonon wavelength gets closer to the dot size), so at some 
point relaxation becomes faster again and Ti decreases 
with field. At 14 T, a 120^s Ti has been observed. At 
still higher fields, the phonon wavelength would become 
shorter than the dot size, and Ti is once more expected 
to go up with field. 

Phase coherence is lost on much shorter timescales. A 
rapid dephasing of the electron spin results from the un- 
certainty in the nuclear field, T 2 * w 10 ns, irrespective of 
Bq. If the uncertainty in the nuclear field is removed or if 
the resulting unknown time evolution is unwound, we re- 
cover T 2 or T ec /j respectively, which are much longer. 
Phase randomization of the electron spin then results 
from the (slow) fluctuations in the nuclear field, which 
occur on a timescale of 100 /is to perhaps seconds, and 
should lead to a Tb or T ec ho of 1 — 100/^s. Indeed, a lower 
bound on T ec / 10 of 1/is was experimentally observed at 
100 mT. If the effect of the nuclear field on the electron 
spin coherence could be suppressed, the spin-orbit inter- 
action would limit T 2 , to a value of 2Ti (to first order 
in the spin-orbit interaction), which is as we have seen a 
very long time. 



We first ignore the spin of the electrons and de- 
scribe the basic electronic properties of double quan- 
tum dots. The properties of "spinless" electrons in 
do uble dots are treated in d etail by Van der Wiel et 
al. ( Van der Wiel et q/ll2003f) . Here, we give all the the- 
ory relevant for electron spins in double dots without 
going into the details of the derivations. 



1. Charge stability diagram 
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FIG. 27 Charge stability diagrams for (a) uncoupled and (b) 
coupled double dots, depicting the equilibrium electron num- 
bers (Ni,N2) in dot 1 and 2 respectively. The lines indicate 
the gate voltage values at which the electron number changes. 
In (b), a finite cross-capacitance between gate 1 (2) and dot 
2 (1) is taken into account. 

Consider two quantum dots, labelled 1 and 2, whose 
electrochemical potentials are controlled independently 
by the gate voltages Va,i and Vg,2, respectively. Fig- 
ure [27^, shows the equilibrium electron numbers (Ni,N2) 
of the quantum dots as a function of Vg,i and Vg,2, for 
the case that the dots are completely uncoupled. Such a 
plot is called a charge stability diagram. The lines indi- 
cate the values of the gate voltages at which the number 
of electrons in the ground state changes. Note that the 
lines are exactly horizontal and vertical, since the elec- 
trochemical potential in either dot is independent of the 
charge on the other dot, and each gate voltage only af- 
fects one of the dots. 

When the dots are capacitively coupled, addition of 
an electron on one dot changes the electrostatic energy 
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of the other dot. Also, the gate voltage Vg,i {Yg,2) gen- 
erally has a direct capacitive coupling to quantum dot 2 
(1). The resulting charge stability diagram is sketched in 
Fig. 127b. Each crosspoint is split into two so-called triple 
points. The triple points together form a hexagonal or 
"honeycomb" lattice. At a triple point, three different 
charge states are energetically degenerate. The distance 
between the triple points is set by the capacitance be- 
tween the dots (the interdot capacitance) C m . At low 
source-drain bias voltage, electron transport through the 
double dot is possible only at these triple points. In 
contrast, a charge sensing measurement will detect any 
change in the electron configuration and therefore map 
out all the transitions, including those where an electron 
moves between the dots (e.g. from (0,1) to (1,0)). 
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FIG. 28 (Color in online edition) Charge sensing data on a 
double dot in the few-electron regime. The dark lines signal 
the addition of a single electron to the double dot sytem. 
The absolute number of electrons in dot 1 and 2 is indicated 
in each region as "Ni, Ns" . (a) The absence of dark lines in 
the lower left region indicates that the dot is empty there. 
(b) Zoom-in of the b oxed region of (a). Data adapted from 
lElzerman et al\ <|2003t) . 



Figure[28]shows charge sensing data in the few-electron 
regime. The absence of charge transitions in the lower 
left corner of Fig. I2"5h indicates that here the double dot 
structure is completely depleted of electrons. This allows 
the absolute number of electrons to be determined unam- 
biguously in any region of gate voltage space, by simply 
counting the number of charge transition lines from the 
(0,0) region to the region of interest. Figure [28b displays 
a zoom-in of the boxed region in Fig. I2"8a . The bright 
yellow lines in between the triple points in Fig. I2"5b are 
due to an electron moving from one dot to the other. 
This changes the number of electrons on each individual 
dot, while keeping the total number of electrons on the 
double dot system constant. 

From now on, we assume the dots are in series, such 
that dot 1 is connected to the source and dot 2 to 
the drain reservoir. From a similar analysis as in Sec- 
tion III.D1 the electrochemical potential of dot 1 is found 
to be: 

Mi(JVi,JV 2 ) = U(N U N 2 ) - U(JVi - l,JV a ) 



_|__Cm (Q D y D _|_ C 22 Vg, 2 + C 2 \Vg \) 

e 



(32) 



where Cy is the capacitance between gate j and dot i, 
Cs {Cd) is the capacitance from dot 1 (2) to the source 
(drain) , Eci is the charging energy of the individual dot 
i and Ec m is the electrostatic coupling energy 8 . The 
coupling energy Ec' m is the change in the energy of one 
dot when an electron is added to the other dot. One can 
obtain fi 2 (Ni, N 2 ) by simply interchanging 1 and 2 and 
also C D V D and C S V S in Eq.El 

The solid lines in Fig. [221 depict the electrochemical 
potentials around the triple points for low source-drain 
bias. The diagrams schematically show the level arrange- 
ment at different positions (ni(Ni, N 2 ) and fJ, 2 (Ni, N 2 ) 
are shown in short form (Ni,N 2 ) in the left and right dot 
respectively). In this case, we have assumed the tunnel 
coupling to be small (i.e. negligible with respect to the 
electrostatic coupling energy). This is called the weak- 
coupling regime. 

When the tunnel coupling t c becomes significant, the 
electrons are not fully localized anymore in single dots 
but rather occupy mole c ular orbitals that span both 
dots dVan" der Wiel et all 120031 ) . The molecular bonding 
orbital, ips, and the antibonding orbital, ipA, are super- 
positions of the single-dot states in the left dot, <fii, and 
the right dot, 4> 2 : 



tpB = CUp! + (3<p 2 , 
lp A = P<f>l ~ U(j>2. 



(33) 
(34) 



When the single-dot states are aligned, the energy of 
bonding orbital is lower by \t c \ than the energy of the 
single-dot orbitals, and the energy of the antibonding or- 
bital is higher by the same amount. 

The tunnel coupling is revealed in the charging di- 
agram by a bending of the honeycomb lines near the 
triple points, as depicted by the dashed lines in Fig. [29] 
The value of the tunnel coupling can be determined ex- 
perimentally fr om such a plot by measuring the bend- 
ing o f the lines ( Huttel et all 120051 : iPioro-Ladriere et all 
2005), or by measuring the charge distribution as a 
funct ion of detuning between left and right dot poten- 
tials (|DiCarlo et aZ.IJ2004 IFetta et aLl . l2004 . 

Note that, when drawing the diagram in Fig. [211 we as- 
sumed the electrons to be "spinless" . Therefore, the first 
electron can enter the molecular bonding orbital which 
takes \t c \ less energy than in the case of weak tunnel 
coupling. The second has to move into the antibonding 
orbital because of the Pauli exclusion principle. Thus, 
the extra energy needed to add the second electron is 



Note that in IVan der Wiel et al. I 1120031) the cross-capacitance 
terms (C12 an d C21) are neglected. However, they generally are 
significant in lateral dots. 
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FIG. 29 Electrochemical potential lines around the triple points for weak tunnel coupling (solid and dotted lines) and strong 
tunnel coupling (dashed lines). Each line outlines the gate voltages where the corresponding electrochemical potential in the 
dots is equal to the electrochemical potential in the reservoirs (which is defined as zero). The level diagrams indicate the 
positions of the electrochemical potentials at several points in gate space. Inset: dotted line indicates the "naive" expectation 
for the electrochemical potentials when spin is included. 



Ec + 2 \t c \ {\t c \ more than for the weak-coupling case) 
and the second triple point is pushed to higher gate volt- 
ages with respect to the weak-coupling case. 

When spin is taken into account, the orbitals become 
doubly degenerate and both electrons can occupy the 
bonding orbital. Thus, it only takes Ec extra energy 
to add the second electron. This changes the charging 
diagram drastically; namely, the dashed line in the (1,1) 
region (in the upper right corner in Fig. I29p moves to the 
other side of the triple point! This scenario is sketched 
in the inset of Fig. 1291 However, e xperiments on dou- 
ble dots with large tunnel coupling (Hii ttel et all l2005t 
iPioro-Ladriere etal\ , 120051 ) reproduce the main diagram 
of Fig. [221 and not the diagram of the inset that includes 
spin. The reason for this is that the Coulomb interaction 
is typically one or two orders of magnitude larger than the 
tunnel coupling. Therefore, when the double dot is oc- 
cupied by two electrons, the electrons are again strongly 
localized. The orbital energy of the two-electron system 
is then equal to the sum of the two single-dot orbitals, 
which is the same as the sum of the bonding and the 
antibonding orbital. When the second electron is added, 
the tunnel coupling energy that was gained by the first 
electron has to be "paid back" , which causes the second 
triple point to appear at higher gate voltages than in the 
weak-coupling case. Therefore, Fig.[5H]is recovered when 
spin is included. 



2. High bias regime: bias triangles 

When the source-drain bias voltage is increased, two 
different types of tunneling can occur. Up to now, we 
have only discussed tunneling between aligned levels, 
where the initial and final electronic state by definition 
have the same energy. This is termed elastic tunneling. 
However, tunneling can also occur when there is an en- 
ergy mismatch between the initial and final state (levels 
are misaligned) , in which case the process is called inelas- 
tic. For inelastic tunneling to take place, energy exchange 
with the environment is required to compensate for the 
energy mismatch, since the process as a whole has to con- 
serve energy. One important example of energy exchange 
is the absorption of one or more photons under microwave 
or radio-fr equency radiation, l e ading to photon-assisted 
tunneling ( Van der Wiel et all . f2003) . Energy emission 
usually takes place through phonons in the surrounding 
lattice. Note that at cryogenic temperatures, the number 
of photons and phonons in thermal equilibrium is usually 
negligle. Since inelastic tunneling is a second-order pro- 
cess, the inelastic tunneling rate is in general much lower 
than the clastic tunneling rate. However, when there are 
no aligned levels elastic tunneling is suppressed and in- 
elastic tunneling dominates the electron transport. 

The rate of inelastic tunneling between the dots is 
highly sensitive to the density of states and the occu- 
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pation probability of photons and phonons. Therefore, 
a double dot system can b e used as a probe of the 
semiconductor enviro nment dVan der Wiel et al\. 



or as a nois e detec tor ( Aguado and Kouwenhoverl 
lOnac et Ml, [2006). The energy window that is being 
probed is determined by the misalignment between the 
levels in the two dots. Since this misalignment is eas- 
ily tuned by gate voltages, a wide range of the energy 
spectrum can be investigated with very high resolution 
(typically of order 1 fxeV). 




FIG. 30 The triple points for an applied source-drain bias 
Vsd and the drain kept at ground. A triangle is formed from 
each triple point. Within this "bias triangle" charge transport 
through the dot is energetically allowed. The grey lines and 
regions in the triangles illustrate the gate voltages at which 
transitions involving excited state levels play a role. Electro- 
chemical potentials corresponding to transitions involving an 
excited state are shown in grey in the level diagrams. 

When the source-drain bias voltage is increased, the 
triple points evolve into "bias triangles" , as depicted in 
Fig. [3D] for weak tunnel coupling. The electron numbers 
refer to the triple points where the first electrons are 
added to the double dot system; however, the following 
discussion is valid for any number of electrons on either 
dot. 

To understand the electron transport within such a 
triangle, we first look at the three legs. Along the base 
leg, iii(Ni + l, N 2 )=^2(Ni, N 2 + l), and elastic tunneling 
occurs. Moving along this same slope anywhere in the 
plot will not change the relative alignment of the levels in 
the two dots, but only change their (common) alignment 
with respect to the source and drain. 

Moving upwards along the left leg of the triangle, fii is 



fixed (/ii(l, 0) is aligned with the source electrochemical 
potential) and only /12 is changed. At the bottom of the 
triangle, the levels corresponding to transitions involving 
only the ground states (//i(l,0) and ^(0,1), the black 
levels in the diagrams) are aligned and elastic tunneling 
is possible. Then, as \x 2 is pulled down the levels become 
misaligned and only inelastic tunneling can take place. 
Generally, this will cause the current to drop. When \ii 
is pulled down so much that a level corresponding to a 
transition involving an excited state (grey level of dot 2 in 
the diagrams) enters the bias window, elastic tunneling 
becomes possible again, leading to a rise in the current. 
When we move from this point into the triangle, along 
a line parallel to the base of the triangle, these levels 
remain aligned. Therefore, a line of elastic tunneling is 
observed parallel to the base of the triangle (depicted as 
a dark grey line). Going even farther up along the left 
leg of the triangle, levels arc again misaligned and only 
inelastic current flows. Beyond the top of the triangle, 
/i2(0, 1) falls below the drain electrochemical potential 
and the system is in Coulomb blockade. 

Moving down the upper leg from the top of the trian- 
gle, [12 is fixed (/^(0, 1) remains aligned with the drain 
electrochemical potential) and [i\ is pulled down. When 
a level corresponding to a transition involving an excited 
state in the left dot is pulled into the bias window (grey 
level in left dot in the diagrams), there are two paths 
available for electrons tunneling from the source onto the 
first dot. A different current can therefore be expected 
in the grey part in the upper right corner of the triangle. 

When the source-drain bias is inverted, the electrons 
move through the dot in the opposite direction and the 
roles of dot 1 and dot 2 are reversed. 

In principle, the different lines and regions of elastic 
and inelastic tunneling allow the full energy-level spec- 
trum to be determined of both dots. However, the visibil- 
ity of lines and regions depends strongly on factors such 
as the relative heights of the three tunnel barriers, the 
efficiency of inelastic tunneling processes (which again 
depends on the environment) and relaxation within the 
dots. For example, if relaxation in the first dot is much 
slower than the typical time for interdot tunneling, elas- 
tic current involving excited states in both dot 1 and 2 
can be observed. Another example: if the tunnel bar- 
rier between the source and dot 1 is much higher than 
the other two, the tunnel process from source to dot 1 
dominates the behaviour of the system and only excited 
states of dot 1 will be resolved in the current spectrum. 
In practice, the system should be tuned such that the pa- 
rameter of interest has the strongest effect on the current 
pattern, while the other factors can be neglected. 

From a high source-drain bias measurement as dis- 
cussed here, all four gate capacitances can be deduced. 
With those, all the energy scales such as charging ener- 
gies, electrostatic coupling energy, tunnel coupling and 
energy level spacing can be calibrated. 

In many of the experiments that are discussed in the 
following sections, the levels in the two dots are detuned 
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FIG. 31 Level diagrams for different detunings e between dot 
1 and dot 2 (dotted line in the bias triangle). Note that the 
average of the levels in the two dots is kept constant, and only 
the difference between the levels is changed. 



with respect to each other, while keeping the average of 
the two at a constant level. This is achieved by changing 
the gate voltages along a line exactly perpendicular to the 
base of the bias triangle. The resulting axis is commonly 
referred to as the detuning axis, which we denote by e. 
Figure [31] displays the level arrangements as a function 

Of £. 



B. Spin states in two-electron double dots 

The physics of one and two-electron spin states in sin- 
gle dots was discussed in section HVl In double quantum 
dots, electrons can be transferred from one quantum dot 
to the other by changing the electrostatic potentials using 
gate voltages. These interdot charge transitions conserve 
electron spin and are governed by spin selection rules, 
leading to a phenomenon called Pauli spin blockade. In 
order to understand this spin blockade, we first examine 
the spin states in the double dot system and the possi- 
ble transitions between these spin states, while neglecting 
processes that lead to mixing of these spin states. Such 
mixing terms will be introduced later, in section fVIILDI 

We focus on the two-electron regime, which has been 
the focus of many recent double dot experiments. We 
work in the region of the charge stability diagram where 
the occupancy of the double dot can be (0,1), (1,1), or 
(0,2). 

For (0,1) and (0,2) spin states, the spin physics is iden- 
tical to the single dot case since the left quantum dot is 
not occupied. We briefly repeat the description of the 
single-dot states, as discussed in Section HVl In the (0,1) 
charge state, the right dot contains a single electron. At 
zero magnetic field, the two spin states are degenerate. A 
finite magnetic field results in a Zeeman splitting between 
the spin- up and spin-down electrons, with E[=E^+Ez, 
where Ez is the Zeeman energy. In the (0,2) charge state, 
there are four possible spin states: the singlet, denoted 
by 5(0,2), and the three triplets T+(0,2), T (0,2) and 



T_(0,2). The spin parts of the wavefunctions of these 
states are: 

5(0,2) = (|T2| 2 )-|4 2 T2»/V2 (35) 
T+(0,2) = | Tats) (36) 
T (0,2) = (\hh) + \hh))/V2 (37) 
T_(0,2) = || 2 | 2 ), (38) 

where the subscript denotes the dot in which the electron 
resides. At zero magnetic field, the triplets T(0, 2) are 
separated by Est from the singlet ground state 5(0,2). 
An in-plane magnetic field Zeeman splits the triplet spin 
states. As in the single dot case, a perpendicular mag- 
netic field tunes Est and also Zeeman splits the triplet 
states. 

In the (1,1) charge state, the two-electron states are 
also spin singlets and triplets, but with the electrons in 
different dots: 



5(1,1) 

r (i,i) 
r_(i,i) 



(iTiia)-UiTa 
|TiTa> 

(ItilaHUita 
II1I2), 



(39) 
(40) 
(41) 
(42) 



The energy difference between the lowest-energy singlet 
and triplet states, J, depends on the tunnel coupling t c 
and the single-dot charging energy Ec- When the single- 
dot levels in the two dots are align e d, J = Atj/Eg in the 
Hubb ard approximation ( Burkardl . l200ll : iBurkard et all 
I1999D . Figure [32a depicts the energies of the two-electron 
spin states as a function of detuning between the two 
dots, for the case of negligibly small tunnel coupling. 
Since the three triplet states are degenerate, we denote 
them here by T(l, 1) and T(0, 2). The diagrams indicate 
the electrochemical potentials in left and right dot for 
three values of e. 

Due to the tunnel coupling the (1,1) and (0,2) charge 
states hybridize. In the case of spinless electrons, this 
would simply result in an avoided crossing between the 
(1,1) and (0,2) charge states that is characterized by a 
tunnel splitting, 2y/2t c . However, since the interdot tran- 
sitions preserve spin, the (1,1) singlet (triplet) states only 
couple to (0,2) singlet (triplet) states. As a result, the 
ground state singlets hybridize at a different value of 
detuning than the triplets, as depicted in Figure l3"2"b . 
At B=0, Est is typically in the range 0.4-1 meV in 
electrostatically defined dots in GaAs. This pushes the 
avoided crossings of the triplets far away from the avoided 
crossing of the singlets, which has two interesting con- 
sequences. First, the singlet-triplet energy difference J 
strongly depends on detuning, allowing simple electrical 
control over J. Second, the charge distribution of the sin- 
glet and triplet states are very different over a wide range 
of detunings. For example, at the value of detuning where 
the singlets have an avoided crossing, the electrons are in 
the charge state (|(1, 1) ) + |(0, 2)))/y/2 in case they form a 
spin singlet, but almost fully in the charge state |(1, 1) ) if 
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FIG. 32 (Color in online edition) Energies of the two-electron 
spin singlet and triplet levels in a double dot as a function of 
detuning e between the levels in the two dots for the case of 
(a) 5=0 and negligible tunnel coupling t c , (b) B=Q but sig- 
nificantly high value for t c , and (c) finite B and significantly 
high value for t c . The electrochemical potentials in the two 
dots are indicated for three values of detuning in the diagrams 
on top, for case (a). Here, (1,1) denotes the electrochemical 
potential of both degenerate states 5(1,1) and T(l,l). Note 
that other, higher-energy single-dot states will lead to avoided 
crossings at even larger values of detuning. 



they form a spin triplet. This spin-dependent charge dis- 
tributio n allows readout of the spin state throu gh charge 
sensing (|Engel et all 12004 iTavlor et all l2005t) . 

In a finite magnetic field, the triplet states are split by 
the Zeeman energy. Figure [32b shows the energy levels 
for a Zeeman splitting that exceeds the tunnel coupling. 
The application of a large magnetic field can be used to 
decouple the T + and T_ triplet states from the Tq triplet 



state, thus confining the relevant state space to 5 and 

The singlet-triplet energy difference J is often referred 
to as the exchange energy. In the strict sense of the word, 
exchange energy refers to the difference in Coulomb en- 
ergy between states whose orbital wavefunctions differ 
only in their symmetry (symmetri c for a spin singlet and 
antis ymmetric for a spin triplet) (|Ashcroft and Merminl . 
I1974D . In the case of two electrons in a double dot, J can 
also include a large contribution due to hybridization be- 
tween the (1,1) and the (0,2) and (2,0) charge states, es- 
pecially near one of the avoided crossings. In this sense, 
one could argue that J is not a true exchange energy. 
However, the double dot spin system can still be de- 
scribed by the Heisenberg spin Hamiltonian H = JS1S2, 
where Sx,2 are the electron spin operators. Therefore, J 
acts as an effective exchange coupling. To avoid confu- 
sion, we minimize reference to J as 'exchange energy' in 
this review. 



C. Pauli spin blockade 

The conservation of spin in electron tunneling leads 
to current rectification in dc transport in the two elec- 
tron regime. This effect, termed spin blockade or Pauli 
blockade, was first obse rved in experimen ts on vertically 
coupled quantum dots (lOno et all 120021 ). Later experi- 
ments in few-electron lateral dots combin ed charge sens- 
ing an d transport to study the effect (jjohnson et all 
l2005bl ). Measurements of transport in the Pauli blockade 
regime provided some of the first indications that the hy- 
pcrfine interaction plays an important role in the electron 
spin dynamics. Pauli blockade has also been utilized to 
implement spin-to-charge conversion for read out of the 
spin state of electrons in double quantum dots. 

The origin of Pauli blockade is schematically illus- 
trated in the insets of Fig. [33] At negative bias elec- 
trons are transferred through the device in the sequence 
(0,1)^(0,2)^(1, 1)->(0,1). In this cycle the right dot al- 
ways contains a single electron. Assume this electron is 
spin-up. Then, in the transition (0,1)^(0,2) the right dot 
can only accept a spin-down electron from the leads due 
to Pauli exclusion, and a 5(0, 2) state is formed. Sim- 
ilarly, only a spin-up electron can be added if the first 
electron is spin-down. From 5(0,2), one electron can 
tunnel to the left dot and then out to the left lead. 

In contrast, when the bias voltage is posi- 
tive charge transport proceeds in the sequence 
(0,1)^(1, l)->(0,2)->(0,l) and the left dot can be 
filled from the Fermi sea with either a spin-up or a 
spin-down electron, regardless of the spin of the electron 
in the right dot. If the two electrons form a singlet state 
5(1,1), the electron in the left dot can transfer to the 
right dot forming 5(0, 2). However, if the electrons form 
one of the triplet states T(l, 1), the electron in the left 
dot will not be able to tunnel to the right dot because 
T(0, 2) is too high in energy. The system will remain 
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FIG. 33 Current (7) measured as a function of source-drain 
voltage (V) in a vertical double dot system. Non-zero current 
is measured over the entire range of negative voltage. For pos- 
itive bias, current is blocked in the range 2<V<7 mV. At bias 
voltages exceeding 7 mV, the (0,2) triplet state becomes ac- 
cessible and Pauli blockade is lifted. Insets: Device schematic 
and energy level configuration a t positive and ne gative bias 
voltages. Data reproduced from lOno et all l|2002T ). 



stuck in a (1,1) charge state until the electron spin 
relaxes. Since the T\ time can approach milliseconds, 
the current in this direction is negligible and the dot 
is said to be in spin blockade. Because it is the Pauli 
exclusion principle that forbids the electrons to make a 
transition from a T(l,l) state to S(0, 2), this blockade 
is also referred to as Pauli blockade. 

The spin blockade effect leads to current rectification 
in dc transport. Figure [33J shows an I-V curve taken 
from a vertical double dot. Non-zero current is ob- 
served for negative voltages. For positive bias voltage, 
spin blockade is observed in the range 2-7 mV. Once the 
bias voltage exceeds the singlet-triplet splitting Est of 
the (0,2) charge state, also the T(0, 2) state is energeti- 
cally accessible from T(l, 1) and the blockade is lifted. A 
theor etical model reproduces the o bserved current pat- 
tern ( Fransson and Rasanderl l2006f) . 



Note that the s pin blockade can be lifted by photon- 

the photon then 



pm 

assisted tunneling ( San chez et all 



supplies the energy needed to make the transition from 
T(l, 1) to T(0, 2). Interestingly, it is predicted that for 
a suitable choice of the applied photon frequency, the 
resulti ng pumped current can have a large s pin polariza- 
tion (|Cota et all , 120051: ISanchez et al\ , 120061 ) . 

Pauli blockade has also been observed in lateral double 
dot systems. In these systems the tunnel rates and offset 
energies are easily tuned. Moreover, devices equipped 
with a charge sensor can be used to measure the average 
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FIG. 34 (Color in online edition) Double dot current mea- 
sured as a function of Vl and Vr in the one and two-electron 
regime. In the one-electron regime (a, c) the finite bias tri- 
angles at negative bias mimic the positive bias data, except 
for an overall change in the sign in the current. However, 
in the two-electron regime, charge transport shows a strik- 
ing asymmetry when the sign of the bias voltage is changed. 
At negative bias in the two-electron regime charge transport 
is blocked, except near the edges of the finite bias triangles, 
where exchange of electrons with the leads lifts the spin block- 
ade. Insets show simple rate equ ation predictions of ch arge 
transport. Data reproduced from iJohnson et al\ l|2005bl ). 



occupancy of the double dot during charge transport (see 
Fig. [5]for a device image). Figure 1341 shows experimental 
data from measuring current as a function of Vl and 
Vr in t he one- and two-electr on regimes for both signs 
of bias ( Johnson et all , l2005bl ) . Apart from a change in 
the sign of current when the voltage is changed, the one- 
electron data are mirror images of each other for positive 
and negative bias. This is in contrast with data acquired 
in the two-electron regime, where current flows freely for 
positive bias but is strongly suppressed at negative bias 
due to Pauli exclusion (the voltage bias conv ention in this 
paper is opposite to Ref. (|Ono et aE|2002), so blockade 
is observed at negative bias). At negative bias, current is 
only observed along the edges of the bias triangles, where 
an electron can be exchanged with the leads lifting the 
spin blockade (see diagrams in Fig. 
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FIG. 35 (Color in online edition) Charge sensor conductance, 
g s measured as a function of Vl and Vr in the one and two- 
electron regimes. The charge sensor conductance in the lower- 
left finite bias triangle in the one electron regime is a weighted 
average of the (0,0), (0,1), and (1,0) charge sensing signals. 
In the two-electron Pauli blockade regime, the charge sensor 
conductance in the finite bias triangles is pinned to the (1,1) 
charge sensing value. This indicates that charge transport 
is blocked by the (1,1)— >(0,2) charge transition. Insets show 
simple rate equation pr edictions of charge senso r conductance. 
Data reproduced from IJohnson" et ad (|2005bT ). 



Charge sensing measurements of the time-averaged oc- 
cupancy of the double quantum dot during transport di- 
rectly demonstrate that the current rectification is due 
to a blocked inter-dot charge transition. Figure l3"5l shows 
the charge sensor conductance measured as a function of 
Vl and Vr in the one- and two-electron regimes for both 
positive and negative bias. For the two-electron case at 
positive bias charge transport in the lower-left bias trian- 
gle proceeds in the sequence (0,l)-^(0, 2)— >(1, 1)— >(0,1). 
The charge sensing signal in the finite bias triangles is a 
weighted average of the (0,1), (0,2) and (1,1) charge sens- 
ing levels. At negative bias in the two-electron regime 
charge transport in the lower-left bias triangle follows 
the sequence (0,l)->(l,l)->(0,2)->(0,l). The data in Fig. 
tJSTd) show that the charge sensing conductance in the 
finite bias triangles is practically identical to the back- 
ground (1,1) charge sensing signal. These data indicate 



that the charge transition from (1,1) to (0,2) is the lim- 
iting step in the current: precisely what is expected for 
a double dot in spin blockade. 



D. Hyperfine interaction in a double dot: Singlet- Triplet 
mixing 

Early experiments in semiconducting heterostructures 
in the quantum Hall regime demonstrated that spin po- 
larized cur rents could b e used to polarize the nuc lei in the 
substrate (jDixon et all Il997t IWald et all Il994l) . These 
measurements gave a clear indication that electronic ef- 
fects can have a strong influence on the nuclear spin sys- 
tem. So far we have ignored the consequences of the 
hyperfine interaction in few-electron quantum dots, but 
early indications of this important interaction were vis- 
ible in the firs t Paul i blockade experiments by Ono et 
al. ( Ono et al\ , I2002T ) . In this section we review several 
recent experiments that have shown that hyperfine ef- 
fect can have profound consequences on the electron spin 
dynamics in GaAs quantum dots. 

In GaAs quantum dots each electron spin is coupled to 
a bath of nuclear spins through the contact hyperfine in- 
teraction (see Section [VII. Bp . The importance of the hy- 
perfine field becomes apparent when considering two spa- 
tially separated electron spins in a double dot structure. 
Each electron has a distinct orbital wavefunction and av- 
erages over a different set of nuclei. As a result, each 
electron experiences a slightly different nuclear field. The 
difference in the nuclear fields, A£?at jZ , couples the singlet 
and triplet spin states. For example, the z-component of 
the nuclear field couples 5(1,1) and To (1,1), with the 
Hamiltonian (in the basis (5(1, 1), Tq(1, 1)): 
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gUB&B 



N,z 







(43) 



Since 5 and To are not eigenstates of this Hamiltonian, 
the off-diagonal terms will drive rotations between 5 and 
T). Similarly, the x-component and the y-component of 
the nuclear field mix T+(l, 1) and T_(l, 1) with 5(1, 1). 

Figure. [3"rjf a) shows measurements of the "leakage cur- 
rent" in the Pauli spin blockade regime in vertical dou- 
ble dots as a function of magnetic field for tw o differ- 
ent sweep directions ( Ono and Taruchal l2004h . Upon 
increasing the magnetic field, the leakage current was 
nearly constant until B—0.5 T, where a sudden increase 
in the leakage current was measured. The leakage current 
decreased suddenly for fields exceeding 0.9 T. Measure- 
ments of the leakage current for the opposite magnetic 
field sweep direction showed hysteretic behavior. The 
amount of hysteresis decreased for slower magnetic field 
sweep rates. In the high leakage current regime (B w 0.7 
T), the leakage current showed surprising oscillations in 
time. The frequency of these oscillations was a sensitive 
function of the external field (see Fig. (b)). By mov- 
ing in and out of the Pauli spin blockade regime using 
gate voltages, Ono et al. determined that the oscillatory 



time dependence of the leakage current developed on ; 
5 minute timescale. Moreover, the leakage current wa 
modified by the application of cw radiation at the 71 G 
or 69 Ga NMR lines. All these aspects indicate that th 
nuclear spins play a major role in the observed behaviour 
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FIG. 36 (a) Pauli blockade leakage current as a function of 
magnetic field for increasing and decreasing magnetic field 
sweeps, (b) Leakage current as a function of time for fields 
in the range of 0.7 T (bottom trace) to 0.85 T (top trace), 
(c) Transient behavior of the leakage current measured by 
moving in and out of t he Pauli blockade regim e using Vsd- 
Data reproduced from lOno and Taruchal (|2004T ). 



The leakage current in the Pauli spin blockade region 
occurs due to spin relaxation from T_(l,l) to 5(1,1) 
and the hysteretic behavior observed in Fig. [36b can 
be explained in terms of triplet-to-singlet relaxation via 
hyperfine-induced flip-flops with the spins of the lattice 
nuclei in the dot. In the measured device, the detuning 
between the two dots corresponds to a point just to the 
right of the avoided crossing between 5(1, 1) and 5(0, 2)in 
Fig. [32b . Here, 5(1, 1) is slightly higher in energy than 
T(l, 1). This energy separation is about 10 /ieV in the 
measured device. At zero magnetic field, this energy mis- 
match makes the "flip-flop" mechanism between electron 
and nuclear spins inefficient. However, the energy differ- 
ence is compensated by the Zeeman energy at a magnetic 
field of about 0.5 T, which is comparable to the mag- 
netic field where a current step is observed (indicated by 
a triangle in Fig. [36b). On approaching this particular 
magnetic field, T(_l,l) and 5(1,1) become degenerate 
(see Fig. 13Tb). Then, the hyperfine-induced T_(l,l)-to- 
5(1, 1) relaxation becomes efficient, because energy as 
well as spin is conserved in "flip-flops" between the elec- 
tronic and nuclear spin systems. Many such flip-flops 
lead to a finite nuclear spin polarization, which acts back 
on the electron as an effective magnetic field (see Sec- 
tion [VILA] Because the nuclear spin has a long lifetime 
(on the order of minutes), a nuclear spin polarization 
accumulates to sustain the T_-5 degeneracy condition 
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FIG. 37 (a) Effect of the Zeeman energy on 5(1, 1) and 
T(l, 1) states, which are separated by an energy J. The 
5(1, 1) and TL(1, 1) become degenerate when the Zeeman en- 
ergy is equal to J. (b) Magnetic field dependence of the leak- 
age current measured at various values of source-drain volt- 
ages, Vsd, in the Pauli blockade region. Each curve is offset 
by 1 p A to the top. Data reproduced from iTarucha et al\ 
(|2006l h 



on sw eeping down the external field ( Ono and Taruchal 
I2004D . An increasing nuclear field thus compensates the 
decreasing external magnetic field; in other words, the 
effective magnetic field resulting from nuclear spin po- 
larization adds to the external field. From the consid- 
erations in Appendix we see that this implies that 
the electron spin is changed by A5 Z =+1; this is con- 
sistent with hyperfine-induced transitions from T_ to 5. 
The hyperfine interactions are thus the origin of the hys- 
teretic loop. We note that the similar effect was well 
studie d in ESR experiments on two-dimensiona l electron 
gases (|Dobers et all Il988t iTeraoka et all , I2004D . 

More detailed experiments on the hysteretic behav- 
ior are perfo rmed for a ve r tical double dot, as shown 
in Fig. [37b (|Tarucha et all 120061 ). The observed hys- 
teretic behavior significantly depends on the source-drain 
voltage, that is, the hysteretic loop becomes small and 
shifts to the lower field for the higher source-drain volt- 
age Vsd ■ This is well understood in terms of the decrease 
of singlet-triplet energy splitting, which is estimated from 
the threshold field (arrows) as a measure: increasing Vsd 
increases the detuning between two dots. As can be seen 
from Fig. [32b . this decreases the energy difference be- 
tween T(l, 1) and 5(1,1) and therefore a smaller mag- 
netic field is needed to compensate for it. 

Further insight into the role of the hyperfine interac- 
tion on the electron spin dynamic s was gained in exper i- 
ments on lateral quantum dots ( Koppens et all l2005h . 
These experiments measured the Pauli spin blockade 
leakage current as a function of the external magnetic 
field and of the exchange splitting separating the (1,1) 
singlet and triplet spin states. Figure[38]explores the tun- 
nel coupling and magnetic field dependence of the Pauli 
blockade in plots of the double dot current as a func- 
tion of Vl and Vr. For strong interdot tunnel couplings 
current rectification due to Pauli blockade is observed 
(Fig. I3"5h). When the tunnel coupling is reduced, the 
Pauli blockade is lifted and a substantial current starts to 
flow, as shown in Fig. 138b . Increasing the magnetic field 
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energy splitting between these states is less than or com- 
parable to the nuclear field scale, £J„ UC ~100 neV. This 
condition is achieved over the entire range of detunings 
for small tunnel coupling but only at large detuning for 
strong tunnel coupling. An external field splits off the 
ms = ±1 triplet states, T+ and T_ by the Zeeman energy 
(Fig. I3"2"b). When B < Bm these states also rapidly mix 
with 5(1, 1) due to the inhomogeneous hyperfine fields. 
However, when B > Bn the T + and T_ states do not mix 
with 5(1, 1) anymore, and the spin blockade is recovered. 
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FIG. 38 (Color in online edition) Transport in the Pauli 
blockade regime as a function of Vl and Vr. (a) In the limit 
of strong tunnel coupling current suppression due to Pauli 
blockade is observed. For weak tunnel coupling, the Pauli 
blockade leakage current displays a striking magnetic field 
dependence. At B—0 mT (b), Pauli blockade is lifted near 
the (l,l)-(0,2) charge transition (near zero detuning). In con- 
trast, for _B=100 mT (c), current is suppressed due to Pauli 
blockade. Data reproduced from iKoppens et al\ (|2005T ) . 



to 100 mT quenches this leakage current (see Fig. l38b). 
In all large current is observed when the volt- 

age bias exceeds the (0,2) singlet-triplet energy difference 
Est- 

These data can be explained by considering the de- 
pendence of the two-electron spin states on magnetic 
field and exchange splitting, as illustrated in Fig. [32l 
(see a lso (Cois h and LossL 2005; Ijouravlev and Nazarovl 
l2006h ) . For small tunnel coupling (Fig. 132a). the singlet 
5(1, 1) and the three triplets T(l, 1) are nearly degener- 
ate over the entire range of detuning. Increasing the tun- 
nel coupling results in a finite exchange splitting between 
5(1, 1) and all T(l, 1) states (Fig. [35b). The inhomoge- 
neous hyperfine fields mix 5(1,1) and T(l,l) when the 
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FIG. 39 (Color in online edition) Charge sensor conductance, 
g 3 , measured as a function of Vl and Vr using the Ti pulse 
sequence. The triangular shaped region in the (0,2) region of 
the charge stability diagram, termed the "pulse triangle", is 
due to spin blocked interdot charge transitions. A relaxation 
time is determined by measuring the decay of this signal as a 
function of time (see (a),(b)). Ti shows a strong dependence 
on magnetic field. This is apparent in the B=0 mT data of 
(c) , where near zero detuning the spin states have completely 
relaxed on 8 /is timescales. For long times, tm=80 /j,s and 
B=0 mT the spin states have completely re laxed and the 

ulse t riangle is absent. Data reproduced from iJohnson et al\ 

2005d ). 



Time-resolved techniques have been used to measure 
the hyperfme-induced relaxation of a sp in triplet state 
in a two-electron double q uantum dot ( Johnson et all 
l2005d: iPetta et all l2005a| ). These experiments used 
pulsed gate techniques to prepare a spin triplet state and 
then measure the decay of that spin state using spin- 
to-charge conversion. The pulse experiment is performed 
near the (l,l)-(0,2) region of the charge stability diagram. 
Gates are set in (0,1) to empty the left dot. A pulse then 
shifts the gate voltages to the (1,1) region of the charge 
stability diagram. A spin-up or spin-down electron enters 
the left dot forming a spin singlet or spin triplet state. A 
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The full dependence of the spin triplet relaxation time 
Ti as a function of magnetic field and detuning is plotted 
in Fig. [40] At small detunings near the interdot charge 
transition, T\ displays a strong dependence on magnetic 
field. Simply increasing the field from to 100 mT ex- 
tends Ti from microsecond to millisecond timescales. At 
larger values of the detuning, T\ is nearly independent of 
magnetic field. This indicates that hyperfine-mediated 
spin relaxation is no longer dominant, but that relax- 
ation is instead due to a coupling to the leads (which is 
independent of magnetic field). Experimental data are fit 
using a simple model of spin relaxation from T(l, 1) to 
5(1, 1) followed by inelastic decay from 5(1, 1) to 5(0, 2). 
The model assumes hyperfine driven spin relaxation as 
well as a spin relaxation contribution from coupling to 
the leads at large detunings. Best fits to the model give 
=2. 8 mT, which is consistent with the estimated num- 
ber of nuclei in the dot. 



FIG. 40 (Color in online edition) Spin relaxation time Ti of 
the double-dot spin triplet plotted as a function of detuning 
for a range of external magnetic fields. At low detunings, a 
strong magnetic field dependence is observed due to hyperfine 
driven spin relaxation. At large detunings, spin relaxation 
occurs due to coupling to the leads, and is independent of 
magnetic field. Data are fit using a simple model of hyper- 
fine driven relaxation and ther mally activated coupli ng to the 
leads. Data reproduced from Ijohnson et al\ (|2005a ). 



spin triplet state is formed 75% of the time. To measure 
the relaxation time T\ of the spin triplet state, a third 
pulse is applied to the device which tilts the double well 
potential so that 5(0, 2) is the ground state. In order 
for the left electron to tunnel to the right dot, the (1,1) 
triplet state must spin relax to 5(1, 1) and then tunnel to 
5(0,2). By measuring the occupancy of the double dot 
as a function of the time spent in the biased configuration 
the spin relaxation time can be determined. 

Representative data are shown in Fig. [39] as a function 
of magnetic field and time in the biased configuration, 
tm . In Fig. [35] (a) g s is plotted as a function of Vl and 
Vr with i?=100 mT and tm= 8 [is. A triangular shaped 
signal ("pulse triangle") appears in the (0,2) region of 
the charge stability diagram, which is indicative of spin 
blocked transitions. For £?=100 mT and r m — 80 fis the 
signal in the pulse triangle reduces to a value approach- 
ing the background 5(0, 2) charge sensing level, indicat- 
ing that tm^tst and the spin blocked triplet states have 
relaxed to (1,1)S and tunneled to 5(0, 2). In addition to 
the observed time dependence a strong magnetic field 
effect is observed. Reducing B from 100 mT to mT 
quenches the triplet state signal near the interdot charge 
transition for the tm=8 [is data, which implies that spin 
relaxation is much faster near zero field at small detun- 
ings. Finally, with B=0 mT and tm =80 [is the signal in 
the pulse triangle is completely absent indicating com- 
plete spin relaxation. 



IX. COHERENT SPIN MANIPULATION 

A. Single-spin manipulation: ESR 

A variety of techniques can be used to coherently 
drive transitions between the Zeeman split levels of a 
single electron. The most well-known approach is elec- 
tron spin resonance (ESR) , whereby a rotating magnetic 
field, B\, is applied perpendicularly to the static field B 
along z, and on-resonance with the spin-flip tra nsition 
energ y (f ac — gfisB/h), as illustrated in Fig. l4Tl|PooleL 
1983). Alternatively, spin rotations could be real- 
ized by electrical or optical excitation. Electric fields 
can couple spin states through the spin-orbit interac- 



tion dDebald and Emarvl 120051: iDobrowolska et al. . 1982 : 
Golovach et all 20061 iKato et all l2003bl : ISchulte etal. 
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20051 ) (see also Section IVII.Ap , by making use of an in 



homogeneous static magnetic field dTokura et all I2006D . 
or by g-tensor modulation ( Kato et all . l2003ah . Op- 



tical ex citation can indu c e spin flips via Raman tran- 
sitions (llmamoglu et al. , 1999) or the optical Stark 
effect ([Gupta et all . l200lh . To date, driven coher- 
ent rotations of a single spin in a solid have only 
been realized using ESR, and only in a few spe- 



cific systems ( Hanson et all 20061: Jelezko et 
iRuear et all 120041: IXiao et all l2004h . including in a 
quantum dot ([Koppens et all 12006). In addition, the 
free precession of an electron spin in a quantum dot has 
been observed with op tical techniques ( Putt et aZI . 120051 : 
iGreilich ef"aZl . l2006aj ). 

The quantu m dot ESR experime nt was realized by 
Koppens et al. ( Koppens et all 12006). and is inspired by 
the idea of Engel and Loss to tune a single quantum dot 
to Coulomb blockade with the electrochemical potential 
alignment as shown in Fig. [42k . such that the Coulomb 
blockad e is lifted when the electro n spin is repeatedly 
flipped ( Engel and LossL [200 ll . |200^ ) . In practice, this re- 
quires excitation in the microwave regime, as the Zeeman 
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a) 




FIG. 41 Motion of the electron spin during a spin resonance 
experiment, (a) The motion as seen in a reference frame that 
rotates about the z axis at the same frequency f ac as the spin 
itself and the resonant rotating magnetic field, B\. Naturally, 
the rotating field B\ lies along a fixed axis in this rotating 
reference frame. An observer in the rotating frame will see 
the spin simply precess about B\, with a rate u)i. (b) An 
observer in the lab reference frame sees the spin spiral down 
over the surface of the Bloch sphere. 



splitting must be well above the thermal energy. Further- 
more, the alternating electric fields that are unavoidably 
also generated along with the alternating magnetic held, 
can kick the electron out of the dot via photon-assiste d 
tunneling (PAT) processes ( Platero and Aguadol . l2004h . 
In early attempts to detect ESR, PAT processes and heat- 
ing of the electron reservoirs lifted the blockade long be- 
fore enough power was applied to lift the blockade by 
ESR ((Hansorj, [2005). Efforts to suppress the electric field 
component while max imizing the magneti c component, 
via optimized cavi ties dSimovic et al\. 120061) or microfab- 
ricated striplines (jKoppens et all 120061 ) , have so far not 
been sufficient to overcome this problem. 



(a) 



(b) 





FIG. 42 (Color in online edition) Schematic diagrams of a 
single quantum dot and a double quantum dot, illustrating 
electrical detection of ESR. In both cases, transport through 
the system is blocked, but the blockade is lifted when the ESR 
condition is satisfied and the spin of the electron is flipped, (a) 
The two electrochemical potential levels shown are the spin- 
up and spin-down levels of the lowest single-electron orbital. 
The system is in Coulomb blockade, (b) The levels shown 
are as in the discussion of spin blockade in double dots in 
Section rVIII. CI The Zeeman sublevels are not shown. 

Instead, ESR detection in quantum dots has been real- 
ized using two quantum dots in series, tuned to the spin 
blockade regime described in section IVIII.CI The two 
dots are weakly coupled, and subject to a static magnetic 
field B, such that the To state is mixed with the singlet 
but the T± states are not. Current is then blocked as 



soon as the double dot is occupied by two electrons with 
parallel spins (one electron in each dot), but the block- 
ade is lifted when the spin in the left or the right dot is 
flipped (Fig. HUb). 

In this double dot ESR detection scheme, the relevant 
transition occurs between the two dots. This transition is 
not affected by temperature broadening of the leads. As 
a result, ESR detection can be done with Zeeman split- 
tings much below the thermal energy, and thus with ex- 
perimentally much more accessible frequencies. Further- 
more, by applying a large voltage bias across the double 
dot structure, photon-assisted tunneling processes can be 
greatly suppressed. 

The ESR response is seen clearly in transport measure- 
ments through the double dot. When the static magnetic 
field is swept, clear peaks in the current develop at the 
resonant field when an AC magnetic field is turned on, 
as seen in Fig. 23] (the alternating magnetic field B ac 
can be decomposed into a component with amplitude 
B\ = B ac /2 rotating in the same direction as the spin 
precession and responsible for ESR, and a component 
rotating the opposite way, which hardly affects the spin 
because it is very far off-resonance). The linear depen- 
dence of the satellite peak location on the RF frequency, 
which is the characteristic signature of ESR, is clearly 
seen. The characteristic signature of ESR is the linear 
dependence of the satellite peak location on the RF fre- 
quency which is clearly seen in the data when the RF 
frequency is varied from 10 to 750 MHz. A linear fit 
through the top of the peaks gives a (/-factor with mod- 
ulus 0.35 ± 0.01, which is similar to the values obtained 
from high-bias transport measurements in single dots (see 
Section |TV}. 

In order to observe coherent single-spin rotations, the 
system is pulsed into Coulomb blockade while B ac is ap- 
plied. This eliminates decoherence induced by tunnel 
events from the left to the right dot during the spin ro- 
tations. The experiment then consists of three stages 
(Fig. I44p : initialization through spin blockade in a sta- 
tistical mixture of tt an d II, manipulation by a RF burst 
in Coulomb blockade, and detection by pulsing back for 
projection (onto 5(0,2)) and tunneling to the lead. If 
one of the electrons is rotated over (2n + l)ir (with inte- 
ger n), the two-electron state has evolved to || (or ||), 
giving a maximum contribution to the current (as before, 
when the two spins are anti-parallel, one electron charge 
moves through the dots). However, no electron flow is 
expected after rotations of 2n7r, where two parallel spins 
are in the two dots after the RF burst. 

The measured dot current oscillates periodically with 
the RF burst length (Fig. [45)) , demonstrating driven, co- 
herent electron spin rotations, or Rabi oscillations. A 
key signature of the Rabi process is a linear dependence 
of the Rabi frequency on the RF burst amplitude, B\ 
(friabi = g^BBi/h). This is verified by extracting the 
Rabi frequency from a fit of the current oscillations of 
Fig. 145b with a sinusoid, which gives the expected linear 
behavior (Fig. 145b . inset). The maximum B\ that could 
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FIG. 43 (Color in online edition) (a) Energy levels of the two- 
electron spin states in the double dot. ESR can drive transi- 
tions between the states with parallel spins to the states with 
anti-parallel spins, thereby lifting spin blockade, (b) Mea- 
sured current through two quantum dots in the spin blockade 
regime in the absence (blue) or presence (pink) of an AC mag- 
netic field. The pink curve is offset by 100 fA for clarity. At 
zero-field, all three triplets are admixed with the singlet, so 
here the current is never blocked. With the AC field turned 
on, two satellite peaks develop at the electron spin resonance 
condition. Inset: the amplitude of the ESR peaks increases 
linearly wi th RF power (oc B^ r ) before saturation occurs, as 
predicted (En gel and Loss! . 1200 if ), (c) Measured current (in 
color-scale) through the two dots as a function of static mag- 
netic field and excita tion frequency. Data reproduced from 
iKoppens efaH (|2006T ). 
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FIG. 44 (Color in online edition) The control cycle for coher- 
ent manipulation of the electron spin via electron spin reso- 
nance. 



be reached in the experiment was ~ 2 mT, correspond- 
ing to 7r/2 rotations of only 25 ns (i.e. a Rabi period of 
~ 100 ns). The main limitation that prevented the use 
of larger Bi's was still photon-assisted tunneling, even 
in this double dot detection scheme. From the spread in 
the nuclear field and the RF field strengths that could 
be applied, a fi delity of 75% was est imated for intended 
180° rotations (jKoppens et all 120061 ). 

The oscillations in Fig. 145b remain visible throughout 
the entire measurement range, up to 1/zs. This is strik- 
ing, because the Rabi period of > 100 ns is much longer 
than the time-averaged coherence time T 2 * of roughly 
25 ns, caused by the nuclear field fluctuations (see sec- 
tion IVII.Bp . The slow damping of the oscillations is 
only possible because the nuclear field fluctuates very 
slowly compared to the timescale of spin rotations and 
because other mechanisms, such as the spin-orbit inter- 
action, disturb the electron spin coherence only on even 
longer timescales. 

Finally, we note that in this first ESR experiment, the 
excitation was on-resonance with either the spin in the 
left dot or the spin in the right dot, or with both, de- 
pending on the value of the random nuclear fields in 
each of the two dots. In all cases, blockade is lifted and 
ESR is detected. In future experiments, controllable ad- 
dressing of the spins in the two dots separately can be 
achieved through a gradient in either the static or the 
oscillating magnetic field. Such gradient fields can be 
created relatively easily using a ferromagnet or an asym- 
metric stripline. Alternatively, the resonance frequency 
of the spins c an be selectively shifted using local g-factor 
engin eering ( Jiang and YablonovitchL l200lb ISalis et all , 
120011) . 



B. Manipulation of coupled electron spins 

It has been shown that single spin rotations com- 
bined with two-qubit operations can be used to cre- 
ate basic quantum gates. For example, Loss and Di- 
Vincenzo have shown that a XOR gate i s implem ented 
by co mbining single-spin rotations with \J SWAP opera- 
tions ( Loss and DiVincenzol . ll998f) . In the previous sec- 
tion experiments demonstrating single-spin manipulation 
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(1,1) - (0,2) charge transition (see Fig. [32]) . The abso- 
lute number of electrons in the double dot is determined 
through charge sensing with the QPC. 
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FIG. 45 (Color in online edition) Coherent single-spin rota- 
tions, (a) The dot current - reflecting the spin state at the 
end of the RF burst - oscillates as a function of RF burst 
length (curves offset by 100 fA for clarity) . The period of the 
oscillation increases and is more strongly damped for decreas- 
ing RF power (P represents the estimated power applied to 
the on-chip stripline). Each measurement point is averaged 
over 15 seconds. The solid lines are obtained from numeri- 
cal computation of the time evolution of the electron spins, 
using a simple Hamiltonian that includes B, B\ and a Gaus- 
sian distribution of nuclear fields in each of the two dots, (b) 
The oscillating dot current (in colorscale) is displayed over 
a wide range of RF powers (the sweep axis) and burst du- 
rations. The dependence of the extracted Rabi frequency 
fRabi on RF power i s show n in the inset. Data reproduced 
from lKoppens et al\ (|2006T ). 



were reviewed. To implement more complicated gate se- 
quences, two-qubit interactions are required. In this sec- 
tion we review experiments by Petta et al. that have 
used fast control of the singlet-triplet e nergy sp litting in a 
double dot system to demonstrate a V SWAP operation 
and implement a singlet-triplet spin echo pu lse sequence , 
leading to microsecond dephasing times ( Petta et all 
12005th . 

A few-electron double quantum dot is used to isolate 
two electron spins (the device is similar to that shown 
in Fig. [5]). The device is operated in the vicinity of the 
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FIG. 46 (Color in online edition) (a) Schematic representa- 
tion of pulse sequence used to measure the singlet state de- 
cay, (b) Energy of the two-electron spin states as a function 
of detuning for the singlet-triplet qubit. Zero detuning is de- 
fined here as the value for which the energies of 5(1,1) and 
5(0,2) are equal. At positive detunings, the ground state is 
5(0,2). For negative detunings, and at finite fields, 5 and 
To are nearly degenerate. At zero magnetic field, the triplet 
states are degenerate, (c) Singlet state probability measured 
as a function of detuning and magnetic field for rs=200 ns 
2>T 2 *. (d) Hybridization of the (1,1) and (0,2) charge states 
results in a gate- voltage-tunable ener gy splitting, J(e). Data 
reproduced from iPetta et al\ (|2005bl ). 

The energy of the two-electron spin states as a function 
of detuning is illustrated in Fig. 1467b) (see also Fig. [32] 
for a zoom-out). At positive detuning the ground state 
is 5(0,2). The triplets, T +A -(0,2) are off-scale in this 
plot (Est = 0.4 meV). For sufficiently negative detun- 
ings, 5(1,1) and Tb(l,l) are nearly degenerate. An exter- 
nal magnetic field splits off T+(l,l) and T_(l,l) by the 
Zeeman energy. Near e=0, the singlet states 5(1,1) and 
5(0,2) are hybridized due to the interdot tunnel coupling 
t c . This hybridization results in an energy splitting J(e) 
between Tb(l,l) and 5(1,1) that is a sensitive function of 
the detuning. 

This energy level diagram can be mapped out experi- 
mentally by measuring the decay of a initially prepared 
singlet state as a function of magnetic field and detuning. 
The pulse sequence is schematically shown in Fig. l46la) 
(see also Fig. [32j). The singlet state, 5(0,2), is prepared 
at positive detuning. A pulse is applied to the device 
which lowers the detuning, so that the two electrons 
forming the spin singlet state are separated (one elec- 
tron in each dot, 5(1, 1)). The spins are then held in the 
separated configuration for a time r s >>T 2 *. At locations 
in the energy level diagram where 5 is nearly degenerate 
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with one of the triplet states fast spin mixing will occur, 
thereby reducing the singlet occupation P s . Fig. HCTc) 
shows Ps as a function of B ext and e. A strong magnetic 
field dependent signal is observed, corresponding to the 
5(1,1)-T + (1,1) degeneracy. For dctunings more negative 
than —1.5 mV, 5(1,1) and To(l,l) are nearly degener- 
ate resulting in a reduced singlet state probability. J(e) 
is extracted from the 5(1,1)-T + (1,1) degeneracy and is 
plotted in Fig. l46Td). As can be seen from this figure, a 
shift in detuning of just a few mVs reduces J from a few 
fieV to well below 100 neV. 

Hyperfine fields were shown in Section IVIII.DI to lead 
to current leakage in the Pauli blockade regime and to en- 
hanced low-field spin relaxation rates. One relevant ques- 
tion for quantum information processing is how long two 
spatially separated electron spins retain coherence in this 
solid state environment. To directly measure this time a 
two-electron spin singlet state is prepared, then the elec- 
tron spins are spatially separated, and finally correlations 
between the electron spins are measured at a later time. 
This experiment is performed using fast electrical con- 
trol of J. In the spatially separated (1,1) configuration 
the electron spins experience distinct hyperfine fields. In 
a semiclassical picture, the electron spins precess about 
the local hyperfine fields. Spatial variations in Bn, ABn 
result in different spin precession rates for the spatially 
separated electron spins. This drives a rotation between 
5(1, 1) and the triplet states. To measure the rotation 
rate in the hyperfine fields the separation time t s is var- 
ied. 

The rotation rate in the presence of the hyperfine fields 
is determined by performing spin-to-charge conversion af- 
ter a separation time ts- Detuning is increased and the 
double well potential is tilted so that 5(0,2) is the ground 
state. A separated singlet state 5(1,1) will adiabatically 
follow to 5(0,2), while the triplets T+,0 _(1,1) will remain 
in a spin blocked (1,1) charge state for a long time, T\. 
A charge sensing signal of (0,2) indicates that the sep- 
arated spins remain in the singlet state, while a charge 
signal of (1,1) indicates that the separated spins rotated 
into a triplet state. 

Figure [47| shows the singlet state probability as a func- 
tion of separation time ts, P(ts) for B=0 and B=100 
mT. For rg^T^ we find Ps~l. Ps exhibits a Gaussian 
decay on a 10 ns timescale and has long time saturation 
values of 0.5 (0.7) for B=0 (£=100 mT). The data are fit 
using a simple semiclassical model of the hyperfine fields 
assuming an average over many nuclear spin configura- 
tions. Best fits to the data give B N =2.3 mT and T 2 * = 10 
ns. The theoretical curves account for a measurement 
contrast of ~60%. Long time Ps values reflect the spin 
state degeneracy at zero and finite fields. The measure- 
ment shows that the separated spins lose coherence in ~ 
10 ns. 

Two electron spins can be manipulated by fast control 
of the singlet-triplet energy splitting J. The Hamiltonian 
of the two-electron system in the basis (|5), |Xb)) can be 
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FIG. 47 (Color in online edition) Singlet state probability, Ps, 
measured as a function of separation time, ts. Data points 
are acquired at B=0 and B=100 mT. Solid lines are best 
fits to the data using a semiclassic al model of t he hyp erfine 
interaction. Data reproduced from iPetta et al\ (|2005bl ). 



approximated for zero and negative detuning by 

H _ ( -J( £ ) g^BAB NtZ 
\ gfi B AB Nz 



(44) 



Note that |5) and |T ) are defined as the lowest-energy 
spin singlet state and spin Tn-triplet state, respectively. 
Whereas the |Tn) state is almost a pure (1,1) orbital state 
in the region of interest, the state |5) has an orbital char- 
acter that changes with detuning due to the hybridization 
of 5(1,1) and 5(0,2) (see Fig. 135Tb)). Since the differ- 
ence between the nuclear fields in the dots, ABn, acts 
on 5(1,1) but not on 5(0,2), the Hamiltonian (gj) is not 
exact. However, it is a very good approximation in the 
regime where t c ^> g/isAiJjv- 9 

To visualize the effects of J and ABm we draw the two- 
electron spin states using a Bloch sphere representation 
in Fig. [321 The effect of J in this representation is to 
rotate the Bloch vector about the z-axis of the Bloch 
sphere. An initially prepared \H) spin state will rotate 
into a ||t) s P m state in a time te=ttTi/ J(s). This is a 
SWAP opera tion. Le aving J "on" for half of this time 
per forms a V SWAP operation. 

V SWAP combined with single-spin rotations can be 
used to create arbitrary quantum gates. In fact, this 
two-spin operation allows universal quantum comput- 
ing by itself, when the logical qubit is encoded in three 
spins (jDiVincenzo et aEEooo). If an inhomogenous ef- 
fective magnetic field is present, encoding a qubit in just 



For J 3> g/igAB^f, \S) and |Tb) are good eigenstates. The 
mixing term g/igABpf only has an effect when J < gfig ABjv ■ 
Thus, if \S) is almost equal to 5(1,1) for J < g/igABjy , the use 
of Hamiltonian 1441) is valid. If t c 2> g/igABpi, the condition 
J < gfigABjy is only met far away from the avoided crossing. 
Here, J is of order t 2 / |e|, and the condition can thus be rewrit- 
ten as t 2 / |e| < gfigABj^f. Since t c S> g/J,gAB]\f, it follows 
that t c / \e\ <C 1. Therefore, the weight of the 5(0,2)-component 

in \S), given by ss t c / (2W e 2 + t§), is indeed negligible in this 
regime. 
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FIG. 48 (Color in online edition) Coherent two-electron spin state rotations, (a) Pulse sequence, (b) Singlet state probability, 
Ps, measured as a function of the pulse time, te and detuning, e. (c) Horizontal cuts through the data in (b) show clear 
oscillations in Ps . (d) By increasin g the tunnel coupling, a fast V SWAP operation time of ~180 ps is achieved. Data 
reproduced from iPetta et al\ l|2005bl ). 



two spins is sufficient for creating any quantum gate us- 
ing just the exchange interaction ( Levvl . I2002T ) . In this 
system, the qubit basis states are the singlet and the 
To triplet state. Note that a Loss-DiVincenzo v SWAP 
operation corresponds to a single-qubit rotation in the 
singlet-triplet basis. 

A SWAP operation has been implemented using fast 
control of J. The pulse sequence is illustrated in Fig. l4~5b . 
The system is prepared at positive detuning in 5(0,2). 
The singlet is then spatially separated by making the 
detuning more negative; this is at first done fast with 
respect to the hyperfine mixing rate T 2 * to avoid mixing 
with the T + state. Once beyond the 5-T + degeneracy, the 
detuning is lowered further, but now slowly with respect 
to T 2 *. This prepares the system in the ground state 
of the hyperfine fields, here defined |fi)- This state is 
an eigenstate of the nuclear fields and is insensitive to 
hyperfine fluctuations. 

To perform coherent two-electron spin rotations a pulse 
is applied to the system which increases the energy split- 
ting J between 5 and To. This drives a z-axis rotation in 
the Bloch sphere representation by an angle 9. The ro- 
tation is then turned off by again lowering the detuning. 
A spin state projection measurement is performed by re- 
versing the initialization process, thereby mapping |TJ.) — * 
|5(1, 1)) and ||T) — * |Tb(l, 1)). Spin-to-charge conversion 
is then used to determine the spin state. 



Figure l48l b) shows the measured singlet state proba- 
bility as a function of the rotation pulse time te and e 
during the rotation pulse. Ps shows clear oscillations as a 
function of both e and te- The period of the oscillations 
agrees well with a theoretical calculation obtained using 
a calibration of J(e) from the 5(1, 1)-T+(1, 1) resonance 
condition. Horizontal cuts through the data are shown in 
Fig. l487 c). By increasing t c and hence J, a fast V SWAP 
operation time of 180 ps is obtained (see Fig. l4"8Td)). 

Fast control of the J can be harnessed to implement 
a singlet-triplet spin echo pulse sequence. As shown in 
Fig. [47l hyperfine fields lead to fast dephasing of the spin 
singlet state. In the Bloch sphere representation, 
drives a random x-axis rotation. Since B n is a fluctu- 
ating quantity, this rotation rate will vary from one ex- 
perimental run to the next. However, since the nuclear 
spin dynamics are much slower than the electron spin dy- 
namics, the hyperfine dephasing can be reversed using a 
spin-echo pulse sequence. 

The spin-echo pulse sequence is illustrated in 
Fig. l49T a). The singlet state, 5(0,2) is prepared at pos- 
itive detuning. The detuning is decreased quickly with 
respect to Bn but slowly compared to t c , creating a (1,1) 
singlet state. Each spin evolves in the presence of the hy- 
perfine fields during the separation time rg, which in the 
Bloch sphere representation corresponds to an x-axis ro- 
tation. An exchange pulse of angle n is applied to the 



4G 




x E (ns) t s +t S ' (us) 

FIG. 49 (Color in online edition) Error correction, (a) Singlet-triplet spin echo pulse sequence, (b) Singlet state probability, 
Ps, measured as a function of pulse time, te and detuning, e. Clear singlet state recover ies are obser v ed for ir, 3-7T, 5it exchange 
pulses, (c) The singlet state recovery persists out to T2=1.2 fjs. Data reproduced from iPetta et al\ l|2005rJ ). 



system, which rotates the Bloch vector about the z-axis of 
the Bloch sphere. Exchange is turned off and the spins 
evolve for a time ts' ■ During this time, the hypcrfmc 
fields rotate the Bloch vector back towards 5(1,1), refo- 
cusing the spin singlet state. 

Figure [49] (b) shows Ps as a function of e and te in 
the spin echo pulse sequence. Ps shows clear oscillations 
as a function of te- For tt, 3 7r, and 5 tt pulses clear 
singlet state recoveries are observed. To determine the 
coherence time we set ts=ts> and vary the total sep- 
aration time t to t=Ts+Ts' ■ Figure l49Tc) shows Ps as a 
function of ts~ts> for increasing t to t ■ A singlet state re- 
covery is observed for t to t exceeding 1 microsecond. A 
fit to the singlet state decay using an exponential form 
leads to a best fit T-x=\.1 microseconds. Remarkably, 
this spin echo pulse sequence extends the coherence time 
by a factor of 100. Experiments are currently being per- 
formed to determine the physical origin of the 1.2 mi- 
crosecond decay. Possible sources of the decay are nu- 
clear spin_jvohvtfon_(see Secti on IVII.Bp , charge dephas- 
ing ( Hu and Das Sarmal 120061) . and decay to T + (l,l). 



X. PERSPECTIVES 

This review has described the spin physics of few- 
electron quantum dots. Much of this work can be 
evaluated within the context of spin-based classical or 
quantum information processing. In this context, the 



state-of-the-art can be best summarized by making a 
comparison with the first five "DiVincenzo criteria" 
( DiVincenzd . l2000h . applied to the Loss-DiVincenzo pro- 
posal for_en£Ojiing_a_Jo^^ in a single electron 
spin ( Loss and DiVincenzd . Il998h . 

1. Have a scalable physical system with well-defined 
qubits. Electron spins are certainly well-defined qubits. 
The Zeeman energy difference between the qubit states 
can be made much larger than the thermal energy. The 
states can be measured using transport spectroscopy (see 
section III). Concerning scalability it is difficult to make 
predictions. In principle, circuits of solid state devices 
are scalable, but evidently many practical problems will 
have to be surmounted. 

2. Be able to initialize to a simple fiducial state such as 
pOOO...). By waiting until relaxation takes place at low 
temperature and in high magnetic field, the many-qubit 
ground state will be occupied with probability close to 1. 
Another option is to use the energy difference between 
the states or the different coupling to the reservoir to 
induce spin-selective tunneling from the reservoir onto 
the dot. 

3. Have long coherence times. The T2-coherence time 
has not been determined extensively, but already, a lower 
bound of ~ 1/j.s has been established at 100 mT. How 
quantum coherence scales with the size of the system is 
an interesting open question. The coherence times of 
qubits can be prolonged by error correction, one of the 
holy grails in this field. This can be done effectively only 
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when many manipulations are allowed before decoher- 
ence takes place. The rule of thumb is that the coherence 
time should be at least 10 4 times longer than the time 
for a typical one- or two qubit operation. 

4. Have a universal set of quantum gates. The Loss 
and DiVincenzo proposal provides two gates which to- 
gether allow for universal quantum computing. Single- 
qubit rotations have been implemented by ESR, with a 
fidelity of ~ 75%, and a duration of 25 ns for a ir/2 rota- 
tion. The two-qubit gate is based on the two-spin SWAP 
operation, which has been demonstrated as well, com- 
bined with single-spin rotations. The SWAP has already 
been operated at sub-ns levels (180 ps for a %/SWAP 
gate), although the fidelity is yet unclear. These are only 
the first experimental results and further improvements 
are expected. 

5. Permit high quantum efficiency, qubit-specific mea- 
surements. The procedure of spin to charge conversion 
and measuring the charge is a highly efficient measure- 
ment of the qubit state. It allows for a single-shot read- 
out measurement with demonstrated fidelities already ex- 
ceeding 90%. An optimization of experimental parame- 
ters can certainly increase this to > 99%. We note also 
that the QPC charge meter is a fairly simple device that 
can be integrated easily in quantum dot circuits. 

We see that qubits defined by single electron spins in 
quantum dots largely satisfy the DiVincenzo criteria. As 
an alternative, it is also possible to encode the logical 
qubit in a combination of spins. For instance, when the 
logical qubit is encoded in three spins instead of a sin- 
gle spin, the exchange interaction b y itself is sufficient 
for u niversal quantum computation ( DiVincenzo et all , 
2000). Adding a difference in Zeeman energy between the 
two dots reduces the number of spins per logical qubit to 
two (iLevvi . 12002ft . Coherent operations on this so-called 
singlet-triplet qubit have already been experimentally 
demonstrated (see Section HX.B[) . These two- and three- 
electron qubit encodings eliminate the need for the tech- 
nologically challenging single-spin rotations. Many more 
variations for encoding qubits in several electron spins 



drawbacks (Bvrd and Lidar, 2002: 


Hanson and Burkard 


2007; Kvriakidis and Pennev. 20051: iMeier et all 2003 


Tavlor et al. 


. l2005t IWu and Lidar. 2002a b). In the end 



the best implementation for a given system will depend 
on many factors that are hard to oversee at this stage. 

In the near future, the natural continuation of the re- 
cent work will be to combine the various components 
(readout, ESR and exchange gate) in a single experiment. 
This may allow for new experiments exploring quantum 
coherence in the solid state, for instance involving non- 
local entanglement and testing Bell's inequalities. As an- 
other example, the precise role of quantum measurements 
may be investigated in this system as well. 

On a longer timescale, the main challenges are scalabil- 
ity and coherence. Scalability is mostly a practical issue. 
The coherence challenge provides a number of very inter- 
esting open questions. The coherence time is currently 



limited by the randomness in the nuclear spin system. 
If this randomness is suppressed the coherence time will 
become longer. Polarization of the nuclei turns out not 
to be very efficient, except for polarizations > 99.9%. As 
an alternative, the nuclear spins coul d be put and kept in 
a particular, known quantum state dGiedke et~al\ . l2006t 
iKlauser efafl . l2006tlStepanenko et a/.l . l2006f) . ~ 

It is yet unknown if nuclear spins can indeed be con- 
trolled up to a high level of accuracy. A completely dif- 
ferent approach would be using a different material. The 
isotopes of the III-V semiconductors all have a non-zero 
nuclear spin. In contrast, the group IV semiconductors 
do have isotopes with zero nuclear spin. If spin qubits 
are realized in a material that is isotopically purified to 
for instance 28 Si or 12 C only, the hyperfine interaction is 
completely absent. 

We believe that the techniques and physics described 
in this review will prove valuable regardless of the type 
of quantum dot that is used to confine the electrons. The 
unprecedented level of control over single electron spins 
will enable exploration of new regimes and pave the way 
for tests of simple quantum protocols in the solid state. 
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APPENDIX A: Sign of the ground state spin and the 
nuclear fields in GaAs 

In this Appendix we derive the sign of the ground state 
of electron and nuclear spins in GaAs and the sign and 
magnitude of the effective magnetic field felt by electrons 
due to thermal and dynamical nuclear polarization. 
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1. Sign of the spin ground states 

We define the spin to be 'up' if it is oriented in the di- 
rection of the externally applied magnetic field B z along 
the z-axis. In other words, an electron with spin S is 
spin-up if the quantum number for the ^-component of 
the spin, S z , is positive. The magnetic moments asso- 
ciated with the electron spin S and the nuclear spin J 
are 



Ms 
A 4 / 



-gs 



N 

2 m, 



lei 

9i^t 



-gsHBS z 



(Al) 
(A2) 



where (1b and fi^ are the Bohr magneton (57.9 /zeV/T) 
and the nuclear magneton (3.15 neV/T), respectively 
(note that in our notation, the spin angular momentum 
along z is given by TiS z ). The difference in the sign of the 
magnetic moments is due to the difference in the polarity 
of the electron and proton charge. The Zeeman energy 
is given by Ez = —fi'B. Since both free electrons and 
protons have a positive g-factor, the spins in the ground 
states of a free electron (spin-down) and a proton (spin- 
up) are anti-parallel to eachother. 

The nuclear ^-factors of the isotopes in GaAs are 
all positive: gi( 69 Ga) =+1.344, gi( 71 Ga) =+1.708 and 
gi( 75 As) =+0.960. The electron (/-factor in GaAs is neg- 
ative (<7s=-0.44). Hence, both the nuclei and the elec- 
trons in the ground state in GaAs have their spin aligned 
parallel to the external field, i.e. they are spin-up. 



2. Sign and magnitude of the thermal nuclear field 

The two Ga-isotopes, 69 Ga (60.11% abundance) and 
71 Ga (38.89% abundance), and 75 As all have nuclear 
spin 3/2. We can calculate the thermal average of the 
spin <I> of each isotope using the Maxwell-Boltzmann 
distribution. For example, at 10 T and 20 mK, < 
i>B9 Go =+0.30, <J>7i (7n =+ 0.38 and <I>75 A 



-0.22. 

Then, following Paget et al. ( Paget et ali\l9Ta i we ap- 
proximate the effective field, generated by the polariza- 
tion of isotope a through the hyperfine contact interac- 
tion, by 



Bn,<x = b N (a) (I) a 



(A3) 



with 6jv( b9 Ga)=-0.91 T, M 71 Ga)=-0. 78 T and 
b N ( 7S A s)=-1.8A T (formula 2.17-19 of iPaget et all 
()1977D ). Since (I) a is always positive in thermal equi- 
librium, we derive from equation IA3I that the thermal 
nuclear field acts against the applied field. 



3. Sign of the dynamic nuclear field 

A nuclear polarization can build up dynamically via 
flip-flop processes, where an electron and a nucleus flip 
their spin simultaneously. Because of the large energy 
mismatch between nuclear and electron Zeeman energy, 



a flip-flop process where an electron spin is excited is very 
unlikely, since the required energy is not available in the 
system [AEz,nuci "Cfcs? 1 <C AEz, e i) ■ Therefore, we only 
consider the flip-flop processes where the electron flips 
its spin from down to up (AS z —+l), thereby releasing 
the Zeeman energy. This brings the nucleus to a differ- 
ent spin state with AI Z —-1. Many of these processes can 
dynamically build up a considerable polarization, whose 
sign is opposite to that of the thermal nuclear field. This 
has already been observed in th e ESR experiments on 
2DEGs (see e.g. iDobers et all (|l988h ). where the ex- 
cited electron spin relaxes via a flip-flop process. The 
external field at which the ESR field is resonant shifts 
to lower values after many of these processes, indicating 
that indeed this nuclear field adds to the external field. 



References 

Abalmassov, V. A., and F. Marquardt, 2004, Phys. Rev. B 
70, 75313. 

Abragam, A., 1961, The Principles of Nuclear Magnetism, 
Clarendon. 

Abragam, A., and B. Bleaney, 1986, Electron Paramagnetic 
Resonance of Transition Ions (Dover Publications). 

Aguado, R., and L. P. Kouwenhoven, 2000, Phys. Rev. Lett. 
84, 1986. 

Akera, H., 1999, Phys. Rev. B 60, 10683. 

Amasha, S., K. MacLean, I. Radu, D. M. Zumbuhl, M. A. 

Kastner, M. P. Hanson, and A. C. Gossard, 2006, cond- 

mat/0607110 . 

Ashcroft, N. W., and N. D. Mermin, 1974, Solid State Physics 

(Saunders, New York). 
Ashoori, R. C, H. L. Stormer, J. S. Weiner, L. N. Pfeiffer, 

S. J. Pearton, K. W. Baldwin, and K. W. West, 1992, Phys. 

Rev. Lett. 68, 3088. 
Astafiev, O., Y. Pashkin, T. Yamamoto, Y. Nakamura, and 

J. S. Tsai, 2004, Phys. Rev. B 69, 180507. 
Atature, M., J. Dreiser, A. Badolato, A. Hogele, K. Karrai, 

and A. Imamoglu, 2006, Science 312, 551. 
Augustine, M. P., and E. L. Hahn, 1997, The Journal of 

Chemical Physics 107, 3324. 
Awschalom, D. D., and M. Flatte, 2007, Nature Physics 3, 

153. 

Berezovsky, J., M. H. Mikkelsen, O. Gywat, N. G. Stoltz, 
L. A. Coldren, and D. D. Awschalom, 2006, Science 314, 
1916. 

Willems van Beveren, L. H, R. Hanson, I. T. Vink, F. H. L. 

Koppens, L. P. Kouwenhoven, and L. M. K. 3, 2005, New 

Journal of Physics 7, 182. 
Bir, G. L., A. G. Aronov, and G. E. Pikus, 1976, Sov. Phys. 

JETP 42, 705. 

Bjork, M. T., A. Fuhrer, A. E. Hansen, M. W. Larsson, L. E. 

Froberg, and L. Samuelson, 2005, Phys. Rev. B 72, 201307. 
Bjork, M. T., C. Thelander, A. E. Hansen, L. E. Jensen, 

M. W. Larsson, L. R. Wallenberg, and L. Samuelson, 2004, 

Nano Letters 4, 1621. 
Bockelmann, U., 1994, Phys. Rev. B 50, 17271. 
Borhani, M., V. N. Golovach, and D. Loss, 2006, Phys. Rev. 

B 73, 155311. 

Bracker, A. S., E. A. Stinaff, D. Gammon, M. E. Ware, J. G. 
Tischler, A. Shabaev, A. L. Efros, D. Park, D. Gershoni, 



49 



V. L. Korenev, and I. A. Merkulov, 2005, Phys. Rev. Lett. 
94, 47402. 

Braun, P. F., X. Marie, L. Lombez, B. Urbaszek, T. Amand, 
P. Renucci, V. K. Kalevich, K. V. Kavokin, O. Krebs, 
P. Voisin, and Y. Masumoto, 2005, Phys. Rev. Lett. 94, 
116601. 

Braunstein, S. L., C. M. Caves, R. Jozsa, N. Linden, 
S. Popescu, and R. Schack, 1999, Phys. Rev. Lett. 83, 1054. 

Bulaev, D. V., and D. Loss, 2005, Phys. Rev. B 71, 205324. 

Burkard, G., 2001, Quantum Computation and Communi- 
cation using Electron Spins in Quantum Dots and Wires, 
Ph.D. thesis, University of Basel. 

Burkard, G., D. Loss, and D. P. DiVincenzo, 1999, Phys. Rev. 
B 59, 2070. 

Bychkov, Y. A., and E. I. Rashba, 1984, JETP Lett. 39, 78. 
Byrd, M. S., and D. A. Lidar, 2002, Phys. Rev. Lett. 89, 
47901. 

Cheng, J. L., M. W. Wu, and C. Lii, 2004, Phys. Rev. B 69, 
115318. 

Ciorga, M., M. Pioro-Ladriere, P. Zawadzki, P. Hawrylak, and 
A. S. Sachrajda, 2002, Appl. Phys. Lett. 80, 2177. 

Ciorga, M., A. S. Sachrajda, P. Hawrylak, C. Gould, P. Za- 
wadzki, S. Jullian, Y. Feng, and Z. Wasilewski, 2000, Phys. 
Rev. B 61, R16315. 

Climente, J. I., A. Bertoni, G. Goldoni, M. Rontani, and 
E. Molinari, 2007, Phys. Rev. B 75, 081303(R). 

Cobden, D. H., M. Bockrath, P. L. McEuen, A. G. Rinzler, 
and R. E. Smalley, 1998, Phys. Rev. Lett. 81, 681. 

Coish, W. A., and D. Loss, 2004, Phys. Rev. B 70, 195340. 

Coish, W. A., and D. Loss, 2005, Phys. Rev. B 72, 125337. 

Cooper, J., C. G. Smith, D. A. Ritchie, E. H. Linfield, Y. Jin, 
and H. Launois, 2000, Physica E: Low-dimensional Systems 
and Nanostructures 6(1-4), 457. 

Cota, E., R. Aguado, and G. Platero, 2005, Phys. Rev. Lett. 
94, 107202. 

Cronenwett, S. M., T. H. Oosterkamp, and L. P. Kouwen- 

hoven, 1998, Science 281, 540. 
Debald, S., and C. Emary, 2005, Phys. Rev. Lett. 94, 226803. 
Dekker, C, 1999, Physics Today 52, 22. 

von Delft, J., and D. C. Ralph, 2001, Physics Reports 345, 
61. 

Deng, C, and X. Hu, 2006, Phys. Rev. B 73, 241303. 
DiCarlo, L., H. Lynch, A. C. Johnson, L. I. Childress, 

K. Crockett, C. M. Marcus, M. P. Hanson, and A. C. Gos- 

sard, 2004, Phys. Rev. Lett. 92, 226801. 
Dickmann, S., and P. Hawrylak, 2003, JETP Letters 77(1), 

30. 

DiVincenzo, D. P., 2000, Fortschritte Der Physik-Progress of 

Physics 48(9-11), 771. 
DiVincenzo, D. P., D. Bacon, J. Kempe, G. Burkard, and 

K. B. Whaley, 2000, Nature 408, 339. 
Dixon, D. C, K. R. Wald, P. L. McEuen, and M. R. Melloch, 

1997, Phys. Rev. B 56, 4743. 
Dobers, M., K. Klitzing, and G. Weimann, 1988, Phys. Rev. 

B 38, 5453. 

Dobrowolska, M., H. D. Drew, J. Furdyna, T. Ichiguchi, 
A. Witowski, and P. A. Wolff, 1982, Phys. Rev. Lett. 49, 
845. 

Dresselhaus, G., 1955, Phys. Rev. 100, 580. 

Duncan, D. S., D. Goldhaber-Gordon, R. M. Westervelt, 

K. D. Maranowski, and A. C. Gossard, 2000, Appl. Phys. 

Lett. 77, 2183. 

Dutt, M. V. G., J. Cheng, B. Li, X. Xu, X. Li, P. R. Berman, 
D. G. Steel, A. S. Bracker, D. Gammon, S. E. Economou, 



et al, 2005, Phys. Rev. Lett. 94, 227403. 
Dyakonov, M. I., and V. Y. Kachorovskii, 1986, Sov. Phys. 

Semicond. 20, 110. 
Dyakonov, M. L, and V. I. Perel, 1972, Sov. Phys. Solid State 

13(12), 3023. 

Dzhioev, R. I., V. L. Korenev, I. A. Merkulov, B. P. Za- 

kharchenya, D. Gammon, A. L. Efros, and D. S. Katzer, 

2002, Phys. Rev. Lett. 88, 256801. 
Ellenberger, C, T. Ihn, C. Yannouleas, U. Landman, K. En- 

sslin, D. Driscoll, and A. C. Gossard, 2006, Phys. Rev. Lett. 

96, 126806. 
Elliott, R. J., 1954, Phys. Rev. 96, 266. 

Elzerman, J. M., R. Hanson, L. H. Willems van Beveren, 
L. M. K. Vandersypen, and L. P. Kouwenhoven, 2004a, 
Appl. Phys. Lett. 84, 4617. 

Elzerman, J. M., R. Hanson, L. H. W. van Beveren, 
B. Witkamp, L. M. K. Vandersypen, and L. P. Kouwen- 
hoven, 2004b, Nature 430, 431. 

Elzerman, J. M., R. Hanson, J. S. Greidanus, L. H. 
Willems van Beveren, S. De Franceschi, L. M. K. Vander- 
sypen, S. Tarucha, and L. P. Kouwenhoven, 2003, Phys. 
Rev. B 67, 161308. 

Engel, H.-A., V. N. Golovach, D. Loss, L. M. K. Vandersypen, 
J. M. Elzerman, R. Hanson, and L. P. Kouwenhoven, 2004, 
Phys. Rev. Lett. 93, 106804. 

Engel, H.-A., and D. Loss, 2001, Phys. Rev. Lett. 86, 4648. 

Engel, H.-A., and D. Loss, 2002, Phys. Rev. B 65, 195321. 

Engel, H.-A., and D. Loss, 2005, Science 309, 586. 

Erlingsson, S. I., and Y. V. Nazarov, 2002, Phys. Rev. B 66, 
155327. 

Erlingsson, S. I., and Y. V. Nazarov, 2004, Phys. Rev. B 70, 
205327. 

Erlingsson, S. I., Y. V. Nazarov, and V. I. Falko, 2001, Phys. 

Rev. B 64, 195306. 
Falko, V., B. Altshuler, and O. Tsyplyatyev, 2005, Phys. Rev. 

Lett. 95, 76603. 
Field, M., C. G. Smith, M. Pepper, D. A. Ritchie, J. E. F. 

Frost, G. A. C. Jones, and D. G. Hasko, 1993, Phys. Rev. 

Lett. 70, 1311. 

Florescu, M., and P. Hawrylak, 2006, Phys. Rev. B 73, 45304. 
Fransson, J., and M. Rasander, 2006, Phys. Rev. B 73, 
205333. 

Freeman, R., 1997, Spin choreography: basic steps in high 

resolution NMR (Spektrum, Oxford). 
Friesen, M., C. Tahan, R. Joynt, and M. A. Eriksson, 2004, 

Phys. Rev. Lett. 92, 37901. 
Fujisawa, T., D. G. Austing, Y. Tokura, Y. Hirayama, and 

S. Tarucha, 2002a, Nature 419, 278. 
Fujisawa, T., D. G. Austing, Y. Tokura, Y. Hirayama, and 

S. Tarucha, 2002b, Phys. Rev. Lett. 88, 236802. 
Fujisawa, T., T. Hayashi, Y. Hirayama, H. D. Cheong, and 

Y. H. Jeong, 2004, Appl. Phys. Lett. 84(13). 
Fujisawa, T., T. Hayashi, and S. Sasaki, 2006a, Rep. Prog. 

Phys 69, 759. 

Fujisawa, T., T. Hayashi, R. Tomita, and Y. Hirayama, 2006b, 

Science 312, 1634. 
Fujisawa, T., T. H. Oosterkamp, W. G. van der Wiel, B. W. 

Broer, R. Aguado, S. Tarucha, and L. P. Kouwenhoven, 

1998, Science 282, 932. 
Fujisawa, T., Y. Tokura, and Y. Hirayama, 2001a, Physica B 

298(1-4), 573. 

Fujisawa, T., Y. Tokura, and Y. Hirayama, 2001b, Phys. Rev. 
B 63, 81304. 

Giedke, G., J. M. Taylor, D. DAlessandro, M. D. Lukin, and 



■50 



A. Imamoglu, 2006, Phys. Rev. A 74, 32316. 
Goldhaber-Gordon, D., H. Shtrikman, D. Mahalu, D. Abusch- 
Magder, U. Meirav, and M. A. Kastner, 1998, Nature 391, 
156. 

Golovach, V. N., M. Borhani, and D. Loss, 2006, Phys. Rev. 

Lett. 74, 165319. 
Golovach, V. N., A. Khaetskii, and D. Loss, 2004, Phys. Rev. 

Lett. 93, 016601. 
Golovach, V. N., A. Khaetskii, and D. Loss, 2007, cond- 

mat/0703427 . 

Greentree, A. D., A. R. Hamilton, L. Hollenberg, and R. G. 
Clark, 2005, Phys. Rev. B 71, 113310. 

Greilich, A., R. Oulton, E. A. Zhukov, I. A. Yugova, D. R. 
Yakovlev, M. Bayer, A. Shabaev, L. E. Al, I. A. Merkulov, 
V. Stavarache, D. Renter, and A. Wieck, 2006a, Phys. Rev. 
Lett. 96, 227401. 

Greilich, A., D. R. Yakovlev, A. Shabaev, A. L. Efros, I. A. 
Yugova, R. Oulton, V. Stavarache, D. Reuter, A. Wieck, 
and M. Bayer, 2006b, Mode Locking of Electron Spin Co- 
herences in Singly Charged Quantum Dots. 

Gueron, S., M. M. Deshmukh, E. B. Myers, and D. C. Ralph, 
1999, Phys. Rev. Lett. 83, 4148. 

Gunther, L., and B. Barbara, 1994, Quantum Tunneling of 
Magnetization (Kluwer, Dordrecht, The Netherlands). 

Gupta, J. A., R. Knobel, N. Samarth, and D. D. Awschalom, 
2001, Ultrafast Manipulation of Electron Spin Coherence. 

Hanson, R., 2005, Electron spins in semicon- 
ductor quantum dots, Ph.D. thesis, Delft Uni- 
versity of Technology, online available from 
http : / /www.library. tudelft .nl /dissertations / . 

Hanson, R., L. H. Willems van Beveren, I. T. Vink, J. M. Elz- 
erman, W. J. M. Naber, F. H. L. Koppens, L. P. Kouwen- 
hoven, and L. M. K. Vandersypen, 2005, Phys. Rev. Lett. 
94, 196802. 

Hanson, R., and G. Burkard, 2007, Phys. Rev. Lett. 98, 
050502. 

Hanson, R., O. Gywat, and D. D. Awschalom, 2006, Phys. 
Rev. B 74, 161203(R). 

Hanson, R., L. M. K. Vandersypen, L. H. Willems van Bev- 
eren, J. M. Elzerman, I. T. Vink, and L. P. Kouwenhoven, 
2004a, Phys. Rev. B 70, 241304(R). 

Hanson, R., I. T. Vink, D. P. DiVincenzo, L. M. K. Vander- 
sypen, J. M. Elzerman, L. H. Willems van Beveren, and 
L. P. Kouwenhoven, 2004b, Proceedings of the 39th Ren- 
contres de Moriond, see also cond-mat/0407793 . 

Hanson, R., B. Witkamp, L. M. K. Vandersypen, L. H. 
Willems van Beveren, J. M. Elzerman, and L. P. Kouwen- 
hoven, 2003, Phys. Rev. Lett. 91, 196802. 

Hayashi, T., T. Fujisawa, and Y. Hirayama, 2003, Physica 
Status Solidi(b) 238(2), 262. 

Heiss, D., M. Kroutvar, J. J. Finley, and G. Abstreiter, 2005, 
Solid State Communications 135, 591. 

Hermann, C, and C. Weisbuch, 1977, Phys. Rev. B 15(2), 
823. 

Herzog, B., and E. L. Hahn, 1956, Phys. Rev. 103, 148. 
Hitachi, K., M. Yamamoto, and S. Tarucha, 2006, Phys. Rev. 
B 74, 161301. 

Hofmann, F., T. Heinzel, D. A. Wharam, J. P. Kotthaus, 
G. B6hm, W. Klein, G. Trankle, and G. Weimann, 1995, 
Phys. Rev. B 51, 13872. 

van Houten, H., C. W. J. Beenakker, and A. A. M. Staring, 
1992, in Single Charge Tunneling, edited by H. Grabert 
and M. H. Devoret, NATO ASI Series B294 (Plenum, New 
York, 1992), sec also cond-mat/0508454. 



Hu, X., and S. Das Sarma, 2006, Phys. Rev. Lett. 96, 100501. 

Hiittel, A. K., S. Ludwig, H. Lorenz, K. Eberl, and J. P. 
Kotthaus, 2005, Phys. Rev. B 72, 81310. 

Hiittel, A. K., J. Weber, A. W. Holleitner, D. Weinmann, 
K. Eberl, and R. H. Blick, 2004, Phys. Rev. B 69, 73302. 

Imamoglu, A., D. D. Awschalom, G. Burkard, D. P. DiVin- 
cenzo, D. Loss, M. Sherwin, and A. Small, 1999, Phys. Rev. 
Lett. 83, 4204. 

Ionicioiu, R., and A. E. Popescu, 2005, New Journal of 
Physics 7, 120. 

Jelezko, F., T. Gaebel, I. Popa, A. Gruber, and J. Wrachtrup, 

2004, Phys. Rev. Lett. 92, 076401. 
Jiang, H. W., and E. Yablonovitch, 2001, Phys. Rev. B 64, 

041307. 

Johnson, A. C, C. M. Marcus, M. P. Hanson, and A. C. 

Gossard, 2005a, Phys. Rev. B 71, 115333. 
Johnson, A. C, J. R. Petta, C. M. Marcus, M. P. Hanson, 

and A. C. Gossard, 2005b, Phys. Rev. B 72, 165308. 
Johnson, A. C, J. R. Petta, J. M. Taylor, A. Yacoby, M. D. 

Lukin, C. M. Marcus, M. P. Hanson, and A. C. Gossard, 

2005c, Nature 435, 925. 
Jouravlev, O. N., and Y. V. Nazarov, 2006, Phys. Rev. Lett. 

96, 176804. 

Jung, S. W., T. Fujisawa, Y. H. Jeong, and Y. Hirayama, 

2004, Appl. Phys. Lett. 85, 768. 
Kane, B. E., 1998, Nature 393, 133. 

Kato, Y., R. C. Myers, D. C. Driscoll, A. C. Gossard, J. Levy, 
and D. D. Awschalom, 2003a, Gigahertz Electron Spin Ma- 
nipulation Using Voltage-Controlled g- Tensor Modulation. 

Kato, Y., R. C. Myers, A. C. Gossard, and D. D. Awschalom, 
2003b, Nature 427, 50. 

Khaetskii, A. V., D. Loss, and L. Glazman, 2002, Phys. Rev. 
Lett. 88, 186802. 

Khaetskii, A. V., D. Loss, and L. Glazman, 2003, Phys. Rev. 
B 67, 195329. 

Khaetskii, A. V., and Y. V. Nazarov, 2000, Phys. Rev. B 61, 
12639. 

Khaetskii, A. V., and Y. V. Nazarov, 2001, Phys. Rev. B 64, 
125316. 

Klauder, J. R., and P. W. Anderson, 1962, Phys. Rev. 125, 
912. 

Klauser, D., W. A. Coish, and D. Loss, 2006, Phys. Rev. B 
73, 205302. 

Klein, D. L., P. L. McEuen, J. E. B. Katari, R. Ross, and 

A. P. Alivisatos, 1996, Applied Physical Letters 68, 2574. 
Kogan, A., S. Amasha, D. Goldhaber-Gordon, G. Granger, 

M. A. Kastner, and H. Shtrikman, 2004, Phys. Rev. Lett. 
93, 166602. 

Konemann, J., R. J. Haug, D. K. Maude, V. I. Falko, and 

B. L. Altshuler, 2005, Phys. Rev. Lett. 94, 226404. 
Koppens, F. H. L., C. Buizert, K.-J. Tielrooij, I. T. Vink, 

K. C. Nowack, T. Meunier, L. P. Kouwenhoven, and 
L. M. K. Vandersypen, 2006, Nature 442, 766. 

Koppens, F. H. L., J. A. Folk, J. M. Elzerman, R. Han- 
son, L. H. W. van Beveren, I. T. Vink, H. P. Tranitz, 
W. Wegscheider, L. P. Kouwenhoven, and L. M. K. Van- 
dersypen, 2005, Science 309, 1346. 

Korkusihski, M., P. Hawrylak, M. Ciorga, M. Pioro-Ladriere, 
and A. S. Sachrajda, 2004, Phys. Rev. Lett. 93, 206806. 

Kouwenhoven, L. P., D. G. Austing, and S. Tarucha, 2001, 
Rep. Prog. Phys. 64, 701. 

Kouwenhoven, L. P., C. M. Marcus, P. L. McEuen, 
S. Tarucha, R. M. Westervelt, and N. S. Wingreen, 1997, 
Kluwer, Series E 345, 16, in Mesoscopic Electron Trans- 



51 



port, edited by L. L. Sohn, L. P. Kouwenhoven and G. 
Schon. 

Krenner, H. J., E. C. Clark, T. Nakaoka, M. Bichler, 

C. Scheurer, G. Abstreiter, and J. J. Finley, 2006, Phys. 

Rev. Lett. 97, 76403. 
Kroutvar, M., Y. Ducommun, D. Heiss, M. Bichler, D. Schuh, 

G. Abstreiter, and J. J. Finley, 2004, Nature 432, 81. 
Kupidura, D., M. C. Rogge, M. Reinwald, W. Wegscheider, 

and R. J. Haug, 2006, Phys. Rev. Lett. 96, 46802. 
Kyriakidis, J., and S. J. Penney, 2005, Phys. Rev. B 71, 

125332. 

Kyriakidis, J., M. Pioro-Ladriere, M. Ciorga, A. S. Sachrajda, 
and P. Hawrylak, 2002, Phys. Rev. B 66, 35320. 

Levy, J., 2002, Phys. Rev. Lett. 89, 147902. 

Lindemann, S., T. Ihn, T. Heinzel, W. Zwerger, K. Ensslin, 
K. Maranowski, and A. C. Gossard, 2002, Phys. Rev. B 66, 
195314. 

Loss, D., and D. P. DiVincenzo, 1998, Phys. Rev. A 57, 120. 
Lu, W., Z. Q. Ji, L. Pfeiffer, K. W. West, and A. J. Rimberg, 

2003, Nature 423, 422. 
Marquardt, F., and V. A. Abalmassov, 2005, Phys. Rev. B 

71, 165325. 

McEuen, P. L., 2000, Physics World 13(6), 31. 
Meier, F., J. Levy, and D. Loss, 2003, Phys. Rev. Lett. 90, 
47901. 

Meier, F., and B. P. Zakharchenya, 1984, Optical orientation. 
Modern problems in condensed matter sciences (North- 
Holland). 

Merkulov, I. A., A. L. Efros, and M. Rosen, 2002, Phys. Rev. 
B 65, 205309. 

Meunier, T., I. T. Vink, L. H. Willems van Beveren, F. H. L. 
Koppens, H. P. Tranitz, W. Wegscheider, L. P. Kouwen- 
hoven, and L. M. K. Vandersypen, 2006, Phys. Rev. B 74, 
195303. 

Meunier, T., I. T. Vink, L. H. Willems van Beveren, K. J. 
Tielrooij, R. Hanson, F. H. L. Koppens, H. P. Tranitz, 
W. Wegscheider, L. P. Kouwenhoven, and L. Vandersypen, 
2007, Phys. Rev. Lett. 98, 126601. 

Naber, W. J. M., T. Fujisawa, H. W. Liu, and W. G. van der 
Wiel, 2006, Phys. Rev. Lett. 96, 136807. 

Nielsen, M. A., and I. L. Chuang, 2000, Quantum Compu- 
tation and Quantum Information (Cambridge Univ. Press, 
Cambridge). 

Ohno, H., 1998, Science 281, 951. 

Onac, E., F. Balestro, L. H. Willems van Beveren, U. Hart- 

mann, Y. V. Nazarov, and L. P. Kouwenhoven, 2006, Phys. 

Rev. Lett. 96, 176601. 
Ono, K., D. G. Austing, Y. Tokura, and S. Tarucha, 2002, 

Science 297, 1313. 
Ono, K., and S. Tarucha, 2004, Phys. Rev. Lett. 92, 256803. 
Paaske, J., A. Rosch, P. Woffle, N. Mason, C. M. Marcus, 

and J. Nygard, 2006, Nature Physics 2, 460. 
Paget, D., 1982, Phys. Rev. B 25, 4444. 
Paget, D., G. Lampel, B. Sapoval, and V. I. Safarov, 1977, 

Phys. Rev. B 15, 5780. 
Park, J., A. N. Pasupathy, J. I. Goldsmith, C. Chang, 

Y. Yaish, J. R. Petta, M. Rinkoski, J. P. Sethna, H. D. 

Abruna, and P. L. McEuen, 2002, Nature 417, 722. 
Petta, J. R., A. Johnson, A. Yacoby, C. Marcus, M. Hanson, 

and A. Gossard, 2005a, Phys. Rev. B 72, 161301. 
Petta, J. R., A. C. Johnson, C. M. Marcus, M. P. Hanson, 

and A. C. Gossard, 2004, Phys. Rev. Lett. 93, 186802. 
Petta, J. R., A. C. Johnson, J. M. Taylor, E. A. Laird, A. Ya- 
coby, M. D. Lukin, C. M. Marcus, M. P. Hanson, and A. C. 



Gossard, 2005b, Science 309, 2180. 
Petta, J. R., and D. C. Ralph, 2001, Phys. Rev. Lett. 87, 
266801. 

Pfeffer, P., 1999, Phys. Rev. B 59, 15902. 

Pioro-Ladriere, M., M. R. Abolfath, P. Zawadzki, J. Lapointe, 

S. A. Studenikin, A. S. Sachrajda, and P. Hawrylak, 2005, 

Phys. Rev. B 72, 125307. 
Platero, G., and R. Aguado, 2004, Phys. Rep. 395, 1. 
Poole, C, 1983, Electron Spin Resonance, 2nd ed. (Wiley, 

New York) . 

Potok, R. M., J. A. Folk, C. M. Marcus, V. Umansky, M. Han- 
son, and A. C. Gossard, 2003, Phys. Rev. Lett. 91, 16802. 

Ralph, D. C, C. T. Black, and M. Tinkham, 1995, Phys. Rev. 
Lett. 74, 3241. 

Rashba, E. I., 1960, Sov. Phys. Solid State 2, 1109. 

Recher, P., E. Sukhorukov, and D. Loss, 2000, Phys. Rev. 
Lett. 85, 1962. 

Reimann, S. M., and M. Manninen, 2002, Rev. Mod. Phys. 
74, 1283. 

Rugar, D., R. Budakian, H. J. Mamin, and B. W. Chui, 2004, 

Nature 430, 329. 
Rushforth, A. W., C. G. Smith, M. D. Godfrey, H. E. Beere, 

D. A. Ritchie, and M. Pepper, 2004, Phys. Rev. B 69, 

113309. 

Salis, G., Y. Kato, K. Ensslin, D. C. Driscoll, A. C. Gossard, 

and D. D. Awschalom, 2001, Nature 414, 619. 
San- Jose, P., G. Zarand, A. Shnirman, and G. Schon, 2006, 

Phys. Rev. Lett. 97, 76803. 
Sanchez, R., E. Cota, R. Aguado, and G. Platero, 2006, Phys. 

Rev. B 74, 35326. 
Sanchez, R., G. Platero, R. Aguado, and E. Cota, 2006, phys. 

stat. sol.(b) 1, 5. 
Sasaki, S., D. G. Austing, and S. Tarucha, 1998, Physica B 

256, 157. 

Sasaki, S., T. Fujisawa, T. Hayashi, and Y. Hirayama, 2006, 

Phys. Rev. Lett. 95, 056803. 
Schleser, R., E. Ruh, T. Ihn, K. Ensslin, D. C. Driscoll, and 

A. C. Gossard, 2004, Appl. Phys. Lett. 85, 2005. 
Schliemann, J., A. V. Khaetskii, and D. Loss, 2002, Phys. 

Rev. B 66, 245303. 
Schoelkopf, R. J., P. Wahlgren, A. A. Kozhevnikov, P. Dels- 

ing, and D. E. Prober, 1998, Science 280, 1238. 
Schulte, M., J. G. S. Lok, G. Denninger, and W. Dietsche, 

2005, Phys. Rev. Lett. 94, 137601. 
Semenov, Y. G., and K. W. Kim, 2004, Phys. Rev. Lett. 92, 

26601. 

Shenvi, N., and K. Rogerio de Sousa, 2005, Phys. Rev. B 71, 
224411. 

Shulman, R. G., B. J. Wyluda, and H. J. Hrostowski, 1958, 

Phys. Rev. 109, 808. 
Sih, V., Y. Kato, and D. D. Awschalom, 2005, Physics World 

18, 33. 

Simovic, B., P. Studerus, S. Gustavsson, R. Leturcq, K. En- 
sslin, R. Schuhmann, J. Forrer, and A. Schweiger, 2006, 
Review of Scientific Instruments 77, 064702. 

Slichter, C. P., 1990, Principles of Magnetic Resonance, 3rd 
ed. (Springer-Verlag, Berlin). 

Snelling, M. J., G. P. Flinn, A. S. Plaut, R. T. Harley, A. C. 
Tropper, R. Eccleston, and C. C. Phillips, 1991, Phys. Rev. 
B 44, 11345. 

de Sousa, R., 2006, cond-mat/0610716 . 

de Sousa, R., and S. Das Sarma, 2003a, Phys. Rev. B 67, 
33301. 

de Sousa, R., and S. Das Sarma, 2003b, Phys. Rev. B 67, 



■52 



33301. 

de Sousa, R., and S. Das Sarma, 2003c, Phys. Rev. B 68, 
155330. 

Sprinzak, D., Y. Ji, M. Heiblum, D. Mahalu, and H. Shtrik- 

man, 2002, Phys. Rev. Lett. 88, 176805. 
Stano, P., and J. Fabian, 2005, Phys. Rev. B 72, 155410. 
Stano, P., and J. Fabian, 2006, Phys. Rev. Lett. 96, 186602. 
Stepanenko, D., G. Burkard, G. Giedke, and A. Imamoglu, 

2006, Phys. Rev. Lett. 96, 136401. 
Su, B., V. J. Goldman, and J. E. Cunningham, 1992, Phys. 

Rev. B 46, 7644. 
Tarucha, S., Y. Kitamura, T. Kodera, and K. Ono, 2006, 

phys. stat. sol. (b) 243, 3673. 
Taylor, J. M., H. A. Engel, W. Dur, A. Yacoby, C. M. Marcus, 

P. Zoller, and M. D. Lukin, 2005, Nature Physics 1. 
Teraoka, S., A. Numata, S. Amaha, K. Ono, and S. Tarucha, 

2004, Physica E: Low-dimensional Systems and Nanostruc- 

tures 21, 928. 

Tokura, Y., W. G. Van der Wiel, T. Obata, and S. Tarucha, 
2006, Phys. Rev. Lett. 96, 047202. 

Vandersypen, L. M. K., and I. L. Chuang, 2004, Rev. Mod. 
Phys. 76, 1037. 

Vandersypen, L. M. K., J. M. Elzerman, R. N. Schouten, L. H. 
Willems van Beveren, R. Hanson, and L. P. Kouwenhoven, 
2004, Appl. Phys. Lett. 85, 4394. 

Vandersypen, L. M. K., R. Hanson, L. H. Willems van Bev- 
eren, J. M. Elzerman, J. S. Greidanus, S. De Franceschi, 
and L. P. Kouwenhoven, 2002, in Quantum Computing 
and Quantum Bits in Mesoscopic Systems, Kluwer Aca- 
demic/Plenum Publishers (New York, 2003) See also quant- 
ph/0207059. 

Viola, L., E. Knill, and S. Lloyd, 1999, Phys. Rev. Lett. 82, 
2417. 

Viola, L., and S. Lloyd, 1998, Phys. Rev. A 58, 2733. 
Wald, K. R., L. P. Kouwenhoven, P. L. McEuen, N. C. van der 

Vaart, and C. T. Foxon, 1994, Phys. Rev. Lett. 73, 1011. 
Wang, Y. Y., and M. W. Wu, 2006, Phys. Rev. B 74, 165312. 
van Wees, B. J., H. van Houten, C. W. J. Beenakker, J. G. 

Williamson, L. P. Kouwenhoven, D. van der Marel, and 



C. T. Foxon, 1988, Phys. Rev. Lett. 60, 848. 
Weinmann, D., W. Hausler, and B. Kramer, 1995, Phys. Rev. 

Lett. 74, 984. 

Weis, J., R. J. Haug, K. Klitzing, and K. Ploog, 1993, Phys. 

Rev. Lett. 71, 4019. 
Wharam, D. A., T. J. Thornton, R. Newbury, M. Pepper, 

H. Ahmed, J. E. F. Frost, D. G. Hasko, D. C. Peacock, 

D. A. Ritchie, and G. A. C. Jones, 1988, J. Phys. C 
21(L209-214), 26. 

Van der Wiel, W. G., S. De Franceschi, J. M. Elzerman, T. Fu- 
jisawa, S. Tarucha, and L. P. Kouwenhoven, 2003, Rev. 
Mod. Phys. 75, 1. 

Witzel, W. M., and S. Rogerio de Sousa, 2005, Phys. Rev. B 
72, 161306. 

Wolf, S. A., D. D. Awschalom, R. A. Buhrman, J. M. 
Daughton, S. von Molnar, M. L. Roukes, C. A. Y., and 
D. M. Treger, 2001, Science 294, 1488. 

Woods, L. M., T. L. Reinecke, and Y. Lyanda-Geller, 2002, 
Phys. Rev. B 66, 161318. 

Wrinklcr, R., 2003, Spin-orbit coupling effects in two- 
dimensional electron and hole systems (Springer- Verlag, 
Berlin) . 

Wu, L. A., and D. A. Lidar, 2002a, Phys. Rev. A 65, 42318. 
Wu, L. A., and D. A. Lidar, 2002b, Phys. Rev. A 66, 62314. 
Xiao, M., I. Martin, E. Yablonovitch, and H. W. Jiang, 2004, 

Nature 430, 435. 
Yafet, Y., 1963, Solid State Physics 14, 1. 
Yao, W., R. B. Liu, and L. J. Sham, 2005, Phys. Rev. B 74, 

195301. 

Zhang, L.-X., J. P. Leburton, R. Hanson, and L. P. Kouwen- 
hoven, 2004, Appl. Phys. Lett. 85, 2628. 

Zumbiihl, D. M., C. M. Marcus, M. P. Hanson, and A. C. 
Gossard, 2004, Phys. Rev. Lett. 93, 256801. 

Zumbiihl, D. M., J. B. Miller, C. M. Marcus, K. Campman, 
and A. C. Gossard, 2002, Phys. Rev. Lett. 89, 276803. 

Zutic, I., J. Fabian, and S. Das Sarma, 2004, Rev. Mod. Phys. 
76, 323. 



